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ABSTRACT 


This  study  examines  seismic  discrimination  between  underground  nuclear  explo¬ 
sions  and  earthquakes  in  the  Southwestern  United  States.  A  thorough  review 
of  theoretical  and  applied  research  on  this  problem,  especially  as  it  relates 
to  that  region,  is  presented  first,  followed  by  a  presentation  of  the  seismic 
discrimination  parameters  computed  for  a  suite  of  events  in  the  Southwest  and 
a  series  of  experiments  with  multi-dimensional  discrimination. 


Review  of  past  work  finds  some  theoretical  support  for  successful  dis¬ 
crimination  based  upon  several  distinct  measurements  on  seismic  recordings. 
Although  negative  first  motion  and  presence  of  pP  both  indicate  a  natural 
earthquake  source,  they  cannot  be  found  confidently  in  many  cases.  The  rela¬ 
tive  level  of  shear-wave  and  Love-wave  phases  should  be  a  good  classification 
parameter.  The  long-period  P  spectrum  of  explosions  should  be  diminished  due 
to  the  surface  interaction  and  consequent  cancellation  by  pP.  The  surface- 

wave  magnitude  M  should  be  only  0.2  higher  for  small,  shallow  earthquakes  com- 
s 

pared  to  equivalent  explosions,  if  the  source  time  functions  are  identical. 

Observed  separation  is  routinely  greater  than  0.2  though,  and  this  fact 

requires  incoherent  rupture  for  earthquakes  and  perhaps  some  overshoot  in 

explosion  source  time  functions.  High  frequency  spectral  decay  of  earthquakes 
-3  -2 

should  be  to  while  that  of  explosions  is  to  .  However,  other  effects,  such 
as  attenuation,  may  overwhelm  this  difference.  For  earthquakes  that  have  com¬ 
plex  ruptures  and  appreciable  focal  depth,  the  complexity  of  the  P-wave  record¬ 
ing  should  be  higher,  but  not  greatly  so.  While  Rayleigh-wave  spectral  shapes 
for  earthquakes  and  explosions  at  shallow  depths  should  not  be  much  different, 
deeper  earthquakes  could  be  identified  by  higher  modes  and  a  relative  decrease 
of  shorter  period  amplitude. 


The  region  under  study  extends  from  California  to  the  southern  Rocky 
Mountains  and  from  roughly  40°N  to  the  Gulf  of  California.  A  region  of  high 
seismicity,  it  is  also  a  complex  region  encompassing  several  tectonic  pro¬ 
vinces,  including  the  Basin-Range  Province  where  the  Nevada  Test  Site  is 
located.  Source  mechanisms  are  diverse  over  the  entire  region  and  even 
within  small  subregions.  Focal  depths  are  almost  all  <  15  km  and  few  depths 
can  be  determined  reliably  from  teleseismic  recordings.  Twelve  earthquakes. 


mostly  selected  because  of  low  Mg  estimates  from  other  studies,  and  eleven 
explosions  comprised  the  sources  for  this  study.  Seven  teleseismic  LRSM  sta¬ 
tions  supplied  the  seismic  data. 

In  general,  seismic  discriminant  parameters  obtained  from  the  recordings 
reflected  the  theoretical  expectations  of  earthquake-explosion  differences. 
Path-station  effects  were  large  for  every  parameter,  especially  for  short- 
period  data.  Love-wave  magnitudes  were  a  superior  discriminant  when  plotted 
versus  m^.  Shear  waves,  if  measurable,  were  also  excellent.  The  common  Mg- 
m^  plot  shows  three  earthquakes  with  anomalously  low  Mg:  the  Denver  earth¬ 
quake,  a  Benham  aftershock,  and  a  Baja  California  earthquake. 

All  the  discrimination  parameters  (10)  were  applied  in  various  combina¬ 
tions  through  a  stepwise  linear  discriminant  program,  which  treated  the  twenty- 
three  events  as  a  training  set.  Using  all  parameters,  the  lowest  a  posteriori 
probability  of  correct  classification  was  .987.  Multi-dimensional  discrimin¬ 
ation  using  network  averages  for  parameters  was  superior  to  using  single-sta¬ 
tion  parameters  or  using  linear  combinations  of  individual  stations'  discrimi¬ 
nant  functions. 
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INTRODUCTION 


Past  studies  in  seismic  discrimination  between  earthquakes  and  under¬ 
ground  nuclear  explosions  have  generally  been  regional  in  scope,  in  which 
attempts  at  discrimination  involved  neighboring  events  on  the  globe.  Although 
they  should  not  be  undervalued,  these  studies  often  failed  to  discern  or  to 
convey  important  regional  properties  of  earthquake  mechanisms  and  earth  struc¬ 
ture  that  may  debilitate  or  enhance  identification  of  events,  either  earth¬ 
quakes  or  explosions,  from  other  areas.  Also,  although  most  studies  have  been 
characterized  as  regional,  in  many  cases  the  region  chosen  is  so  large  that 
distinct  subregions  with  peculiar  properties  exist  which  significantly  affect 
discrimination  parameters.  Evernden  and  Filson  (1971)  and  Marshall  and  Basham 
(1972),  in  the  case  of  M^,  and  Douglas  et  al.  (1973),  in  the  case  of  P-wave 
complexity,  provide  examples  of  attempts  to  overcome  regional  effects  and  to 
present  discrimination  parameters  in  a  unified  sense.  Identification  and 
estimation  of  path  and  receiver  effects  on  discrimination  parameters  are 
essential  if  a  sound  physical  basis  for  distinguishing  natural  and  artificial 
seismic  events  is  to  be  established  not  only  in  regions  where  they  have  been 
studied,  but  also  in  any  unstudied  region  where  seismic  events  may  occur  in 
isolated  instances. 

Magnitude  determination  shows  the  importance  of  removing  or  suppressing 
path-receiver  effects.  Here,  investigators  have  distinguished  between  shallow 
earthquakes  and  explosions  whose  theoretical  separation  on  the  Mg  scale 
(given  equal  m^)  amounts  to  only  a  half  order  of  magnitude  based  on  infini¬ 
tesimal  source  theory  (Gilbert,  1973).  They  accomplish  this  by  using  signal 


Evernden,  J.  F.,  and  J.  Filson  (1971).  Regional  dependence  of  surface-wave 
versus  body-wave  magnitudes.  J,  Geophys.  Res.,  27,  3303-3308. 

Marshall,  P.  D. ,  and  P.  W.  Basham  (1972).  Discrimination  between  earthquakes 
and  underground  explosions  employing  an  improved  Ms  scale.  Geophys.  J. , 
28,  431-458. 

Douglas,  A.,  P.  D.  Marshall,  P.  G.  Gibbs,  J.  B.  Young,  and  C.  Blarney  (1973). 
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amplitude  estimates  known  typically  to  range  over  nearly  two  orders  of  magni¬ 
tude  for  a  well-recorded  seismic  event.  A  casual  approach  to  magnitude  deter¬ 
mination  results  in  explosion  points  ranging  over  almost  two  units  of  magni¬ 
tude  in  Mg  for  a  given  m^,  and  vice-versa,  when  data  from  several  test  sites 
are  plotted  together  (Liebermann  and  Pomeroy,  1969),  Although  some  of  the 
observed  variation  is  undoubtedly  source-related,  the  Rayleigh-wave  spectral 
ratios  of  von  Seggern  and  Lambert  (1970)  and  the  short-period  spectral  ratios 
of  Lacoss  (1969)  provide  additional  examples  of  severe  regional  effects  on 
discrimination. 

While  most  previous  studies  produced  valuable  detailed  coverage  of  one 
aspect  of  discrimination,  they  failed  to  study  several  facets  of  the  data  in 
one  region,  an  omission  that  impairs  evaluation  of  various  discrimination 
parameters  and  limits  efforts  to  improve  discrimination  using  parameters 
graphically  and  statistically,  either  alone  or  in  combination.  Only  a  few 
authors,  together  comprising  only  a  small  portion  of  the  relevant  literature, 
have  used  two  or  more  distinct  discriminants  to  separate  natural  and  artificial 
seismic  events.  Among  them  are  Booker  and  Mitronovas  (1964),  Lacoss  (1969), 
Press  et  al. ,  (1963),  Evernden  (1969),  Rasmussen  and  Lande  (1968), 
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Anglin  (1971),  and  Dahlman  et  al.  (1974).  Even  though  high  reliability  of  the 
Mg  versus  m^  discriminant  has  been  establisned  (Evernden  et  al.,  1971;  Ericsson, 
1970;  Pasechnik  et  al.,  1970;  Weichert  and  Basham,  1974;  Evernden,  1975),  the 
recognition  of  certain  earthquakes  generating  anomalously  low  amplitude  surface 
waves  (Landers,  1972;  Der,  1973)  and  the  frequent  inability  to  detect  surface  waves 
from  explosions  compels  us  to  consider  multiple-variable  discrimination.  The 
importance  of  this  mode  of  surveillance  is  reinforced  by  the  fact  that  anoma¬ 
lous  earthquakes  in  those  studies  which  were  marginally  identified  or  even  uni¬ 
dentified  by  Mg  versus  m^  were  quite  clearly  identified  by  means  of  other  dis¬ 
criminants.  (Whether  anomalous  Mg  earthquakes  in  all  regions  of  the  globe  can 
be  so  easily  disposed  of  remains  to  be  seen.)  Another  reason  to  consider  sev¬ 
eral  discriminants  rather  than  one  is  the  possibility  of  clandestine  nuclear 
blasts  under  test-ban  treaty  controls.  Evaders  could  either  disguise  the  test 
to  appea*  like  an  earthquake  or  hide  the  test  in  the  coda  signals  of  a  large 
earthquake.  However,  such  schemes  may  not  be  successful  if  a  sophisticated, 
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multiple-discriminant  effort  is  employed  (Evernden,  1976). 

The  purpose  of  this  study  is  to  overcome  these  deficiencies  through  a 
comprehensive,  coherent  study  of  explosions  and  earthquakes  in  a  selected 
region.  This  study  attempts  to  identify  and  to  quantify  the  numerous  factors 
at  the  source,  along  the  travel  path,  and  at  the  receivers  which  affect  dis¬ 
crimination  for  a  particular  region  of  the  globe.  Discrimination  parameters 
are  to  be  related  to  the  nature  of  the  sources  and  the  structure  of  the  region 
so  that  why  they  do  or  do  not  work  can  be  understood.  Although  this  is  a 
regional  study,  the  authors  intend  to  emerge  with  enough  knowledge  that  regional 
effects  could  be  satisfactorily  removed,  if  necessary,  from  discrimination  para¬ 
meters  estimated  here.  The  authors  hope  to  refine  discrimination  techniques 
enough  in  the  selected  area  to  permit  meaningful  application  to  other  regions. 

The  choice  of  an  initial  region — the  Western  United  States  (WUS) — for  an 
integrated  study  was  a  natural  one.  Specifically,  the  region  is  composed 
largely  of  the  Central  Cordillera  and  the  Pacific  Ranges,  plus  a  small  north¬ 
western  section  of  the  Southern  Cordillera  (in  all,  roughly  30°-40°N  and  104°- 
124°W) .  This  area  is  not,  in  any  sense,  a  small  or  homogeneous  region,  but  it 
was  deliberately  chosen  to  be  large  enough  to  include  earthquakes  and  explosions 
of  differing  charactistics  and  to  present  some  structural  contrasts.  At  the 
same  time  it  is  small  enough  to  allow  analysis  using  a  common  seismic  network 
and  to  avoid  drastic  path  effects  which  might  otherwise  hinder  development  of 
conclusive  or  significant  results. 

Other  reasons  for  selecting  the  Western  United  States  are: 

(1)  Data,  in  the  form  of  LRSM  (Long  Range  Seismic  Measurement)  film 
and  tape  recordings  and  film  chips  for  many  WWSSN  sites,  are 
readily  available  and  easily  analyzed. 

(2)  This  region  contains  several  underground  nuclear  test  sites  and 
sites  of  peaceful  nuclear  explosions.  The  large  explosions  fired 

Evernden,  J.  (1976).  Study  of  seismological  evasion.  Part  III.  Evaluation 
of  evasion  possibilities  using  codas  of  large  earthquakes.  Bull.  Seism. 
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in  the  Western  United  States  are  well  documented,  and  many  have 
‘  •  been  studied  in  detail.  Parameters  such  as  detonation  depth, 

yield,  and  medium  are  known,  and  the  seismic  signals  have  been 
used  for  previous  discrimination  studies  and  for  investigations 
into  earth  structure. 

>  3)  The  Western  United  States  is  a  region  of  pronounced  seismicity, 

and  many  of  the  larger  earthquakes  in  this  region  are  well  docu¬ 
mented.  Focal  mechanism  solutions,  fault  length,  source  spectral 
shape,  as  well  as  other  characteristics,  have  been  presented  for 
literally  hundreds  of  events. 

4)  Because  some  earthquakes  apparently  are  only  marginally  discrimi¬ 
nated  by  Mg-m^  diagrams  (Basham,  1969;  Peppin  and  McEvilly,  1974; 
Peppin,  1976;  Savino  et  al.,  1971)  or  by  short-period  spectral 
ratios  (Bakun  and  Johnson,  1970),  real  discrimination  problems 
exist  for  events  in  the  Western  United  States.  The  one  element 
that  most  confuses  identification  of  explosions  at  the  Nevada 
Test  Site  is  the  large  amount  of  tectonic  strain  release  accom¬ 
panying  many  explosions,  especially  those  of  larger  yield, 

(Toksoz  and  Kehrer,  1972a),  which  causes  LQ/LR  ratios  for  those 
explosions  to  equal  or  surpass  those  of  earthquakes. 
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5) 


The  geology  and  geophysics  of  the  Western  United  States  (from  the 
surface  down  through  the  upper  mantle)  is  known  as  well  as,  or 
better  than,  that  of  any  other  global  region.  Such  knowledge 
will  aid  considerably  in  determining  source  and  path  effects  on 
seismic  waves  because  one  should  be  able  to  relate  peculiarities 
in  seimic  observations  to  known  or  hypothetical  regional  struc¬ 
tural  elements. 

6)  Because  of  the  abundance  of  background  and  supporting  material  for 
the  Western  United  States,  the  objectives  of  this  study  can  be 
more  easily  accomplished  than  if  an  isolated  and  uncharted  region 
were  chosen. 

7)  Lastly,  since  so  much  discrimination  literature  on  the  Western 
United  States  has  already  been  published,  this  study  provides  an 
appropriate  forum  from  which  to  review  it,  to  synthesize  it,  and 
to  integrate  it  with  any  new  results  from  this  study. 

This  study  is  divided  into  three  major  sections.  First,  published  the¬ 
ories  and  results  pertinent  to  discrimination,  especially  those  pertaining 
to  events  in  the  Western  United  States,  are  reviewed.  Here,  the  physical 
bases  of  discrimination  will  be  emphasized  and  illustrated.  Then,  by  provid¬ 
ing  a  description  of  the  present  structure  of  the  Western  United  States,  a 
specific  framework  for  this  study  will  be  established.  A  review  of  the  pre¬ 
sent  tectonics  will  follow,  including  active  crustal  movements  and  recent 
seismicity.  After  this  background  information,  the  choice  of  data  for  the 
regional  discrimination  study  will  be  listed  and  described.  One  part  of  the 
study  will  be  signal  analysis,  leading  to  descriptions  and  estimations  of 
source,  path,  and  receiver  effects  and  their  impact  on  discrimination.  A 
second  part  of  the  study  will  involve  use  of  discrimination  parameters  in  a 
stepwise  multiple-discrimination  scheme.  Finally,  the  last  section  will  pre¬ 
sent  conclusions  and  recommendations  for  further  study. 


GENERAL  REVIEW  OF  SEISMIC  DISCRIMINATION 


Some  of  the  earliest  discussions  (e.g.  VESIAC,  1962)  of  discrimination 
between  earthquakes  and  underground  nuclear  explosions  focused  on  the  dis¬ 
tinctness  of  the  explosion  source  from  the  earthquake  source  in  terms  of 
symmetry,  compactness,  or  absence  of  rotational  components.  Ironically,  how¬ 
ever,  in  these  same  volumes  are  reports  of  numerous  observations  that  serve 
only  to  confound  discrimination  based  upon  this  simple  concept.  Out  of  this 
confusion  emerged  a  vigorous  program  of  research  that  investigated  the  nature 
and  effects  of  the  two  sources.  Most  of  the  earlier  discrimination  work 
involved  explosions  and  earthquakes  in  the  Western  United  States  because  of 
the  numerous  detonations  at  the  Nevada  Test  Site,  the  high  seismicity  of  this 
region,  the  proximity  of  existing  seismograph  stations,  and  the  addition  of 
many  new  recording  sites  in  the  form  of  LRSM  (Long  Range  Seismic  Measurement) 
stations  and  VELA  Uniform  observatories.  The  results  of  the  empirical,  as  well 
as  the  theoretical,  studies  are  reviewed  below.  For  the  purpose  of  this  review, 
and  the  study  in  general,  a  discn  inant  is  defined  as  any  measurable  informa¬ 
tion  on  seismograms  that  can  be  used,  on  a  statistical  or  other  basis,  to  deter¬ 
mine  the  source  of  a  seismic  event. 

Location,  Including  Depth 

This  discriminant  is  unique  because  it  ignores  the  physical  nature  of  the 
source.  When  an  event  has  been  located  in  an  aseismic  area,  reference  to  num- 
ber-versus-magnitude  recurrence  curves,  which  are  presumably  based  on  natur¬ 
al  events,  may  suggest  the  events  are  artificial.  Since  a  possibility  exists 
for  large  events  in  even  aseismic  areas,  this  information  does  not  provide 
conclusive  evidence  of  an  explosion,  but  it  is  sufficient  to  separate  out  the 
event  for  further  study.  The  boundaries  of  aseismic  areas  are  often  sharply 
defined,  a  situation  creating  the  need  for  high-quality  seismic  network  loca¬ 
tions  for  monitoring  any  test-ban  agreement.  Neglecting  epicenter  bias,  such 
quality  locations  could  be  attained  (95%  confidence  ellipses  of  <500  square 
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km)  for  events  down  to  m^  ~  4  over  almost  the  entire  globe  (Evernden,  1971) 
with  as  few  as  twenty-five  good  stations.  However,  the  epicenter  bias  from 
lateral  inhomogeneity  of  the  earth  can  severely  undermine  the  utility  of  these 
confidence  regions,  unless  relative  location  using  master  events  is  employed 
(Evernden,  1971),  though  serious  bias  is  probably  limited  to  descending  plate 
areas  (Davies  and  McKenzie,  1969).  Note  that  because  of  the  relatively  high 
seismicity  of  the  Western  United  States,  especially  in  the  areas  surrounding 
the  Nevada  Test  Site,  discrimination  based  upon  areal  location  would  be  impos¬ 
sible,  even  if  confidence  ellipses  were  precise  and  small. 

In  contrast,  events  in  seismic  areas  can  often  be  confirmed  as  earthquakes 
if  the  site  itself  would  preclude  an  underground  nuclear  explosion  (for  exam¬ 
ple,  beneath  deep  water  at  sea  or  within  the  boundaries  of  a  nation  known  not 
to  be  capable  of  testing  nuclear  explosions).  Again,  high-quality  location  is 
the  key  to  making  these  decisions. 

The  depth  variable  in  location  can  also  be  considered  a  discriminant.  An 
event  located  at  a  depth  that  would  be  too  deep  for  a  nuclear  detonation,  say 
>5  km,  could  be  classified  as  a  natural  earthquake.  However,  limitations  of 
a  teleseismic  network  in  thic  regard  are  severe;  for  with  teleseismic  location 
methods,  depth  cannot  he  shown  by  consideration  of  location  error  ellipses 
alone  to  be  beyond  very  shallow,  except  for  events  located  by  travel-time 
inversion  deeper  than  roughly  30  to  40  km  (Evernden,  1969).  Again,  bias  in 
depth  equally  serious  to  that  of  areal  location  can  arise  due  to  lateral  inho¬ 
mogeneities  in  the  earth  such  as  dipping  plates,  ocean  ridges,  and  continental 
areas  with  anomalously  thick  crust;  a  well-known  example  is  the  location  of 
LONG  SHOT  roughly  60  km  deep  (Lambert  et  al.,  1969).  Davies  and  McKenzie 
(1969)  suggested  that  travel  time  will  be  badly  biased  relative  to  the 
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standard  curve  in  dipping  plate  areas.  In  the  Western  United  States,  almost 
all  earthquakes  are  shallow  C <15  km);  therefore,  a  teleseismic  network  alone 
could  not  discriminate  between  earthquakes  and  explosions  by  routine  location 


methods  unless  pP  could  be  definitely  identified. 


Using  the  spectra  of  seismic  events  provides  another  approach  to  depth 
determination.  However,  results  obtained  from  this  method  are  questionable 
because  no  general  statistical  measure  exists  comparable  to  the  error  ellipses 
used  with  travel  time.  In  fact,  this  method  has  produced  little  fruitful 
results  in  the  field  of  discrimination  because  all  events  where  spectral  inver¬ 
sion  was  applied  were  already  known  to  be  either  earthquakes  or  explosions. 

At  this  time,  depth  information  derived  from  spectra  can  be  considered  only  of 
diagnostic  value;  that  is,  only  a  very  weak  discriminant.  Indeed,  this  method 
can  admit  gross  errors  in  estimates  of  event  depth  because  of  the  many  other 
factors  that  influence  the  shape  of  seismic  spectra,  including  source  time 
function,  lateral  inhomogeneity,  crustal  and  upper  mantle  layers,  and  attenu¬ 
ation.  Also,  the  possibility  exists  that  the  inversion  can  converge  on  a 
completely  wrong  solution.  Some  workers  resist  these  criticisms  and  contend 
that  solution  for  source  depth  of  earthquakes  is  feasible  with  spectra  of 
either  surface  waves  (Keilis-Borok  and  Yanovskaya,  1962;  Tsai  and  Aki,  1970; 
Tsai  and  Aki,  1971;  Weidner  and  Aki,  1973;  Canitez  and  Toksoz,  1971)  or  body 
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waves  (Cuha  and  Stauder,  1970;  Langston,  1976).  In  the  case  of  surface  waves, 

Canitez  and  Toksoz  claim  a  depth  precision  of  a  few  kilometers  can  be  attained 

using  ratios  of  LQ/LR  spectra,  which  is  a  means  of  cancelling  some  of  the  pro- 

* 

pagation  effects  on  spectra.  However,  in  response  Masse  et  al.  (1973)  pointed 
out  that  large  errors  in  depth  are  still  possible  using  the  LQ/LR  spectral 
ratios,  especially  if  the  earth  structure  is  not  accurately  known. 

Carrying  this  argument  further,  note  that  many  investigators  have  esti¬ 
mated  detonation  depths  of  explosions  from  short-period  P  waves  or  their  spectra 
(Cohen,  1970;  Kulhanek,  1971;  Douglas  et  al.,  1972a;  Frasier,  1972;  Manchee  and 
Hasegawa,  1973)  and  that  the  explosion  depth  estimates  seem  fairly  accurate  with 
this  method  because  of  the  good  agreement  between  different  authors  using  dif¬ 
fering  stations  for  common  events.  However,  these  authors  knew  ,a  priori  that 
the  events  were  explosions,  and  they  are  not  able  to  suggest  that  the  events  can 
be  classified  as  explosions  on  the  basis  of  shallow-focus  estimates.  Discrim¬ 
ination  is  not  really  possible  with  this  technique  because  earthquakes  also 
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can  be  near  the  surface  and  because  earthquakes  at  any  depth  may,  through  multi¬ 
ple-path  effects  or  s  complicated  source  function,  show  spectral  shaping  similar 
to  that  which  was  used  to  estimate  explosion  depth  (Cohen  et  al.,  1972). 

First  Motion 

Because  only  natural  earthquakes,  not  explosions,  are  ideally  able  to  gen¬ 
erate  initial  dilatational  motions,  the  P-wave  first  motions  on  seismic  record¬ 
ings  of  an  event  should  determine  whether  the  event  was  an  explosion  or  an 
earthquake.  However,  several  difficulties  exist  with  this  simple  approach. 
First,  the  spheroidal  shape  and  layering  of  the  earth  may  combine  with  tele- 
seismic  networks  such  that  the  two  quadrants  of  a  double-couple  mechanism  which 
generate  the  initial  dilatational  motions  appear  not  to  be  well-sampled  when 
raypaths  are  traced  back  to  the  source.  Carpenter  (1964)  estimated  that  between 
10%  and  25%  of  all  earthquakes  have  no  teleseismic  dilatations,  depending  upon 
assumptions  made  for  network  distribution,  earth  structure,  and  geometry  of  the 
focal  mechanisms.  Evernden  (1969)  pointed  out  some  concrete  examples  of  where 
first  motion  criteria  would  fail.  Second,  dilatational  motion  has  been  observed 
for  explosions  worldwide  (Enescu  et  al.,  1973)  where  obviously  ideal  source  con¬ 
ditions  are  not  satisfied  and  where  some  nonuniform  forces  are  at  work.  This 
phenomenon  is  usually  associated  with  the  generation  of  S  waves,  another  process 
that  cannot  occur  with  an  ideal  explosion.  The  explanation  may  be  one  or  a  com¬ 
bination  of  the  following:  Taylor  instabilities  in  the  inelastic  region  result¬ 
ing  from  nonuniformities  in  the  medium,  simultaneous  release  of  tectonic  strain 
within  the  inelastic  zone,  or  triggering  actual  earthquakes  near  the  explosion. 
These  hypotheses  will  be  discussed  in  a  later  section  dealing  with  shear-wave 
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discriminants.  Finally,  even  if  ideal  conditions  at  the  source  are  assumed, 
seismic  background  noise  hampers  the  use  of  this  discriminant  by  masking  or 
confusing  first  motion  direction  creating  a  rather  high  effective  thres¬ 
hold  for  its  application,  roughly  m^  ~  5  1/2  according  to  Evernden  (1969). 
Still,  in  spite  of  drawbacks,  this  discriminant  is  often  irrefutable  in 
identifying  earthquakes  either  when  many  clear  dilatations  are  recorded  or 
in  identifying  explosions  when  complete  azimuthal  coverage  of  an  event  at 
regional  distances  fails  to  show  any  dilatations. 

Shear  Waves 

Assuming  that  an  ideal  explosion  is  purely  a  sphere  of  outward  compression 
and  that  an  earthquake  is  essentially  a  shear  fracture,  identification  of 
explosions  should  be  based  upon  the  absence  of  shear  waves  on  seismic  record¬ 
ings.  However,  numerous  observations  of  explosion  signals  show  that  they 
almost  always  include  S-wave  motion  (e.g.,  Kisslinger  et  al.,  1961;  Willis  et 
al.,  1963;  Press  et  al. ,  1963;  Booker  and  Mitronovas,  1964;  Geyer  and  Martner, 
1969;  Hattori,  1972;  von  Seggern,  1972;  von  Seggern,  1973;  and  Blandford  and 
Clark,  1974).  Vertically-polarized  S  can  be  explained  by  conversion  of  P  at 
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interfaces  in  the  earth,  but  the  studies  under  review  here  also  usually 
report  horizontally-polarized  S  waves  which,  because  of  their  arrival  time, 
are  likely  to  be  generated  only  at  the  source.  Observations  of  LQ  surface 
waves  from  explosions  (e.g.,  Press  et  al . ,  1963;  Lambert  et  al.,  1972;  Toksoz 
and  Kehrer,  1972b)  are  also  prevalent  and,  based  upon  travel  time,  further 
attest  to  a  shear-wave  generating  mechanism  at  or  near  the  explosion. 

Any  one  of  several  hypotheses  can  explain  horizontally-polarized  shear 
waves  from  an  explosion.  Among  these  are  mode  conversion,  radiation  asymmetry, 
radial  cracking,  relaxation  of  non-uniform  prestress,  stable  sliding  on  exist¬ 
ing  faults,  and  earthquake  triggering.  Of  these,  the  first  three  could  be 
operative  in  a  stress-free  or  hydrostatically  stressed  environment,  while  the 
others  require  some  degree  of  shear  stress  in  the  immediate  area  of  the  explo¬ 
sion.  Evidence  from  various  explosion  sites  indicates  that  no  single  explana¬ 
tion  is  usually  correct  and  that  explanations  may  vary  with  magnitude,  depth, 
and  exact  location  of  the  detonation.  (In  addition  to  our  discussion  of  these 
effects,  see  Press  and  Archambeau  (1962),  Aki  and  Tsai  (1972),  or  ToksBz  et 
al.  (1971)  who  may  provide  additional  insight.)  In  regard  to  mode  conversion 
of  P  or  SV  at  irregular  interfaces  along  the  travel  path,  Oliver  et  al. 

(1960)  early  suggested  that  this  mechanism  may  be  responsible  for  observed 
LQ  waves  from  underground  explosions.  Press  et  al.  (1963)  also  presented 
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evidence  for  mode  conversion  in  the  case  of  nuclear  detonations  in  the 
air.  LQ  waves  observed  in  the  collapse  of  the  MILROW  explosion  on  Amchitka 
Island  (von  Seggern,  1973)  also  indicated  some  degree  of  mode  conversion 
since  this  type  of  event  is  presumed  to  be  symmetrical  (Smith,  1963).  For 
the  Nevada  Test  Site,  however,  Brune  and  Pomeroy  (1963)  noted  minimal  LQ 
waves  from  several  collapses  associated  with  underground  detonations  where 
LQ  was  seen.  In  all  cases  the  collapse  LQ/LR  ratios  were  much  less  than  those 
for  explosions,  indicating  that  mode  conversion,  if  operative,  accounts  for 
only  a  small  portion  of  the  SH  motion  recorded  for  explosions. 

Perhaps  the  most  convincing  argument  against  mode  conversion  as  a  cause 
of  explosion  SH  and  LQ  waves  is  the  success  in  inferring  double-couple  mechan¬ 
isms  for  strain  release  or  earthquake  triggering  by  underground  nuclear  deto- 

II 

nations  (e.g.,  Toksoz  and  Kehrer,  1972a;  Aki  and  Tsai  1972).  Random  genera- 
tion  of  SH  motion  by  mode  conversion,  well  away  from  the  source  along  the 
path  to  the  detecting  station,  could  not  be  expected  to  result  in  the  type  of 
patterned  data  that  has  been  reliably  inverted  in  numerous  cases  to  get  these 
mechanisms.  Other  evidence  in  the  form  of  dilatational  first  motions  and 
effects  on  Rayleigh-wave  amplitudes  (to  be  discussed  later)  integrate  well 
with  observed  SH  and  LQ  motion  to  prove  that  it  is  generated  primarily  at  the 
source  and  not  by  mode  conversion. 

Commenting  on  effects  within  the  source  itself,  Wright  and  Carpenter 
(1962)  offered  an  explanation  for  the  generation  of  SH  by  showing  the  asymme¬ 
tric  plaster  casts  of  cavities  formed  by  small  detonations  in  soft  clay.  They 
tentatively  explained  that  the  asymmetry  was  the  result  of  instabilities  that 
arise  when  a  light  material  (explosion  gases)  accelerated  a  heavy  one  (surround- 
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ing  solid)  as  described  by  Taylor  (1950).  The  theoretical  work  of  Enescu  et 
al.  (1973)  attempted  to  model  this  type  of  phenomenon  by  adding  three  unequal, 
mutually  perpendicular  forces  to  the  three  mutually  perpendicular  and  equal 
dipoles  ordinarily  used  for  an  explosion.  Lacking,  however,  was  evidence  for 
the  magnitude  of  this  effect  in  large-scale  explosions,  and  it  is  most  likely 
small  in  comparison  to  other  remaining  explanations  for  SH  generation. 

One  remaining  explanation  of  SH  generation  not  requiring  non-uniform  pre¬ 
stress  around  the  explosion  would  be  the  effects  of  radial  cracking  that 
Kisslinger  et  al.  (1961)  suggested  in  reporting  on  a  series  of  explosions  in 
soil.  Grover  (1973),  on  the  basis  of  model  experiments,  supported  radial 
cracking  as  the  cause  of  significant  SH  motion  from  explosions.  Again,  evi¬ 
dence  of  such  a  mechanism  is  missing  for  large-scale  nuclear  detonations  where 
observable  effects  in  support  of  some  mechanism  that  requires  shear  or  non- 
uniform  prestress  are  dominant.  Indeed,  explanations  for  SH  which  invoke  no 
prestress  cannot  account  for  a  level  of  SH  generation  that  will  hinder  dis¬ 
crimination.  Only  mechanisms  associated  with  a  release  of  non-uniform  or 
shear  prestress  can  generate  the  high  SH  and  LQ  amplitudes,  equivalent  to  or 
greater  than  P  and  LR,  which  are  observed  on  explosion  recordings  such  as 
for  PILE  DRIVER,  HARDHAT,  BILBY,  GREELEY,  or  BENHAM  at  the  Nevada  Test  Site. 

Press  and  Archambeau  (1962)  offered  the  first  of  these  prestress  explana¬ 
tions.  It  involved  relaxation  of  ambient  shear  stress  within  a  volume  sur¬ 
rounding  the  explosion  cavity  where  rigidity  vanishes  ("cavity"  includes  all 
of  the  inelastic  zone).  However,  actual  calculations  showed  that  only  liberal 
values  of  the  elastic  radius  (300  ra)  and  of  the  ambient  stress  level  (^10  bars) 
could  account  for  the  observed  SH  motion  from  RAINIER.  Later,  Archambeau 
(1972)  calculated  that  observed  SH  motion  could  be  predicted  by  increasing 
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the  stress  level  by  roughly  an  order  of  magnitude  and  including  stress  relax¬ 
ation  within  a  volume  having  a  radius  roughly  four  times  the  elastic  radius 
of  several  explosions  at  the  Nevada  Test  Site.  Toksoz  et  al.  (1965)  found 
that  HAYMAKER  and  SHOAL  LQ  could  be  accounted  for  by  stress  relaxation. 
Archambeau  and  Sammis  (1970)  calculated  that  the  BILBY  anomalous  radiation 
field  could  be  satisfied  with  the  stress  relaxation  theory,  if  420  m  is 
assumed  for  the  elastic  radius  and  75  bars  is  assumed  for  the  ambient  stress 
level.  All  these  calculations  assumed  complete  stress  relaxation. 

Archambeau  did  not  specify  the  actual  mechanism  of  stress  relaxation; 
it  could  be  either  isotropic  or  along  specific  cracks  as  Toksoz  et  al.  (1971) 
have  suggested.  A  further  confinement  of  the  stress  relaxation  about  a  single 
plane  is  a  possibility.  Separation  of  this  plane  from  the  cavity  brings  us  to 
the  "triggered  earthquake"  model  in  which  the  shear  stresses  of  an  explosion 
trigger  a  distinct  tectonic  event.  This  situation  appears  to  be  the  case  for 
HARDHAT  (Toksoz  et  al.,  1965)  and  for  BENHAM  (Aki  et  al.,  1969).  In  both  cases 
the  anomalous  radiation  field  is  so  large  that  it  dominates  the  explosion  field. 
Such  explosions  are  characterized  by  high  "F"  factors  (ratio  of  tectonic-to- 

If 

explosion  source-function  amplitudes  — see  Toksoz  and  Kehrer,  1972a).  Archam¬ 
beau  (1972)  demonstrated  the  equivalence  of  the  stress  relaxation  theory  to  a 
shear  dislocation  commonly  used  in  models  of  earthquake  source  mechanisms. 

Thus  this  theory,  corrected  by  Snoke  (1976),  is  probably  compatible  with  most 
observations  of  SH  waves  from  explosions,  regardless  of  the  actual  release 
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mechanism.  Consequently,  the  usual  double-couple  model  of  earthquake  source 
mechanisms  can  be  applied  to  SH  and  LQ  observations  from  explosions. 

A  considerable  number  of  studies  support  the  accuracy  of  the  tectonic 
strain  release  model,  and  among  them  are  numerous  successful  interpretations 
of  SH  and  Love  waves  from  explosions  (Toksoz  et  al.,  1965;  Brune  and  Pomeroy, 
1963;  Nuttli,  1969;  Lambert  et  al.,  1972;  Hirasawa,  1971;  Toksoz  et  al.,  1971) 
Much  evidence  also  exists  to  support  the  generation  of  patterns  like  those 
observed  from  underground  nuclear  explosions  from  model  experiments  with  explo 
sions  in  prestressed  sheets  of  material  (Kim  and  Kisslinger,  1967;  Grover, 

If 

1973;  and  Toksoz  et  al.,  1971).  Observations  of  fault  displacements  (McKeown 
and  Dickey,  1969)  near  the  explosion  site  further  support  tectonic  strain 
release  as  the  cause  of  observed  SH  and  LQ  waves  from  explosions.  In  addition 
the  periodic  and  natural  release  of  accumulated  strain  in  the  crust  at  the 
Nevada  Test  Site  by  shallow  earthquakes  (Molnar  et  al.,  1969;  Slemmons  et  al., 
1965)  also  supports  tectonic  strain  release  by  explosions  in  this  area. 

However,  exceptions  to  the  theory  exist,  notably  the  Amchitka  Island 
test  site  which  is  aseismic  (considering  only  shallow  crustal  seismicity  at 
Amchitka).  Toksoz  and  Kehrer  (1972b)  calculated  F  factors  for  the  MILROW  and 
CANNIKIN  explosions  equal  to  ones  for  typical  large  explosions  at  the  Nevada 
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Test  Site.  The  cause  of  the  large  LQ  waves  from  Amchitka  shots  is  not  clear. 
Toksoz  and  Kehrer  (1972a)  also  inferred  strain  release  at  the  Kazakh  and 
Novaya  Zemlya  test  sites.  Although  the  F  factors  for  these  regions  are  in 
the  same  range  as  NTS  shots,  strain  release  is  unexpected;  and  other  mechanisms 
which  do  not  require  shear  prestress  are  probably  responsible  for  SH  generation. 

One  other  model  to  be  considered  is  that  of  Andrews  (1973),  who  studied 
stable  sliding  on  pre-existing  faults.  This  sliding  resulted  from  reduced 
normal  stresses  on  faults  caused  by  passage  of  compressional  waves  or  the 
interaction  of  the  surface  reflection  with  the  primary  P  wave,  like  for  spal¬ 
ling  occurrences  (Eisler  et  al.,  1966).  This  model  is  also  a  trigger  model, 
but  differs  from  Toksoz  and  Kehrer's  because  only  compressional  waves  are 
required  and  because  tension  forces,  not  shear  forces,  are  considered  as  critical 
on  the  fault.  Andrews  shows  the  feasibility  of  his  model  with  numerical  finite- 
difference  calculations. 

Despite  the  less  than  ideal  behavior  of  explosions  in  terms  of  S  waves, 
discrimination  work  has  often  successfully  focused  on  these  waves.  Willis  et 
al.  (1963)  reported  that  S/P  amplitude  ratios  for  explosions  at  the  NTS  were 
roughly  one-third  of  that  for  earthquakes  in  the  surrounding  region.  Press  et  al. 
(1963)  and  Booker  and  Mitronovas  (1964)  also  showed  that  SH  and  LQ  are  lower  for 
explosions  than  earthquakes  in  this  region.  Evernden  (1969),  in  examining 
recordings  from  all  known  nuclear  explosions  through  1966,  found  only  two  long- 
period  S  waves  and  concluded  that  the  difference  in  generation  of  long-period 
S  for  equivalent  m^  earthquakes  and  explosions  must  be  at  least  an  order  of 
magnitude.  Nuttli  (1969),  in  a  comprehensive  analysis  of  LRSM  and  WWSSN  record¬ 
ings  of  HALFBEAK  and  GREELEY,  reported  131  long-period  S  waves,  which  gave  mag¬ 
nitudes  roughly  one  unit  less  than  the  short-period  P-wave  magnitudes  for  these 
two  events.  Blandford  and  Clark  (1974)  reported  40  long-period  S  measurements 
from  several  large  shots  at  NTS  and  Amchitka  Island  and  suggested  that  long- 
period  S-wave  magnitude  versus  Mg  is  a  good  discriminant.  However,  because  m^ 
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versus  Mg  is  an  excellent  discriminant,  m^  versus  long-period  S-wave  magnitude 
should  be  as  good  or  even  better.  By  reviewing  Seismic  Data  Laboratory  reports 
on  NTS  shots  and  LRSM  and  VELA  Observatory  earthquake  bulletins,  histograms  of 
S-wave  amplitude  over  P-wave  amplitude  for  earthquakes  and  explosions  have  been 
constructed  as  shown  in  Figure  1  (details  of  the  data  base  can  be  found  in  von 
Seggern,  1972).  Although  the  explosion  sample  is  very  small,  because  of  the 
limitation  of  background  noise,  the  sample  is  not  biased.  In  fact,  the  true 
density  may  lie.  toward  lower  ratios.  Distribution  for  a  subset  of  large  earth¬ 
quakes  is  shown  separately  to  prove  that  the  earthquake  histograms  are  unbiased 
by  threshold  problems.  Roughly  an  order  of  magnitude  separation  appears  for 
both  cases  between  explosions  and  earthquakes.  Although  noise  obviously  hampers 
application  of  this  discriminant,  recent  advances  in  the  theory  of  negative 
discriminants  (identification  by  absence  of  an  observation)  should  extend  its 
usefulness  to  lower  magnitudes. 

Because  LQ  waves  are  representative  of  the  long-period  SH  spectrum, 
they  will  also  serve  as  useful  discriminants.  Pertinent  to  this  supposition, 
another  histogram  from  von  Seggern  (1972)  is  shown  in  Figure  2.  The  explo¬ 
sions  are  again  from  NTS  where  known  tectonic  strain  release  repeatedly  gener¬ 
ates  significant  LQ  waves,  and  evidently  little  separation  exists  between  obser¬ 
ved  LQ/LR  ratios  for  earthquakes  and  explosions  from  this  area.  However,  since 
LR,  in  terms  of  Mg,  is  an  excellent  discriminant,  LQ  is  inferred  to  be  as  good 
or  better,  especially  for  shots  where  tectonic  strain  release  is  not  as  great 
as  at  NTS.  Savino  et  al.  (1971)  demonstrated  this  theory  by  estimating  Mg  for 
many  events  from  both  the  LQ  and  the  LR  waves.  Further  discriminating  infor¬ 
mation  exists  in  the  spectra  of  the  LQ  waves  because  apparently  there  is 
relatively  more  long-period  energy  in  the  LQ  waves  than  the  LR  waves  from  NTS 
events  (Lambert  et  al.,  1972;  Toksoz  et  al.,  1971),  indicating  a  tectonic 
origin  for  the  LQ  waves  as  opposed  to  the  dominantly  explosive  origin  for  the 
LR  waves.  As  for  S  waves,  noise  hampers  application  of  LQ  discriminants  and 
suggests  that  negative  discrimination  should  be  employed. 

von  Seggern,  D.  H.  (1972).  Seismic  shear  waves  as  a  discriminant  between  earth¬ 
quakes  and  underground  nuclear  explosions.  Seismic  Data  Laboratory 

Report  No.  295,  Teledyne  Geotech,  Alexandria,  Virginia. 


Figure  2  LQ/LR  amplitude  ratios  for  NTS  explosions  and  global  earthquakes. 
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This  discriminant  not  only  gives  excellent  results  in  separating  earth¬ 
quakes  and  explosions,  with  a  low  probability  of  misidentif ication  (Ericsson, 
1970;  Weichert  and  Basham,  1973),  but  it  is  also  easy  to  estimate  routinely. 
Normal  application  involves  calculating  m^  from  the  short-period  P-wave  ampli¬ 
tude  near  the  first  arrival  and  Mg  from  the  long-period  LR  amplitude  near  a 
period  of  20  seconds.  A  number  of  variations  have  been  investigated,  includ¬ 
ing:  1)  AR —  a  measure  of  the  envelope  area  of  the  LR  wave  (Brune  and  Pomeroy, 

1963) ;  2)  ERZ —  an  integral  of  the  energy  spectrum  of  the  LR  wave  (Lambert  et 
al.,  1969);  3)  Mq —  a  measure  of  the  source  moment  from  long-period  signals 

(Hanks  and  Thatcher,  1972);  4)  M  (40) —  an  estimate  of  M  at  a  period  near 

s  s 

40  seconds  rather  than  20  seconds  (Savino  et  al.,  1971).  Since  these  varia¬ 
tions  are  nearly  equivalent  to  the  normal  M^-versus-m^  discriminant  in  that 
they  all  represent  the  generation  of  long-period  LR  waves  relative  to  short- 
period  P  waves,  they  will  not  be  discussed  separately. 

Reasons  for  the  success  of  Mg  versus  m^  can  be  organized  into  four  cate¬ 
gories:  1)  source  geometry;  2)  source  location;  3)  time  dimensions  of  the 
source;  and  4)  spatial  dimensions  of  the  source.  Table  I  lists  these  cate¬ 
gories  along  with  some  of  the  more  important  references  contributing  to  the 
understanding  of  each. 

Beginning  with  elementary  concepts,  first  considered  are  the  effects  on 
Hs-versus-m^  of  the  geometry  of  simple  force  systems  representative  of  explo¬ 
sions  and  earthquakes.  The  uniform  compression  of  an  explosion  in  a  cavity 
can  be  modeled  with  three  mutually  perpendicular  and  equal  couples  without 
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TABLE  I 


Cause 

Source 

Equivalent 

Forces 


Source 

Time 

Function 


Source 

Depth 


Source 

Dimensions 


Possible  Causes  of  the  Mg-mb  Separation 

Additional 

Proponents  Discussions 


Gilbert  (1973) 

Douglas  et  al.  (1971) 
Leet  (1962) 


Davies  and  Smith  (1968) 
Marshall  (1970) 

Molnar  (1971) 

Peppin  and  McEvilly  (1974) 


Hudson  and  Douglas  (1975) 


Press  et  al.  (1963) 
Douglas  et  al.  (1972b) 
Tsai  and  Aki  (1971) 
Rodean  (1971) 


Peppin  (1976) 

McEvilly  and  Peppin  (1972) 
Liebermann  and  Pomeroy  (1969) 
Savino  et  al.  (1971) 

Tsai  and  Aki  (1971) 

Aki  (1972) 

Muller  (1973) 

Aki  et  al.  (1974) 

Kogeus  (1968) 

Toksbz  et  al.  (1964) 


Douglas  et  al.  (1971) 
Douglas  et  al.  (1972b) 


Liebermann  and  Pomeroy  (1969)  Peppin  (1976) 
Liebermann  and  Pomeroy  (1970)  Aki  (1972) 
Wyss  et  al.  (1971) 

Evernden  (1975) 


moment  (Love,  1944;  Maruyama,  1963).  An  earthquake  can  be  modeled  with  two 
such  couples  with  equal,  but  opposite,  moment  (Maruyama,  1963;  Burridge  and 
Knopoff,  1964).  The  far-field  amplitudes  for  these  sources  will  be  expressed 
in  terms  of  the  quantity  moment. 


The  relative  excitation  of  P  waves  for  the  two  types  of  sources  will  be 
computed  on  the  basis  of  the  long-period  portion  of  the  spectra.  Thus,  the 
result  will  have  validity  only  for  m,  measured  at  frequencies  lower  than  the 

D 

corner  frequencies  that  characterize  each  spectrum.  For  an  earthquake,  the 
source  modeled  (as  described  above)  gives  the  far-field  P-wave  spectral  ampli¬ 
tude  as  (Randall,  1973;  Dahlen,  1974) 
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in  a  homogeneous  infinite  medium,  where  is  the  moment  of  one  of  the  two 
equivalent,  orthogonal  couples,  p  is  density,  a  is  compressional  velocity,  r 
is  distance,  and  =  sin  20  cossj>  is  the  radiation  pattern  with  0  and  i|>  the 
polar  and  meridian  angles  of  a  spherical  coordinate  system  with  cf>  in  the 
plane  of  the  fault.  For  an  explosion  model,  the  far-field  spectral  amplitude 
is  (Muller,  1973) 
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where  is  the  strength  of  one  of  three  mutually  perpendicular  dipoles.  If 
the  moments  of  the  two  sources  are  equal  a  simple  ratio  is  defined: 


GQ  (w*0) 

_£ _ 

(ur*0) 
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(pot3)Q 


Of  special  interest  is  the  rms  value  of  R  ,  that  is,  a  value  representing  a 

fairly  uniform  sampling  of  the  focal  sphere  by  a  network  of  stations  estimat- 

2 

ing  m^  for  the  earthquake.  Thus,  R  is  integrated  over  the  surface  of  a 
sphere  (Appendix  I)  to  get  R2  =  4/15  for  P  waves 

edly  overestimated  because  stations  near  nodes  do  not  detect  the  P  wave 
the  result  for  events  of  equal  moment  is 


In  practice  R  is  undoubt- 

0<P 

Thus, 


uQ  (w*0)  =  • 52  ((tf+O)  -,P,0t  /~Q 

P  p  (pa3)^ 

The  difference  in  Rayleigh-wave  excitation  for  the  two  sources,  due  to 
source  type  alone,  can  be  calculated  by  taking  the  theoretical  excitation 
expression  for  a  double-couple  from  Ben-Menahem  and  Harkrider  (1964) 
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where  W  are  vertical  displacements  on  a  plane-layered  earth  model,  k  is  wave- 
number,  A^  is  the  spectral  amplitude  response  of  a  layered  medium,  and 
|x(9,A,6,h)|  is  the  radiation-pattern  factor  at  a  particular  azimuth  0  for  a 
fault  with  slip  angle  A  and  dip  angle  6  at  a  depth  h.  The  corresponding 
explosion  expression,  given  by  Harkrider  (1964) ,  can  be  approximated  for  a 
source  near  the  surface  by 
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where  P  is  the  residual  pressure,  a  is  the  radius  of  the  inelastic-elastic 
o 

boundary  near  the  source,  and  is  Rayleigh-wave  ellipticity  at  the  shallow 
source  depth  h.  The  elementary  relation  between  pressure  Pq  in  the  cavity  of 
radius  a  and  the  moment  of  the  explosion  is  given  by  Love  (1944,  p.  187) 
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Thus,  we  have  for  the  explosion 
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Again,  assuming  equal  for  the  two  types  of  sources,  the  ratio  of  Rayleigh- 
wave  amplitudes  is 
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This  relation  is  given  in  equivalent  form  for  a  homogeneous  half-space  by 
Douglas  et  al.  (1974a).  Just  as  for  the  earthquake  P  wave,  the  rms  of  the 
radiation  pattern  factor  xR  is  adopted,  which  is  given  in  Ben-Menahem  and 
Harkrider  (1964).  To  get  rms  xD  one  integrates  the  expression  for  all  values 
of  slip  angle  A  and  dip  angle  6,  as  well  as  0;  from  Appendix  I  the  result  of 
this  integration  for  a  surface  source  at  a  period  of  20  seconds  where  is 
commonly  estimated  is 
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Substituting  this  in  equation  (2),  using  a  =  6.03  and  8  =  3.53  for  a  Gutenberg 
earth  model,  and  letting 
earth  model,  one  obtains 


earth  model,  and  letting  =  .67  for  20-second  Rayleigh  waves  in  a  Gutenberg 
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This  equation  is  an  approximation  for  near-surface  sources  modeled  by  simple 
equivalent  force  systems,  and  it  does  not  include  depth  effects  or  effects 
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due  to  extended  time  or  space  dimensions  for  the  source.  Combining  equations 
(1)  and  (3)  and  taking  logarithms,  the  theoretical  Mg  -  difference  is  pro¬ 
duced  for  events  having  arbitrary  moments,  not  necessarily  equal,  and  near 
the  surface  in  the  same  medium: 

(Mg  -  <VQ  =  (Ms  '  mb)X+  ,22 
or  for  events  of  equal  nr 


M  Q  •  M  X  +  .22 

s  s  (4) 

This  result  differs  from  Gilbert's  (1973),  who,  in  a  derivation  that  was  elegant 

but  still  dependent  upon  several  crucial  but  not  obviously  reasonable  approx- 

Q  X 

imations,  predicted  M  s  M  +  .52.  The  derivation  here  is  considered  more 

s  s 

appropriate  because  it  takes  into  account  the  exact  theoretical  excitation  for 
Rayleigh  waves  in  a  layered  earth  model  due  to  simple  point  sources.  Lxolosions 
generate  nearly  as  large  Rayleigh  waves  (factor  of  1.7  less)  as  do  typical  shal¬ 
low  earthquakes  because,  although  the  earthquake  mechanism  produces  additional 
S  waves  and,  therefore,  presumably  more  Rayleigh  waves,  much  of  this  additional 
S  is  instead  converted  to  Love  waves.  If  equation  (4)  is  approximately  true, 
then  a  need  exists  to  find  through  other  considerations  realistic  causes  of 
the  large  Mg  difference  which  is  usually  observed  to  be  roughly  one  magnitude 
unit  or  more  between  earthquakes  and  explosions  of  equivalent  m^.  Note  that 
an  alternative  is  to  consider  the  cause  of  the  m^  difference  at  equivalent  Mg. 
Although  the  procedure  followed  in  the  derivation  of  equation  (4)  was  to  take 
the  logarithm  of  the  rms  amplitude  whereas,  in  practice,  M&  and  m^  are  cal¬ 
culated  by  taking  the  arithmetic  mean  of  the  logarithms  of  the  amplitudes,  the 
result  should  not  be  affected  greatly  by  considering  instead  this  second  method 
of  calculation,  which  must  be  treated  numerically  rather  than  analytically. 

Effects  of  Source  Depth  on  Mg  versus  m^ 

To  evaluate  how  depth  of  source  affects  the  Mg-  m^  relation  through  the 
body-wave  and  Rayleigh-wave  amplitude,  equations  (1)  and  (2)  can  be  generalized 
to  any  depth  by  using  the  appropriate  depth-dependent  parameters.  The  explo¬ 
sion  source  will  be  kept  at  the  surface  while  evaluating  | x„(9,X,6,h) j ,  given 
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in  Appendix  I,  for  some  discrete  depths  in  a  Gutenberg  earth  model  used  by 
Ben-Menahem  and  Harkrider  (1964).  The  rms  radiation  will  be  used  again  for 
Rayleigh  u  *es  since  interest  lies  in  gross  dependence  on  depth.  Also  evalu¬ 
ated  is  the  rms  Love-wave  radiation  as  expressed  in  Appendix  I.  According  to 
Harkrider,  the  actual  excitation  ratio  of  Love  and  Rayleigh  waves  for  a  double¬ 
couple  point  source  is  dependent  on  the  factors  A^,  A^,  k^,  and  k^,  thus 
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The  results  of  these  computations  for  a  period  of  20  seconds  are  listed  in 

1/2  1/2 

Table  II.  The  factor  k  A. /k  A  is  .86  for  this  case.  Figures  3,  4,  and 

L  L  K  K 

5  illustrate  the  essential  Mg  -  m^  character  as  the  earthquake  becomes  more 
deeply  located.  Figures  3  and  4  show  that  the  near  equivalence  of  Mg  values 
for  a  surface  explosion  and  a  surface  earthquake  of  equal  moment  still  holds  at 
earthquake  depths  of  30  km  and  that  the  LQ/LR  ratio  for  earthquake  is 
on  the  order  of  one  near  the  surface  but  nearer  two  at  depth.  This  result 
provides  evidence  that  much  of  the  S-wave  excitation  of  earthquake  sources  is 
converted  into  Love  waves,  making  it  unavailable  to  increase  Rayleigh-wave 
amplitude  relative  to  that  of  an  explosion.  Figure  4  illustrates  the  strong 
influence  of  the  pa3  factor  on  P-wave  amplitude.  Bouchon  (1976)  and  Hudson  and 
Douglas  (1975)  have  emphasized  this  effect  in  regards  to  Mg  -  m^  discrimination 
and  yield  estimation.  Figure  5  illustrates  the  predicted  effect  of  depth  on 
the  Mgf -  m^  discriminant  in  a  realistic  Gutenberg  earth  model.  All  explosion 
calculations  here  were  made  for  zero  depth,  but  evaluation  for  this  earth 
model  of  the  full  explosion  factor  of  Harkrider  (1964)  to  depths  of  4  or  5 
kilometers  would  not  alter  this  figure's  predictions  by  more  than  20%.  Figure 
5  indicates  that  earthquake  depth  does  not  improve  the  Mg  -  m^  discriminant 
much  and  that,  in  fact,  for  a  certain  range  of  shallow  depths,  the  rms  Ray¬ 
leigh-wave  excitation  of  earthquakes  is  diminished.  Note  that  the  Gutenberg 
model  contains  no  low-velocity  sediment  layers  near  the  surface.  (Later  in 
this  report  a  model  appropriate  to  the  Southwestern  United  States,  that  has 
such  layers,  will  be  examined.) 


Bouchon,  M.  (1976).  Teleseismic  body  wave  radiation  from  a  seismic  source  in 
a  layered  medium.  Geophys.  J.,  47,  515-530. 

Hudson,  J.  A.,  and  A.  Douglas  (1975).  On  the  amplitudes  of  seismic  waves. 
Geophys.  J.,  42 ,  1039-1044. 
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TABLE  II 


Theoretical  Calculations  of  Surface-Wave  and 
Body-Wave  Excitation  in  a  Gutenberg  Earth  Model 


QUAKE 

DEPTH 

1 

2 

3 

4 

5 

6 

(km) 

0. 

0.516 

0.354 

0.395 

0.894 

0.858 

0.221 

1. 

0.516 

0.354 

0.330 

0.072 

0.727 

0.149 

2. 

0.516 

0.354 

0.283 

1.250 

0.629 

0.086 

3. 

0.516 

0.354 

0.246 

1.440 

0.550 

0.028 

4. 

0.516 

0.355 

0.223 

1.595 

0.501 

-0.013 

5. 

0.516 

0.356 

0.216 

1.647 

0.484 

-0.028 

6. 

0.516 

0.357 

0.220 

1.627 

0.489 

-0.024 

7. 

0.516 

0.359 

0.233 

1.542 

0.515 

-0.001 

8. 

0.516 

0.360 

0.254 

1.418 

0.559 

0.034 

9. 

0.516 

0.362 

0.274 

1.323 

0.593 

0.060 

10. 

0.516 

0.365 

0.297 

1.227 

0.632 

0.088 

11. 

0.516 

0.367 

0.323 

1.136 

0.676 

0.117 

12. 

0.516 

0.370 

0.352 

1.051 

0.724 

0.147 

13. 

0.516 

0.373 

0.375 

0.994 

0.752 

0.163 

14. 

0.516 

0.377 

0.401 

0.941 

0.782 

0.180 

15. 

0.516 

0.381 

0.428 

0.891 

0.812 

0.197 

16. 

0.516 

0.386 

0.457 

0.844 

0.844 

0.213 

17. 

0.516 

0.391 

0.481 

0.813 

0.859 

0.221 

18. 

0.516 

0.397 

0.506 

0.784 

0.874 

0.229 

19. 

0.516 

0.403 

0.534 

0.756 

0.889 

0.236 

20. 

0.516 

0.411 

0.564 

0.729 

0.905 

0.244 

21. 

0.479 

0.412 

0.564 

0.731 

0.873 

0.260 

22. 

0.446 

0.413 

0.565 

0.731 

0.844 

0.277 

23. 

0.415 

0.415 

0.568 

0.730 

0.816 

0.294 

24. 

0.387 

0.416 

0.572 

0.727 

0.791 

0.311 

25. 

0.361 

0.418 

0.579 

0.723 

0.767 

0.328 

26. 

0.349 

0.422 

0.591 

0.713 

0.753 

0.334 

27. 

0.349 

0.428 

0.610 

0.701 

0.745 

0.330 

28. 

0.349 

0.434 

0.630 

0.690 

0.738 

0.325 

29. 

0.349 

0.442 

0.652 

0.678 

0.731 

0.321 

30. 

0.349 

0.450 

0.676 

0. 666 

0.723 

0.317 

31. 

0.349 

0.460 

0.703 

0.654 

0.716 

0.312 

32. 

0.349 

0.470 

0.732 

0.642 

0.709 

0.308 

33. 

0.349 

0.483 

0.765 

0.631 

0.702 

0.303 

34. 

0.317 

0.469 

0.733 

0.640 

0.646 

0.309 

35. 

0.288 

0.457 

0.704 

0.650 

0.595 

0.315 

36. 

0.263 

0.447 

0.678 

0.659 

0.549 

0.320 

37. 

0.240 

0.439 

0.656 

0.669 

0.506 

0.324 

38. 

0.220 

0.431 

0.636 

0.678 

0.468 

0.328 

39. 

0.202 

0.425 

0.618 

0.687 

0.433 

0.332 

40. 

0.185 

0.419 

0.603 

0.696 

0.400 

0.334 

41. 

0.185 

0.418 

0.609 

0.687 

0.390 

0.324 

39- 


TABLE  II  (Continued) 


Theoretical  Calculations  of  Surface-Wave  and 
Body-Wave  Excitation  in  a  Gutenberg  Earth  Model 


42. 

0.185 

0.417 

0.616 

0.677 

0.381 

0.313 

43. 

0.185 

0.416 

0.624 

0.667 

0.371 

0.302 

44. 

0.185 

0.415 

0.633 

0.657 

0.361 

0.290 

45. 

0.185 

0.414 

0.642 

0.645 

0.352 

0.279 

46. 

0.185 

0.413 

0.652 

0.634 

0.342 

0.267 

47. 

0.186 

0.412 

0.660 

0.624 

0.333 

0.254 

48. 

0.187 

0.413 

0.668 

0.618 

0.324 

0.240 

49. 

0.188 

0.413 

0.676 

0.612 

0.315 

0.226 

50. 

0.189 

0.414 

0.684 

0.605 

0.306 

0.211 

1 

P-WAVE  SPECTRAL 

AMPLITUDE 

RATIO  uQ/uX 

P  P 

2  ROOT  MEAN  SQUARE  OF  20-SECOND  LOVE  WAVE  RADIATION  PATTERN  xL 

3  ROOT  MEAN  SQUARE  OF  20-SECOND  RAYLEIGH  WAVE  RADIATION  PATTERN  X 

R 

4  RATIO  (rms  xL/rn>s  Xp) 

5  RAYLEIGH  WAVE  SPECTRAL  AMPLITUDE  RATIO  (rms  w^/rms  wX) 

R  R 

6  M  -m  DIFFERENCES:  (M  -m.  )Q  -  (M  -mv)X 

sb  sb  sb 

Note:  Explosion  is  near  surface. 


-40- 


RMS  EXCITATION  RATIO 


10  20  30  40  50  *0 

EARTHQUAKE  DEPTH  (km) 


Figure  4  Ratio  of  earthquake  surface-wave  and  body-wave  rras  excitation  to 
that  of  an  equivalent-moment  explosion  near  the  surface  versus 
depth  in  the  Gutenberg  earth  model. 
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Figure  5  Predicted  Ms-mt,  difference  versus  depth  between  earthquakes  and 
explosions  of  equal  moment  in  a  Gutenberg  earth  model. 


-43- 


The  effect  of  the  pP  interaction  on  m^  must  also  be  considered,  but  only 
for  very  shallow  events.  This  is  true  because  is  typically  measured  with¬ 
in  the  first  five  seconds  of  a  P-wave  signal,  and  so  the  free-surface  effect 
on  m^  measured  from  short-period  recordings  can  be  neglected  for  events  deeper 
than  roughly  ten  kilometers.  Assuming  zero-degree  incidence  for  explosions, 
the  free-surface  reflection  is  simply  the  P  signal  reversed  in  polarity  and 
equal  in  amplitude.  Clearly,  since  the  P  radiation  pattern  is  uniform  from 
an  explosion  and  the  basic  signal  pulse  is  only  one  cycle,  the  pP  can  raise 
or  lower  only  if  it  arrives  within  this  first  cycle.  Douglas  et  al.  (1971) 
have  adequately  discussed  the  effect:  for  zero  time  delay  of  the  pP  (explo¬ 
sion  right  at  the  surface),  the  P  signal  will  be  cancelled.  For  a  time  delay 
equal  to  roughly  one-half  cycle,  the  P  signal  will  be  reinforced  and  nearly 
doubled  in  amplitude.  Figure  6  shows  a  progression  of  cases  using  a  signal 
synthesized  for  a  10  kt  shot  as  recorded  by  a  short-period  Benioff  seismometer 
over  a  path  with  t*  =  .74.  Because  of  inelastic  effects  near  the  detonation, 
the  cancelling  of  P  motion  for  zero-depth  explosions  is  unrealistic;  and,  in 
fact,  the  equivalent  source  approaches  a  downward  point  force  applied  at  the 
surface,  rather  than  a  uniform  compress ional  force  (Kisslinger  et  al.,  1961). 
Douglas  et  al.  (1971)  also  indicated  the  effect  of  the  free  surface  on  earth¬ 
quake  m^.  Here  the  situation  is  more  complicated  because,  for  events  shallow 
enough  to  have  pP  recorded  within  the  first  five  seconds  of  the  direct  P,  the 
measurement  will  undoubtedly  reflect  the  larger  of  the  two  waves.  Thus, 
receivers  on  or  near  a  node  of  the  P-wave  radiation  pattern  may  be  situated  so 
that  the  pP  is  large  and  thus  still  report  a  normal  or  high  m^  for  the  earth¬ 
quake.  Therefore,  free-surface  reflection  will  produce  some  anomalies  in  m^ 
for  shallow  earthquakes.  Averaged  over  a  well-distributed  network,  the  sur¬ 
face  reflection  effect  for  earthquakes  should  be  to  increase  m^  at  most  one 
or  two  tenths,  and  at  worst  it  will  do  no  more  than  increase  the  explosion 
by  three  tenths.  In  summary,  depth  effects  on  P-wave  shape  cannot  signifi¬ 
cantly  contribute  to  separation  of  earthquakes  and  explosions  on  Mg  -  plots. 


Douglas,  A.,  J.  A.  Hudson,  and  V.  K.  Kembhavi  (1971).  The  relative  excita¬ 
tion  of  seismic  surface  and  body  waves  by  point  sources.  J,  Geophys. 

Res.,  23,  451-460. 

Kisslinger,  C. ,  E.  J.  Mateker,  and  T.  V.  McEvilly  (1961).  SH  waves  from  explo¬ 
sions.  Nature,  253,  242-245. 
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Effect  of  detonation  depth  on  signal  amplitude — perfect  reflection 


Effect  of  Source  Dimension  on  H  -  m 
_ s  b 

This  section  focuses  on  the  effects  which  the  finite  size  of  the  source 
has  on  the  spectra  of  seismic  signals  from  explosions  and  earthquakes. 

Foremost  among  these  effects  is  the  corner  frequency  in  the  spectrum  indicative 
of  the  spatial  extent  of  the  source.  Appropriate  expressions  derived  for  both 
cases  must  contain  source  dimension  so  that  their  effect  on  discrimination  is 
readily  perceived.  Thus,  the  explosion  and  earthquake  expressions  should  be 
obtained  in  as  nearly  an  equivalent  form  as  possible. 


For  the  explosion,  the  solution  for  far- field  displacement  is  easily 
attained  because  of  the  radial  symmetry  of  the  problem.  Jeffreys  (1931), 

Blake  (1952),  Sharpe  (1942),  and  Rodean  (1971),  among  others,  have  already 
performed  a  solution  of  the  wave  equation  subject  to  Cauchy  boundary  condi¬ 
tions  at  the  cavity  radius  (or  equivalent  elastic  radius).  The  common  solu¬ 
tion  is  in  terms  of  the  pressure  on  the  cavity  radii,  and  so  it  is  not  direct¬ 
ly  comparable  to  earthquake  models  which  usually  express  the  far-field  dis¬ 
placement  in  terms  of  displacement  at  the  source.  As  shown  in  Appendix  II, 
the  explosion  result  can  also  be  expressed  in  terms  of  a  step  source  dis¬ 
placement 


a  D 


.  (r  ,oj)  |  = 


U2 


„2 


(5) 


With  a  corner  frequency  at  w  -  a/a,  the  high-frequency  asymptote  is  w 
low  frequencies  the  value  of  the  flat  portion  of  the  spectrum  is 


!  u(r,0) 


azD 


ar 


At 


(6) 


Jeffreys,  H.  (1931).  On  the  cause  of  oscillatory  movements  in  seismograms. 
Monthly  Notices  of  the  Royal  Astronomical  Society,  8^,  408-416. 

Blake,  F.  C.  (1952).  Spherical  wave  propagation  in  solid  media.  J.  Acoust. 
Soc.  Amer . ,  24,  211-215. 

Sharpe,  J.  A.  (1942).  The  production  of  elastic  waves  by  explosion  pressures. 
Part  I.  Geophys ,  ,  _7,  144-154. 

Rodean,  H.  C.  (1971).  Nuclear-Explosion  Seismology.  USAEC  Division  of  Tech¬ 
nical  Information,  Oak  Ridge,  Tennessee. 
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Letting  u2  =  X  +  2p,  this  can  be  rearranged  to  give 

47rpa3rjC(r,0)|  =  M 

'  O 


where  Mq  =  4na2Do  pa2,  an  expression  agreeing  with  Aki  et  al.  (1974)  and 
Muller  (1973). 


For  an  earthquake,  a  finite  source  will  be  introduced  as  a  circular  plane 

where  the  slip  occurs  instantaneously  as  a  step  function  of  time.  The  formal 

expression  for  the  spectrum  of  the  far-field  P-wave  displacement  is  given  by 

Molnar  et  al.  (1973),  which,  for  future  convenience,  is  herein  divided  by  a 
o 

factor  of  sin  0: 


|Q(r,w)|  = 


2pa3r (wa*  sin0/a) 


2  Lwa*sin0 


sin 


■  wa*  sin0 


)- 


COS 


/(pa*  sin  Q 
'  a 


i] 


Where  the  polar  axis  for  0  is  perpendicular  to  the  circular  fault  plane  and 
R  ^  is  the  radiation  pattern  for  P  waves  from  a  double  couple  without  moment. 

O  9 

Here,  displacement  is  assumed  to  be  variable  over  the  fault  rather  than  con¬ 
stant  in  an  effort  to  more  closely  model  the  actual  mechanism  according  to 


(a2-a-*2)  2 


H(t) 


where  a'  (<a)  is  the  radius  at  which  D  „  (t)  is  calculated.  The  long-period 

3. 

asymptote  of  the  earthquake  spectrum  is 


,  URQ .  D  a2 

|Q(r,0)  j  =  — §4-  °  .  „ 

6pa3r 


Using  A=tra  for  the  fault  area,  this  can  be  rearranged  to  obtain 


4irpa3r  |  Q(r,0) 

Red> 


pAD 

o 


Aki,  K. ,  M.  Bouchon,  and  P.  Reasenberg  (1974).  Seismic  source  function  for 
an  underground  nuclear  explosion.  Bull.  Seism.  Soc.  Am.,  64,  131-148. 
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The  right-hand  side  is  exactly  the  moment  of  each  couple  of  the  equivalent 

2 

double-couple  point  source  if  we  replace  — Dq  by  Dq,  the  average  displacement 
on  the  fault,  which  can  be  obtained  by  integrating  equation  (7)  over  the  fault 
plane.  This  relation  between  the  far-field  long-period  portion  of  the  spec¬ 
trum  and  the  seismic  moment  holds  for  any  fault  shape  whose  area  is  A  and  any 
fault  displacement  whose  average  value  is  Dq  (Savage,  1972;  Randall,  1973; 
Dahlen,  1974).  The  corner  frequency  occurs  at  roughly  oja-sinO/ct  ~  1.  The 
minimum  corner  frequency,  when  sin0  =  1,  is  w  ~  ;  and  it  increases  as  sin  0-H) 

3 

until  the  spectrum  becomes  white,  corresponding  to  an  impulsive  P  wave  observed 
on  a  line  perpendicular  to  the  fault.  For  high  frequencies,  the  spectrum  is 
roughly  given  by: 

p 

I  f, ,  \  |  ,ua  .  ,ioa. 

|u(r,a>)  |  a  cos  (—)/(—) 

-2 

The  high-frequency  asymptote  of  the  spectrum  is  thus  proportional  to  w  .  For 

the  even  simpler  model  with  constant  displacement  D  over  the  entire  fault, 

-3/2  ° 

the  spectrum  can  be  shown  to  decay  as  w  at  high  frequencies.  Recall  that 
the  corner  frequency  of  the  explosion  model  with  step  displacement  is  u/a  and 
that  the  high-frequency  asymptote  is  u  \ 

The  effects  of  a  finite  source  on  explosion  and  earthquake  far-field 
P-wave  spectra  have  been  partially  determined.  The  finite  source  models 
already  examined  specified  a  step  function  for  displacement  over  the  source 
area,  normal  to  a  spherical  surface  for  the  explosion  and  parallel  to  a  cir¬ 
cular  surface  for  an  earthquake.  The  spectra  of  the  finite  sources  have  flat 

ot 

low-frequency  portions  and  corner  frequencies  at  roughly  w  =  — ,  the  main  dif- 

-2  -3/2  a-l 

ference  being  an  u  to  w  falloff  for  earthquakes  versus  u  for  explo¬ 
sions.  This  difference  can  be  attributed  to  the  geometrical  shape  of  the  two 
sources  that  govern  the  rate  of  growth  of  the  far-field  P  displacement  pulse. 
The  spherical  shape  of  the  explosion  model  leads  to  a  higher-order  discontinu¬ 
ity  in  the  far-field  pulse  than  the  circular  shape  of  the  earthquake  model. 


Savage,  J.  C.  (1972).  Relation  of  corner  frequency  to  fault  dimensions. 
J,  Geophys.  Res.,  77 ,  3788-3795. 
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Other  fault  shapes  have  been  considered  which  give  different  falloffs  (Has¬ 
kell,  1964;  Savage  1974;  Savage,  1966);  as  little  as  u  ^  falloff  can  be  pre¬ 
dicted  from  a  rectangular  fault,  when  the  observing  station  is  aligned  along 
a  parallel  to  one  edge.  Thus,  the  shape  of  the  source  is  important  in  pre¬ 
dicting  the  asymptotic  behavior  of  the  far-field  spectrum,  but  spectral  ampli¬ 
tude  difference  between  the  earthquake  model  with  instantaneous  rupture  and 
the  explosion  model  are  not  great  enough  to  enable  effective  discrimination 
to  be  based  upon  this  aspect  alone  because  attenuation  in  the  earth  has  a 
large  effect  on  shaping  the  higher-frequency  portion  of  body-wave  spectra. 
However,  such  attenuation  effects  could  be  removed  in  a  strictly  regionalized 
study,  where  both  types  of  events  were  located  within  a  homogeneous  province. 

Effect  of  Source  Time  Function  on  M  -  m, 

_ S  D 

Next  considered  are  source  time  functions  other  than  a  step.  A  ramp 
function  to  a  final  static  value  related  to  the  seismic  moment  is  probably 
a  satisfactory  first-order  approximation  for  earthquakes.  Brune  (1970),  Dah- 
len  (1974),  and  Haskell  (1964),  among  others,  have  all  considered  this  type  of 
function.  Its  effect  is  to  introduce  another  u  1  falloff  in  the  spectrum  of 
the  P-wave,  with  a  corner  frequency  equal  roughly  to  the  inverse  of  the  rise 
time  at  the  source.  That  this  rise  time  should  relate  to  the  fault  dimension 
divided  by  the  rupture  velocity  is  theoretically  appealing  because  the  motion 
of  the  fault  should  continue  until  the  rupture  reaches  the  edges  and  dies  out. 
Archuleta  and  Brune  (1975),  in  model  experiments,  demonstrated  that  such  a 


Haskell,  N.  A.  (1964).  Total  energy  and  energy  spectral  density  of  elastic 
wave  radiation  from  propagating  faults.  Bull,  Seism.  Soc.  Am.,  54, 
1811-1841. 

Savage,  J.  G.  (1974).  Relation  between  P-  and  S-wave  corner  frequencies  in 
the  seismic  spectrum.  Bull.  Seism.  Soc.  Am.,  64,  1621-1628. 

Savage,  J.  C.  (1966).  Radiation  from  a  realistic  model  of  faulting.  Bull. 
Seism.  Soc.  Am. ,  56 ,  577-592. 

Brune,  J.  N.  (1970).  Tectonic  stress  and  the  spectra  of  seismic  shear  waves 
from  earthquakes.  J,  Geophys,  Res.,  75,  4997-5010. 

Archuleta,  R.  J.,  and  J.  N.  Brune  (1975).  Surface  strong  motion  associated 
with  a  stick-slip  event  in  a  foam  rubber  model  of  earthquakes.  Bull. 
Seism.  Soc.  Am.,  65,  1059-1072. 
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time  function  occurs  in  a  simply  fractured  and  stressed  material.  Enough 
variations  can  arise  in  considering  different  rupture  velocities  and  fault 
shapes,  including  position  of  the  observer  relative  to  the  fault  plane,  that 
making  general  predictions  from  more  realistic  models  becomes  difficult. 
Dahlen's  (1974)  theoretical  model  is  especially  well-constructed  though;  and 
it  offers  some  simplicity  in  its  final  form,  predicting  that  the  corner  fre¬ 
quency  due  to  rise  time  and  the  corner  frequency  due  to  finiteness  of  the 
source  are  identical.  Taking  the  rupture  velocity  to  be  equal  to  V  everywhere 
and  assuming  a  Poisson  medium,  the  corner  frequency  for  Dahlen's  work  is 


1/3 


where  D  is  the  velocity  of  particle  movement  between  the  fault  planes  or  slope 

-3 

of  the  ramp  time  function.  The  high-frequency  asymptote  is  to  for  his  model. 

•  “  2  ”  2 

Using  D  =  .5VD  /a,  Mq  =  pA*D,  A  =  tt  a  ,  and  D  =  —D^,  one  can  reduce  this 
expression  to  roughly  1.4V/a  <  gj^  <  1.9V/a,  depending  on  0,  or  ioc  = 

For  the  explosion,  we  take  our  spectrum  given  by  equation  (5)  for  a  step 
function  at  the  source  and  generalize  to  any  time  function  by  replacing  Dq 
by  toDQ  |  D(to)  |  where  |D(to)|  is  the  amplitude  spectrum  of  an  arbitrary  source¬ 
time  function  D(t)  whose  final  value  is  scaled  to  unity. 

Now  specify  D(t)  to  be  a  ramp  function  similar  to  our  earthquake  model, 
with  rise  time  a/a.  Rodean  (1971,  p.  58)  has  shown  that  this  form  is  a  reason¬ 
able  first-order  approximation  to  the  source  function.  This  form  is  reasonable 
because  it  is  theoretically  appealing  to  have  each  point  on  the  cavity  wall 
(or  elastic  radius  limit)  move  outward  until  effects  from  all  other  points  of 
the  source  boundary,  but  which  are  not  occluded  by  the  cavity  itself,  pass  by 
after  traveling  a  distance  a  at  velocity  a.  With  this  assumption 

u  Dq  I  D(u)  |  ->  Dq  w  •>  0 

u  D  |d(<d)  |  ■*  D  /(au>/oO  w  -*•  00 
o  o 

Substituting  this  equation  into  equation  (5)  for  Dq  shows  that  the  low-fre¬ 
quency  asymptote  is  again  related  to  the  moment,  as  in  equation  (6),  and  that 
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Che  high-frequency  asymptote  is  proportional  to  to  and  that  the  corner  fre¬ 
quency  to  -  a/a  is  identical  to  that  of  the  simple  earthquake  model.  Although 

f  the  analytical  approach  was  different  in  each  case,  Mueller  and  Murphy  (1971) 

-2 

and  von  Seggern  and  Blandford  (1972)  have  previously  predicted  the  to  spec¬ 
trum  for  far-field  P  waves  from  explosions. 

Thus  far  rather  idealized  sources  have  been  employed,  at  least  in  the 

case  of  the  earthquake  model.  The  major  theoretical  findings  relevant  to  Mg 

-  discrimination  show  that,  given  similar  ramp  displacement  functions  at 

the  source,  the  earthquake  spectrum  will  fall  off  more  rapidly  beyond  the 

-3  -2 

*  corner  frequency  than  the  explosion  one,  to  versus  to  .  This  fact,  however, 

will  not  alter  the  Mg  -  predictions  if  the  corner  frequency  is  in  fact 

larger  than  roughly  1  cps,  which  is  the  typical  frequency  of  teleseismic  P 

waves  as  recorded  by  the  WWSSN  and  past  and  present  VELA  stations.  The 

at  which  the  corner  frequency  is  1  cps  can  be  found  by  using  the  following 

-1  1/3 

data:  Mueller  and  Murphy  (1971)  showed  empirical  evidence  that  a  2  10  Y 

Qt 

km  for  explosions  in  hard  rock  where  Y  is  in  kilotons.  Assuming  to  y  — , 
r  c  'v  a 

•  the  corner  frequency  can  be  expressed  as 

X  m  /vl/3 
to  *  10a/Y 
c 

Assume  a  2  3.5  km/sec  as  a  velocity  for  a  typical  hard  rock  detonation  medium; 

1/3  1/3 

this  gives  to  -  35/Y  or  f  ~  6/Y  .  Therefore,  the  corner  frequency  for 

c  t 

1  explosions  lies  above  1  cps  until  roughly  Y  =  200  kt  or  roughly  m^  =  6  as 

determined  empirically  (Evernden,  1970);  and  up  to  that  yield  m^  is  measured 
on  the  flat  portion  of  the  spectrum.  That  the  Mg  -  relation  is  nearly  linear 
with  slope  1.0  (Springer  and  Hannon,  1973)  from  small  yields  to  one  megaton 
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confirms  that  the  explosion  corner  frequency  is  beyond  1  cps  until  very  large 
yields.  Thus,  a  point-source  theory  for  explosions  is  a  good  approximation 
in  studying  the  Mg  -  difference  for  all  but  the  largest  detonations.  Sim- 
ilarily,  if  the  simple  coherent  fault  model  which  has  been  introduced  was 
applicable,  then  the  earthquake  corner  frequency  would  be  beyond  1  cps  until  a 
similarly  high  magnitude.  A  point-source  theory  would  also  apply  to  -  m^ 
determinations  for  earthquakes,  resulting  in  little  difference  between  them 
and  explosions.  Thus,  the  observed  large  Mg  -  difference  between  earthquakes 
and  explosions  at  m^  <  6  has  not  been  explained  by  requiring  finite  sources 
and  simple  dynamics. 

At  this  point,  two  possibilities  remain  for  enhancing  the  theoretical 
difference  between  earthquakes  and  explosions:  1)  in  contrast  to 
earthquakes,  the  source  time  function  for  explosions  may  have  a  very  short 
rise  time  and  may  have  a  significant  peak  before  it  reaches  its  static 
value  and  2)  the  source  dimensions  for  an  earthquake  may  be  significantly 
larger  than  those  of  an  explosion  of  the  same  mb>  The  first  explanation  is 
most  often  referred  to  as  differences  in  the  source  time  function,  and  the 
second  as  differences  in  source  dimensions. 

Several  researchers  have  utilized  the  source  time  function  explanation 
(see  Table  I).  However,  in  these  cases  separating  this  phenomenon  from  the 
effect  of  source  dimensions  that  actually  control  the  rise  time  of  the  source 
function  is  conceptually  difficult.  One  prediction  related  to  the  time 
function  states  that  if  the  rupture  propagation  on  a  fault  is  lowered  causing 
the  slope  of  the  ramp  time  function  to  be  lower  while  the  ramp  is  lengthened, 
this  subsonic  faulting  will  tend  to  lower  the  earthquake  corner  frequency  and 
enhance  Mg  -  m^  discrimination. 

Another  aspect  of  the  source  time  function  is  the  possibility  of  a 
significant  peak  in  the  explosion  spectrum.  This  possibility  is  suggested 
by  many  close-in  observations  of  displacements  around  an  underground  nuclear 
detonation  (Werth  et  al.,  1963; 


Werth,  G.  C. ,  R.  F.  Herbst,  and  D.  L.  Springer  (1962).  Amplitudes  of  seismic 
arrivals  from  the  M  discontinuity.  J.  Geophys.  Res.,  67,  1587-1610. 
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Perret,  1972;  Healy  et  al.,  1971).  These  close-in  observations  often  showed 
overshoots  in  the  displacement  that  are  several  times  the  apparent  static  or 
residual  value.  Since  cable  breaks  or  instrument  failure  typically  truncate 
the  recordings  within  two  seconds  of  detonation,  residual  displacement  cannot 
be  inferred  accurately.  This  is  unfortunate  because  shape  of  the  time  func¬ 
tion  out  to  tens  of  seconds  is,  of  course,  important  for  predicting  Rayleigh- 
wave  amplitudes.  In  cases  of  purely  elastic  response,  overshoots  of  such 
magnitude  are  not  predicted  by  theoretical  calculations  for  a  cavity  in  a 
whole  space,  but  are  apparently  allowable  for  cavities  near  the  surface  of 
a  half-space  (Holzer,  1965;  Aboudi,  1972).  Caution  must  be  taken  to  separate 
that  overshoot  part  of  the  source  function  which  is  due  to  the  surface  reflec¬ 
tion  effect  from  that  due  to  inelastic  behavior  since  the  first  should  affect 
both  shallow  explosions  and  earthquakes.  The  limit  of  inelastic  effects  in 
the  case  of  a  shallow  explosion  leads  to  cratering,  and  the  corresponding 
close-in  displacements  exhibit  large  overshoots.  However,  even  for  normal 
depth  of  burial  scaled  to  the  yield,  overshoot  by  a  factor  of  two  may  be 
realistic  in  certain  media.  For  example,  note  the  source  time  function  for 
the  NTS  MAST  explosion  calculated  with  a  finite-difference  stress  code  by 
Barker  et  al.  (1976).  It  is  suggested  here  that  observed  overshoots  in  close- 
in  data  are  due  largely  to  inelastic  behavior,  if  the  measurements  are 
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near  the  surface,  plus  an  elastic  surface-reflection  contribution  which  is  best 
described  by  numerical  computer  codes.  Thus,  the  effect  on  Mg  -  of  possible 
overshoots  in  the  source  time  function  for  explosions  is  most  likely  less  than 
0.3.  The  very  linear  Mg  -  relation  for  explosions  (Evernden  et  al.,  1971; 

Springer  and  Hannon,  1973;  Peppin  and  McEvilly,  1974)  seems  to  reinforce  this 
argument  because  if  a  large  overshoot  occurred,  it  would  result  in  a  large 
hump  in  the  far-field  spectrum  that  would  pass  through  the  narrow  frequency 
band  in  which  is  measured  as  the  yield  was  increased,  thus  leading  to  a 
noticably  non-linear  Mg  -  m^  relation. 

Note,  however,  that  considerable  empirical  evidence  exists,  mostly  from 

Rayleigh  waves  (Toksoz  et  al.,  1964;  Helmberger  and  Harkrider,  1972;  Aki  et 

al.,  1974),  for  the  presence  of  a  large  overshoot  in  the  source  function. 

Most  of  these  investigators  suggest  a  decaying  time  function  in  the  form 
b  — c  t 

kt  e  .  However,  the  decay  rate  that  they  required  cannot  be  accommodated 
by  theoretical  calculations  or  observations  for  the  pressure  inside  a  post¬ 
shot  cavity  (Olsen,  1967,  1970).  This  is  true  because  the  physical  properties 
do  not  allow  a  significant  pressure  decay,  unless  there  is  venting  or  a  cavity 
collapse.  Contrary  to  the  Rayleigh-wave  studies  just  mentioned. 
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Tsai  and  Aki  (1971)  have  concluded  from  Rayleigh  waves  that  the  source  time 
functions  for  explosions  and  earthquakes  are  nearly  identical.  Also,  Peppin's 
(1976)  close-in  spectra  for  many  NTS  explosions  showed  no  apparent  overshoot. 
Numerous  corrections  are  made  in  reducing  far-field  signal  spectra  back  to  the 
source  however,  and  all  investigators  must  admit  the  possibility  that  these 
corrections  are  not  precisely  known. 

Another  aspect  of  the  source  time  function  is  that  earthquake  sources 
are  coherent  only  for  periods  longer  than  their  fault  length  divided  by  the 
wave  velocity.  Spectral  amplitudes  for  periods  less  than  this  can  be  predict¬ 
ed  only  by  statistical  or  stochastic  methods,  and  they  are  necessarily  greater 
for  an  incoherent  model  of  faulting.  Haskell  (1964,  1966)  first  presented  this 
concept  in  an  analytical  fashion,  and  Aki  (1967)  worked  from  Haskell's  model 
to  show  that  such  incoherent  behavior  results  in  predicted  spectra  for  earth¬ 
quakes  that  agree  with  several  types  of  empirical  data  from  large  magnitude 
events.  The  phenomenon  of  incoherence  should  affect  Mg  -  m^  discrimination 
for  events  whose  source  dimensions  are  on  the  order  of  5  km  or  larger  since 
this  is  the  dimension  related  to  m^  measurement  at  1  cps.  Liebermann  and 
Pomeroy  (1970)  gave  varied  evidence  that  this  magnitude  lies  at  4.5  to  5.0. 
However,  on  the  basis  of  the  success  of  the  M  -  m,  discriminant  alone,  one 
should  be  reluctant  to  admit  any  significant  effect  for  incoherence  since 
its  effect  is  to  raise,  not  lower,  the  spectral  high-frequency  level. 

The  question  of  whether  earthquake  and  explosion  Mg  -  m^  lines  converge  at 

low  magnitudes  is  a  disputed  one  which  usually  rests  on  empirical  observations 

at  very  low  signal-to-noise  ratios.  It  has  been  shown  here  that  no  theoretical 

reason  exists  to  support  a  M  -  m  difference  between  explosions  and 

s  b 
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earthquakes  that  is  more  than  roughly  one-half  magnitude  unit  at  low  magni¬ 
tudes.  To  support  a  larger  difference,  one  must  hypothesize  a  combination  of 
subsonic  faulting  for  an  earthquake  and  perhaps  some  significant  overshoot  in 
the  explosion  source-time  function.  Thus,  Mg  -  m^  discrimination  is  a  tenuous 
proposition  for  small  events,  a  contention  that  Peppin's  (1976)  data  supports 
for  low-magnitude  explosions  and  earthquakes  at  the  NTS  since  it  shows  incom¬ 
plete  separation  of  source  type  by  Pg  vs.  LR  measurements. 

Surface-Wave  Spectra 

Derr  (1970)  and  von  Seggern  and  Lambert  (1970),  who  studied  a  large  num¬ 
ber  of  explosions  and  earthquakes,  have  shown  that  the  discrimination  capabil¬ 
ity  in  surface  waves  alone  is  minimal.  Even  assuming  significant  differences 
in  the  source  time  function  between  the  two  types  of  seismic  events,  the  pro¬ 
pagation  and  radiation  pattern  corrections  are  so  variable  that  differences  in 
observed  spectra  cannot  reliably  be  attributed  to  one  source  type  or  another. 
Source  depth  can,  however,  have  a  persistent  effect  on  observed  Rayleigh-wave 
spectra  if  the  earthquakes  are  all  deeper  than  explosions  in  the  same  area 
because  generation  of  longer  periods  is  enhanced  relative  to  shorter  periods 
as  a  source  is  moved  downwards.  Tsai  and  Aki  (1971)  and  Marshall  (1970)  dis¬ 
cussed  this  possible  means  of  discrimination.  The  Mg  difference  between  20- 
second  and  40-second  estimates  has  been  shown  by  Savino  et  al.  (1971)  to  be 
larger  for  earthquakes  than  for  explosions  in  general.  These  results  are 
most  likely  due  to  depth,  and  Marshall  and  Basham  (1972)  even  employ  the 
40-sec  versus  20-sec  Rayleigh-wave  ratio  to  determine  a  depth  correction  for 
Mg.  Tsai  and  Aki  (1971)  showed  that  deep  crustal  earthquakes  in  the  Southwest 
United  States  can  be  identified  as  such  on  the  basis  of  Rayleigh-wave  spectra. 
At  best,  Rayleigh-wave  spectra  can  only  be  described  as  an  aid  to  discrimina¬ 
tion  rather  than  as  a  reliable  discriminant  in  itself. 
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Although  discrimination  on  the  basis  of  short-period  P-wave  spectra  has 

been  demonstrated  by  several  studies  (Frantti,  1963;  Bakun  and  Johnson,  1970; 

Weichert,  1971;  Briscoe  and  Walsh,  1967;  Anglin,  1971;  Manchee,  1972;  Lacoss, 

1969;  Dahlman  et  al.,  1974),  physical  reasons  for  success  of  this  discriminant 

have  only  recently  become  apparent.  Because  pP  will  cancel  P  at  low  frequencies 

-2  -3 

and  because  the  spectrum  falls  off  as  to  for  explosions  versus  to  for  earth¬ 
quakes,  the  amplitude  at  high  frequencies  of  explosion  P  waves  is  enhanced  re¬ 
lative  to  that  of  earthquakes.  The  definition  of  spectral  ratio,  such  as  that 
of  Lacoss  (1969)  for  instance,  takes  advantage  of  this  inherent  difference. 
However,  this  discrimant  may  not  rest  on  a  firm  base  because  the  high-frequency 
portion  of  the  spectrum  (>  1  cps)  is  greatly  affected  by  various  propagation 
effects  which  may  dominate  source  characteristics.  Von  Seggern  and  Blandford 
(1976)  surveyed  spectral  ratios  from  globally-distributed  earthquakes.  Their 
results  are  reproduced  in  Figure  7,  which  shows  nearly  as  much  scatter  as 
typical  Mg  -  m^  plots  for  earthquakes.  The  few  presumed  explosions  are  not 
clearly  separated  in  this  work.  Therefore,  unless  events  are  in  proximity,  there 
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is  no  certainty  that  observed  differences  are  indicative  of  source  type. 

This  problem  relates  to  the  problem  of  large  scatter  in  m^,  which  can  be  attri¬ 
buted  either  to  Q  differences  (Ward  and  Toksoz,  1971;  Solomon,  1972a;  and 
Douglas  et  al. ,  1974a)  or  to  the  focusing  effects  of  lateral  inhomogeneities 
(Davies  and  Julian,  1972).  If  is  so  highly  dependent  on  path,  then  spectral 
ratios  also  may  be,  especially  if  higher  frequencies  are  employed. 

A  better  approach  with  short-period  spectral  ratios  is  to  emphasize  the 
frequencies  below  1  cps.  Here,  the  explosion  spectrum  should  decrease  rapidly 
as  a  result  of  the  effect  of  the  free  surface,  which  theoretically  applies  a 
factor  of  w  to  the  P-wave  spectrum  up  to  a  frequency  characteristic  of  the 
source  depth  (Douglas  et  al.,  1971).  The  notably  small  long-period  P  waves 
from  explosions  (von  Seggern,  1972;  Evernden,  1969)  are  evidence  of  this  effect. 
The  effect  on  the  earthquake  P  spectrum  is  more  complex  because  of  the  radia¬ 
tion  pattern.  However,  if  the  earthquakes  are  generally  deeper  than  the  explo¬ 
sions,  then  their  spectra  are  expected  to  remain  at  the  level  predicted  by  a 
whole-space  model  down  to  lower  frequencies  than  that  for  shallow  explosions. 
Molnar  (1971)  and  Wyss  et  al.  (1971)  observed  this  difference  in  explosion 
and  earthquake  spectra  and  suggested  it  as  a  discriminant.  Such  a  discriminant 
is  related  to  depth  determination  by  spectral  methods  (mentioned  above)  and 
will  be  subject  to  possible  pitfalls  associated  with  that  procedure.  In  gen¬ 
eral,  a  P-wave  spectral  ratio  employing  as  wide  a  frequency  band  as  possible 
should  be  a  good  discriminant,  but  less  reliable  than  Mg  -  m^» 
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Complexity  of  the  Recorded  Signal 

Early  comparisons  of  time  traces  for  explosions  and  earthquakes  revealed 
the  earthquake  signals  to  be  more  complex  in  appearance.  Carpenter  (1964) 
offered  an  explanation  in  terms  of  the  crustal  transfer  function  at  the  source 
which,  theoretically,  should  provide  more  numerous  and  stronger  secondary  arri¬ 
vals  within  thirty  seconds  after  the  initial  P  arrival  for  an  earthquake  at 
depth  than  for  an  explosion  near  the  surface.  In  opposition  to  this  idea, 
Douglas  et  al.  (1974a)  showed  that  a  thick  low-velocity  layer  in  the  crust 
is  required  to  generate  complex  earthquakes.  Douglas  et  al.  (1973),  in 
a  thorough  discussion  of  the  reasons  for  complexity,  concluded  that  Q  along 
the  path  may  be  the  major  influence  and  that  differences  in  complexity  between 
explosion  and  earthquake  signals  are  most  often  explained  by  their  emanating 
from  regions  of  greatly  different  Q,  that  is,  stable  versus  tectonic  provinces. 
Citing  LONG  SHOT  in  their  analysis,  Davies  and  Julian  (1972)  show  that  shadow¬ 
ing  effects  of  dipping  plates  can  cause  explosion  signals  to  appear  highly 
complex.  Overall,  for  explosions  and  earthquakes  in  the  same  area,  complexity 
is  expected  to  be  a  poor  discriminant.  In  fact,  Evernden  (1969)  and  Ericsson 
(1970)  showed  it  to  be  of  little  value  in  comparison  to  Mg  -  m^.  In  studies 
of  short-period  discrimination  which  employ  spectral  ratios  and  complexity 
together  (Anglin,  1971;  Dahlman  et  al.,  1974),  little  of  the  discriminating 
power  can  be  attributed  to  the  complexity  measure.  Greenfield  (1971)  pointed 
out  the  complexity  of  explosion  signals  from  Novaya  Zemlya,  attributing  it  to 
mode  conversion  from  Rayleigh  waves  to  P  waves  in  the  vicinity  of  the  source. 
Clearly,  complexity  caused  in  this  manner  relates  to  the  source  region  and 
not  to  the  nature  of  the  source  and  thus  offers  no  discrimination  capability. 
Complexity,  then,  is  by  far  the  most  tenuous  discriminant. 

Higher-Mode  Surface  Waves 

Shurbet  (1969)  reported  on  the  difference  in  higher-mode  Rayleigh- 
wave  excitation  of  an  explosion  at  NTS  and  a  nearby  earthquake. 


Greenfield,  R.  J.  (1971).  Short-period  P-wave  generation  by  Rayleigh-wave 
scattering  at  Novaya  Zemlya.  J.  Geophys,  Res,,  76,  7988-8002. 
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5339-5341. 
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Forsyth  (1976)  studied  Asian  events  and  found  that  higher  modes  could  aid  in 
identifying  earthquakes.  Observations  such  as  these  could  be  important  for 
discrimination,  but  thus  far  no  detailed  study  has  been  undertaken  for  events 
in  the  United  States.  This  omission  is  due  largely  to  the  higher  detection 
threshold  for  Rayleigh  modes  other  than  the  fundamental  mode  (on  the  order  of 
one  magnitude  unit);  such  thresholds  are  in  agreement  with  theory  (Harkrider, 
1970)  that  predicts  that  the  excitation  of  higher  modes  is  only  a  fraction  of 
the  fundamental  mode  excitation  for  shallow  focal  depths  in  the  Southwestern 
United  States,  as  shown  in  Figure  8. 
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Figure  8  Theoretical  excitation  of  fundamental-mode  and  higher-mode  Rayleigh 
waves  versus  depth  of  focus  in  the  Southwestern  United  States. 
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Figure  8  (cont.)  Theoretical  excitation  of  fundamental-mode  and  higher-mode 
Rayleigh  waves  versus  depth  of  focus  in  the  Southwestern 
United  States. 


Figure  8  (cont.)  Theoretical  excitation  of  fundamental-mode  and  higher-mode 
Rayleigh  waves  versus  depth  of  focus  in  the  Southwestern 
United  States. 


NATURE  OF  THE  CRUST  AND  UPPER  MANTLE 
IN  THE  WESTERN  UNITED  STATES 


General  Geologic  History 

The  geological  history  of  the  Western  United  States  is  as  complex  as  that 
of  any  region  on  the  globe.  For  purposes  of  this  report,  this  history  will  be 
described  only  briefly,  focusing  only  on  those  generally  broad  current  fea¬ 
tures  that  may  be  relevant  to  discrimination  of  events  within  the  region. 

Specific  and  limited  geological  features  that  can  affect  signals  emanating 
from  or  propagating  through  them  will  be  discussed  later  when  selected  events 
are  actually  analyzed. 

A  detailed  account  of  the  region's  geological  history  can  be  found  in 
King  (1969).  The  presence  of  Precambrian  and  Paleozoic  rocks  throughout  the 
area  indicate  that  this  has  long  been  part  of  a  continental  plate,  even  though 
shallow  seas  have  often  covered  it  extensively  in  the  past.  The  most  conspic¬ 
uous  recent  Cenozoic  activities  have  been  a  broad  upwarping  of  the  whole  region, 
major  block  faulting  from  the  Sierra  Nevadas  eastward  through  the  Colorado 
Rockies,  and  strike-slip  faulting  associated  with  the  San  Andreas  transform 
fault  (Hamilton  and  Meyers,  1966).  Figure  9  shows  the  major  physiographic 
regions  of  the  Western  United  States.  The  area  under  study  here  is  roughly 
bounded  by  30°N-42°N  and  104 °W  to  the  California  coast.  In  the  following  sec¬ 
tions  each  of  the  five  major  physiographic  provinces  shown  in  this  figure  are 
discussed  in  terms  of  their  general  static  structure.  A  discussion  of  tectonic 
movements,  seismicity,  and  fault  mechanisms  follows  in  another  section. 

Basin  and  Range 

Knowledge  of  the  present  tectonics  of  the  Basin-and-Range  (B-R)  activity  holds 
the  key  to  understanding  the  surrounding  provinces.  This  province  is  anomalous 
in  its  surface  features,  but  it  is  even  more  so  at  depths  where  the  existence 
of  an  upper  mantle  akin  to  that  under  oceanic  ridges  or  that  behind  island  arcs 
is  firmly  established  by  geophysical  evidence. 


King,  P.  B.,  (1969).  The  tectonics  of  North  America  -  a  discussion  to  accom¬ 
pany  the  tectonic  map  of  North  America.  Geological  Survey  Professional 
Paper  628.  U.S.  Geological  Survey,  Dept,  of  Interior,  Washington,  D.  C. 
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Figure  9  Physiographic  regions  of  the  Southwestern  United  States  and 
schematic  of  tectonic  forces. 


The  crust  is  relatively  thin  for  continental  areas;  refraction  surveys 
show  it  to  be  roughly  30  km  or  less  (Prodehl,  1970;  Pakiser,  1963).  The  crust 
thins  toward  the  Gulf  of  California  to  rougly  20  km  at  the  northern  end  (Thatcher 
and  Brune,  1973).  The  eastern  border  of  the  province  is  fairly  well  defined 
along,  or  just  east  of,  the  Wasatch  Front  where  the  crust  begins  to  thicken 
rapidly  to  40  km  or  greater  beneath  the  Colorado  Plateau  (Braile  et  al. , 

1974).  The  western  border  is  marked  by  the  Sierra  Nevada  orogeny,  where  the 
crust  exceeds  40  km  in  depth  (Oliver  et  al.,  1961;  Mikumo,  1965;  Eaton,  1963). 

The  thinness  of  the  crust,  widespread  block  faulting,  and  crustal  exten¬ 
sion  in  the  Basin  and  Range  province  are  intimately  associated  with  an 
anomalous  upper  mantle  concentrated  beneath  this  region.  The  salient  character¬ 
istics  of  this  anomalous  upper  mantle  are  low  velocities,  high  attenuation, 
high  heat  flow,  and  high  conductivity.  Early  evidence  of  anomalously  low 
velocities  under  the  Basin  and  Range  province  came  from  surface-wave  group- 
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velocity  measurements  of  Press  (1960)  and  Alexander  (1963).  Biswas  and  Knopoff 
(1974)  recently  confirmed  this  using  phase  velocities  of  surface  waves,  and 
their  data  also  showed  that  the  mantle  is  of  nearly  normal  continental  charac¬ 


ter  eastward  of  the  Basin  and  Range  province.  Early  evidence  from  low  vel¬ 


ocities  (Herrin  and  Taggart,  1962)  also  indicated  a  deeper  anomalous  zone  under 
the  Basin  and  Range  province.  Other  important  evidence  for  low  velocities  in 
the  upper  mantle  here  came  from  studies  of  P-wave  travel-time  delays  (Cleary 
and  Hales,  1966;  Herrin  and  Taggart,  1968)  and  S-wave  travel-time  delays  (Hales 
and  Roberts,  1970;  Yasar  and  Nuttli,  1974).  The  P-wave  delays  are  not  so  great 
as  the  S-wave  delays,  but  they  do  indicate  a  low-velocity  zone  for  compres- 
sional  waves  which  is  not  well  established  by  using  surface  waves.  Many  models 
of  the  Basin-Range  upper-mantle  velocities  have  been  derived  from  detailed 
studies  of  P-wave  amplitudes  and  arrival  times  (Archambeau  et  al.,  1969; 
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Julian  and  Anderson,  1968;  Wiggins  and  Helmberger,  1973;  Masse  et  al.,  1972), 
from  apparent  velocities  of  P  (Niazi  and  Anderson,  1965;  Johnson,  1967),  and 
from  apparent  velocities  of  S  (Kovach  and  Robinson,  1969),  These  models  exhi¬ 
bit  significant  differences  in  the  exact  structure  of  the  low-velocity  zone. 
York  and  Helmberger  (1973)  attempted  to  map  the  lateral  structure  by  using  P-PL 
travel  times.  Altogether,  the  data  delineates  a  low-velocity  zone  (LVZ)  for 
both  P  and  S  waves  which  begins  at  or  very  near  the  Moho  (there  is  very  little 
or  no  high-velocity  "lid"  over  the  LVZ)  and  extends  to  a  depth  of  between  150 
and  250  km.  S-wave  velocities  may  be  as  low  as  3.85  km/sec  in  this  zone  (Yasar 
and  Nuttli,  1974),  although  a  range  of  4.0-4. 3  km/sec  is  usually  inferred,  com¬ 
pared  with  normal  upper-mantle  velocities  of  4. 4.-4. 6  km/sec.  P-wave  veloci¬ 
ties  are  not  so  anomalous  as  the  S  velocities,  but  they  are  typically  inferred 
to  be  .2  km/sec  less  than  normal  upper-mantle  P  velocities  under  continents. 
These  velocity  anomalies  abate  considerably  to  the  east  under  the  Colorado 
Plateau  and  where  the  LVZ  become  thinner  and  has  a  thicker  lid  (Archambeau  et 
al.,  1969). 

Other  evidence  of  an  anomalous  zone  in  the  Basin  and  Range  province  is  the 
relatively  high  heat  flow  (Roy  et  al.,  1968;  Sass  et  al. ,  1971;  Warren  et  al., 
1969)  and  relatively  high  electrical  conductivity  (Porath  and  Gough,  1971; 
Gough,  1973),  both  of  which  have  been  associated  roughly  with  the  depths  of 
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the  LVZ.  Gravity  and  magnetic  surveys  show  no  appreciable  anomalies  over 
this  region. 

For  this  study,  the  most  important  feature  is  the  pronounced  low  Q  zone 

associated  with  the  LVZ  in  the  Basin  and  Range  province.  Mapping  of  P-wave 

magnitude  residuals  (Cleary,  1967;  Evernden  and  Clark,  1970)  showed  lower 

amplitudes  in  the  Western  United  States,  but  they  were  not  necessarily  confined 

to  the  Basin  and  Range.  Der  et  al.  (1975),  in  a  study  using  short-period  P  and 

S  waves,  found  high  attenuation  in  the  Western  United  States  relative  to  that  in 

the  Eastern  United  States.  The  same  gross  pattern  is  seen  from  long-period 

•  1 

P  and  S  waves  (Solomon  and  Toksoz,  1970).  Using  S^  waves,  Molnar  and  Oliver 
(1969)  mapped  attenuation  in  the  uppermost  mantle  of  the  Western  United  States 
and  found  that  the  Basin  and  Range  province  has  relatively  high  attenuation. 
Mitchell  (1975)  showed  higher  attenuation  in  the  Western  United  States  using 
surface  waves.  The  values  of  Q  estimated  for  the  upper  mantle  under  the  Basin 
and  Range  province  varied  considerably,  depending  upon  how  much  of  the  actual 
path  lay  in  this  anomalous  region.  Der  and  McElfresh  (1977)  estimated  an 
average  Q  of  100-200  for  the  uppermost  200  km  of  the  mantle  under  the  Western 
United  States,  as  compared  to  roughly  1600  in  the  Eastern  United  States  where 
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a  normal  continental  upper  mantle  exists.  Also  with  short-period  P  waves, 
Veith  and  Clawson  (1972)  postulated  a  Q  of  110  for  a  LVZ  of  100  km  thickness 
while  Helmberger  (1973)  postulated  an  even  lower  Q  on  the  order  of  50  in  the 
LVZ. 

To  account  for  the  presence  of  such  low  Q,  partial  melting  of  the  asthen- 
osphere  has  been  suggested  by,  among  others,  Solomon  (1972a)  and  Archambeau  et 
al.  (1969).  The  cause  of  this  low-Q  character  and  other  associated  phenomena 
(described  above)  has  been  set  forth  in  a  variety  of  theories,  most  of  which 
relate  to  active  tectonics  associated  with  moving  plates.  Menard  (1960),  for 
example,  suggested  that  the  active  oceanic  ridge  in  the  Gulf  of  California 
continues  into  the  B-R  province.  Cook  (1969)  elaborated  on  this  concept, 
pointing  out  the  physiographic  similarities  of  the  B-R  to  oceanic  rifts  and 
the  similarity  of  low  P  velocities  and  high  attenuation  between  the  two. 

Scholz  et  al.  (1971)  thought  this  hypothesis  was  incompatible  with,  among 
other  things,  the  volcanic  pattern  throughout  the  region  and  suggested  instead 
that  the  B-R  province  is  more  akin  to  many  interarc  basins  around  the  Pacific 
(Barazangi  et  al.,  1975)  where  attenuation  is  also  high  and  velocities  low. 

The  crustal  extension  seen  in  the  B-R  province  is  characteristic  of  both  ocean 
rifts  and  interarc  basins.  This  interarc  basin  hypothesis  requires  a  subduc- 
tion  zone.  While  no  convincing  evidence  exists  of  a  contemporary  active  zone 
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west  of  the  B-R  province,  a  prior  active  one,  which  has  been  overridden  by 
westward  movement  of  the  North  American  plate,  has  ample  support  from  a  num¬ 
ber  of  different  workers  (Atwater,  1970;  Hamilton  and  Meyers,  1966;  Shaw  et 
al.,  1971). 

Pacific  Coast  and  Sierra  Nevada 

The  Pacific  Coast  (P-C)  province  is  separated  from  the  B-R  province  by  the 
Sierra  Nevada  orogeny  along  the  eastern  border  of  middle  to  southern  Cali¬ 
fornia.  The  demarcation  becomes  less  clear  physiographically  in  southern 
California,  but  King  (1969)  has  asserted  that  the  Sierra  Nevada  orogeny 
extends  down  into  the  Baja  Peninsula,  thus  providing  the  necessary  division. 

(To  the  north  the  division  with  the  B-R  province  becomes  even  less  clear,  but 
this  area  is  beyond  our  study.)  The  Sierra  Nevada  has  a  crustal  "root"  down 
to  at  least  40  km.  In  contrast  the  crustal  structure  of  the  Pacific  Coast 
is  generally  shallow,  about  20-25  km  (Press,  1960;  Mikumo,  1965;  Healy,  1963). 
This  thickness  is  slightly  thinner  than  that  of  the  B-R  crust;  but  it  does  not 
have  the  low  P^  velocities  of  the  B-R  province  (Herrin  and  Taggart,  1962;  Paki- 
ser,  1963),  indicating  that  at  least  the  lid  over  the  LVZ  has  been  restored  in 
this  region  relative  to  the  B-R.  Evidence  for  a  LVZ  under  the  Pacific  Coast 
province  is  not  so  complete  as  for  the  B-R  province,  but  the  high  attenuation 
inferred  from  short-period  waves,  as  referenced  above  for  the  B-R  province, 
apparently  continues  nearly  unabated  into  the  P-C  province.  However,  P-wave 
delay  times,  as  referenced  above  for  the  B-R  province,  disappear  westward 
toward  the  coast. 
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Colorado  Plateau 

At  the  Wasatch  Front,  the  crust  begins  to  thicken  eastward  (Keller  et  al., 
1975)  to  values  of  40-45  km  for  the  Colorado  Plateau  (Prodehl,  1970).  The 
Wasatch  Front,  however,  is  not  the  boundary  for  high  attenuation,  high  heat 
flow,  or  high  conductivity,  because  these  characteristics  persist  eastward 
under  the  Colorado  Plateau  and  into  the  southern  Rockies.  The  upper-mantle 
velocity  profile  for  the  Colorado  Plateau  is  distinct  from  the  B-R  province 
because  P^  velocities  are  normal  (Herrin  and  Taggart,  1962)  and  there  is  a 
thicker  lid  and  less  pronounced  velocity  anomalies  in  the  upper  mantle  (Archam- 
beau  et  al.,  1969;  Biswas  and  Knopoff,  1974;  York  and  Helmberger,  1973). 

Whether  those  anomalous  features  are  a  remnant  of  past  orogenic  activity  (King, 
1969),  the  first  westward  indications  of  the  overriding  of  an  oceanic  rift 
(Cook,  1969),  or  the  outermost  effects  of  the  formation  of  the  ensialic  inter¬ 
arc  basin  (Scholz  et  al.,  1971)  remains  an  open  question.  The  area  as  a  whole 
is  relatively  aseismic  (Smith  and  Sbar,  1974),  a  fact  that  further  attests  to 
its  stability  and  cohesiveness. 

Southern  Rockies 

This  province  is  the  easternmost  border  of  the  tectonic  cycles  that  have 
churned  the  Western  United  States  and,  as  such,  it  is  the  most  eastward  expres¬ 
sion  of  the  anomalous  upper  mantle  that  characterizes  provinces  to  the  west. 

The  crustal  thickness  is  about  50  km  (Pakiser,  1963),  which  is  only  slightly 
greater  than  for  the  Colorado  Plateau.  The  anomalous  geophysical  aspects  of 
high  conductivity  and  heat  flow  seem  to  persist  eastward  from  the  B-R  province 
into  this  province  (Reitzel  et  al.,  1970;  Porath  and  Gough,  1971).  The  nature 
of  the  LVZ  under  this  province  is  not  well  defined  although  it  is  certainly 
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extant  to  some  degree  in  the  southern  portion,  which  is  a  northward  extention 
of  the  Rio  Grande  rift  system  with  its  anomalous  upper-mantle  velocities 
(York  and  Helmberger,  1973). 
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EARTHQUAKES  IN  THE  SOUTHWESTERN  UNITED  STATES 


Tectonic  Forces 

The  development  of  plate  tectonics  as  the  unifying  theory  for  interpreting 
the  global  distribution  of  seismicity  in  relation  to  prominent  physiographic 
features  within  and  around  the  ocean  margins  (Isacks  et  al.,  1968)  has  gener¬ 
ated  reinterpretations  of  the  present  structure  and  activity  of  the  Western 
United  States.  At  least  four  relatively  distinct  interpretations  of  the  ano¬ 
malous  zone  centered  under  the  B-R  province,  all  associated  with  plate  tecton¬ 
ics,  have  been  recently  presented:  1)  Cook  (1969)  set  forth  the  rift  theory 
holding  that  the  North  American  plate  has  overridden  the  East  Pacific  rise 
which  now  supposedly  lies  beneath  the  Basin  and  Range  province;  2)  Atwater 
(1970)  argued  that  the  Basin  and  Range  province  is  a  result  of  broad  right- 
lateral  movement  in  the  Western  United  States,  centered  in  the  San  Andreas  and 
associated  faults,  but  not  entirely  accommodated  there;  3)  Shaw  et  al.  (1971) 
suggested  that  earth  tidal  power  generated  a  thermal  source  behind  an  subduc- 
tion  zone  on  the  West  Coast,  which  has  now  been  assimilated,  leaving  the  rem¬ 
nant  thermal  source  under  the  Basin  and  Range  province;  4)  Scholz  et  al. 

(1971)  claimed  that  release  of  compressive  stress,  due  to  the  termination  of 
an  active  subduction  zone  on  the  West  Coast,  enabled  a  mantle  diapir  created 
by  the  subduction  zone  to  initiate  crustal  extension  in  the  Basin  and  Range 
province  in  a  way  similar  to  several  interarc  basins  in  the  Pacific  Ocean.  All 
these  hypotheses  are  satisfactory  in  explaining  observed  seismicity  of  the 
Southwestern  United  States  because  each  recognizes  the  two  dominant  patterns 
shown  on  Figure  9  —  transform  motion  along  the  San  Andreas  fault  zone  and 
crustal  extension  in  the  Basin  and  Range  province.  The  presentation  of  Scholz 
et  al.  (1971),  as  well  as  that  of  Smith  and  Sbar  (1974)  and  Bolt  et  al.  (1968), 
gives  a  particularly  thorough  interpretation  of  present  seismic  activity  in  the 
region  of  this  study. 

Seismicity 

Figure  10  shows  the  distribution  of  the  seismic  activity  in  the  South¬ 
western  United  States  over  the  past  fifteen  years,  as  recorded  in  the  NEIS 
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epicenter  file.  All  earthquakes  with  m^  >  3.5  have  been  plotted,  except  with¬ 
in  the  Nevada  Test  Site,  where  most  earthquakes  are  known  to  be  associated 
with  explosions  (Hamilton  et  al.,  1972).  Explosions  at  NTS  have  also  been 
eliminated,  along  with  the  few  other  scattered  events  known  to  be  explosions 
in  the  Southwestern  United  States.  Three  concentrated  seismic  zones  appear 
to  exist.  The  first  zone  is  the  NW-SE  lineation  associated  with  the  trans¬ 
form  motion  of  the  San  Andreas  fault  in  western  California.  A  second  occurs 
in  eastern  California  and  western  Nevada.  A  third  prominent  zone  is  the 
interraountain  seismic  belt  running  N-S  along  the  Wasatch  Front  in  Utah. 

The  first  zone  along  western  California  contains  events  with  primarily 
right-lateral  strike-slip  motion  (Allen  et  al. ,  1965;  Bolt  et  al.,  1968). 

This  observation  is  consistent  with  Wilson's  (1965)  hypothesis  that  the  San 
Andreas  system  is  a  transform  fault  linking  the  East  Pacific  Rise  and  the 
Gorda  Ridge.  This  motion  reflects  an  average  slip  rate  of  roughly  4  mm/year 
(Bolt  et  al.,  1968),  which  is  compatible  with  geodetic  determinations  (Sav¬ 
age  and  Burford,  1973)  and  with  known  spreading  rates  for  the  East  Pacific 
Rise  and  Gorda  Rise.  The  only  major  deviation  from  this  right- lateral  strike- 
slip  trend  is  along  the  Transverse  Ranges  and  the  Garlock  Fault  (see  Figure 
9),  where  thrust  and  steep  reverse  faults  with  a  left-lateral  component  have 
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been  identified  (Alien  et  al.,  1965).  For  southern  California,  Thatcher 
and  Hanks  (1973)  and  Wyss  and  Brune  (1971)  provided  evidence  for  a  consider¬ 
able  range  of  source  mechanisms,  dimensions,  and  stress  so  that  a  generali¬ 
zation  about  this  area  is  not  possible  though  localized  similarities  do  occur. 
Northern  Baja  California  and  the  northern  Gulf  of  California  are  similarly 
varied.  However,  Thatcher  (1972)  pointed  out  gross  differences  between  most 
Baja  and  Gulf  earthquakes  in  this  area.  The  differences  are  smaller  source 
dimensions  and  moments  for  the  Baja  events  and  larger  stress  drops  for  them, 
one  or  two  orders  of  magnitude  perhaps.  Thatcher's  work  shows  that  several 
of  these  Baja  earthquakes  are  difficult  to  discriminate  by  M^  -  m^  from  NTS 
explosions.  It  also  shows  that  there  are  a  number  of  exceptions  where  Baja 
events  have  large  source  dimensions  like  the  Gulf  events. 

The  second  zone  of  concentrated  seismicity,  the  "Nevada  seismic  zone," 
reflects  both  the  crustal  extension  produced  by  the  anomalous  upper  mantle 
under  the  B-R  province  and  the  broad  right-lateral  plate  motion  of  the  North 
American  plate  against  the  Pacific  plate.  This  motion  interacts  with  the 
fairly  rigid  Sierra  Nevada  batholith  (Gumper  and  Scholz,  1971)  producing  a 
varied  stress  pattern.  Nearly  all  earthquake  mechanisms  east  and  north  of 
this  seismic  zone  are  of  dip-slip  nature  (Smith  and  Sbar,  1974;  Gumper  and 
Scholz,  1971;  Ryall  and  Malone,  1971;  Stauder  and  Ryall,  1976)  which  is  com¬ 
patible  with  simple  block  faulting  in  the  B-R  province.  The  crustal  extension 
required  for  this  phenomenon  is  affirmed  by  strain  and  geodetic  measurements 
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(Meister  et  al. ,  1968;  Priestly,  1974).  Overall,  focal  mechanisms  and  these 
physical  measurements  suggest  that  the  extension  is  in  the  NW-SE  direction, 
with  the  E-W  component  becoming  progressively  larger  eastward.  The  total 
extension  in  the  late  Cenozoic  era  is  at  least  100  km  and,  according  to  Scholz 
et  al.  (1971),  it  may  be  as  much  as  300  km. 

In  the  Nevada  Seismic  Zone,  but  toward  the  California-Nevada  border,  more 
complex  tectonics  exist.  Gumper  and  Scholz  (1971)  determined  that  in  addi¬ 
tion  to  normal  faulting,  left-lateral  strike-slip  motion  existed  over  a  broad 
area  near  38°  N,  118°-119°  W.  Gumper  and  Scholz  asserted  that  this  area  was 
a  transform  between  two  lines  of  crustal  extension  to  the  north  and  south  of 
roughly  38°N,  with  the  northern  seismicity  offset  to  the  east.  In  addition, 
Gumper  and  Scholz  found  a  right-lateral  component  to  several  of  the  predomi¬ 
nant  dip-slip  mechanisms  of  this  area.  Although  many  fault  mechanisms  have 
been  determined  in  this  area,  data  on  stress  drops  and  source  dimensions  for 
earthquakes  in  this  zone  are  sparse. 

The  third  concentrated  zone  in  the  region  under  study  is  the  "intermoun¬ 
tain  seismic  belt,"  described  in  detail  by  Smith  and  Sbar  (1974)  and  Sbar  et 
al.  (1972).  This  zone  encompasses  the  physiographic  boundary  known  as  the 
Wasatch  Front  in  central  Utah  and  continues  northward  beyond  the  region  under 
study.  Dip-slip  motion  on  steeply-dipping  fault  planes  is  indicated  by  com¬ 
posite  focal-mechanism  solutions  for  various  parts  of  this  zone.  Smith  and 
Sbar  assert  that  this  can  be  explained  either  as  differential  vertical  move¬ 
ment  of  the  B-R  province  relative  to  the  Colorado  Plateau,  a  manifestation  of 
crustal  extension  characteristic  of  the  B-R  province  to  the  west,  or  as  evi¬ 
dence  that  the  B-R  province  is  a  subplate  moving  westward  relative  to  North 
America.  They  have  also  correlated  a  number  of  earthquake  swarms  in  this 
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zone  with  areas  of  high  heat  flow  and  characteristic  geothermal  features. 

The  westward-trending  lower  part  of  this  zone  in  southern  Utah  and  eastern 
Nevada  appears  to  have  left-lateral  strike-slip  motion,  which  is  further  evi¬ 
dence  of  westward  motion  of  the  B-R  province  relative  to  the  Colorado  Pla¬ 
teau.  Almost  no  data  is  available  on  source  dimensions  and  stress  drops  for 
earthquakes  in  this  zone. 

In  addition  to  the  natural  seismicity  in  the  Western  United  States,  exam¬ 
ples  exist  of  artificial  seismicity.  The  artificial  seismicity  associated 
with  NTS  shots  is  well  known.  Also,  earthquakes  have  been  associated  with 
fluid  injection  near  Denver,  Colorado  (described  by  Healy  et  al.,  1968). 
Effective  stress  and  drops  for  these  events  have  been  determined  by  Wyss  and 
Molnar  (1972).  Their  calculations  suggest  small  stress  drops  and  normal 
source  dimensions  relative  to  values  expected  for  their  magnitude  (M^) . 

Depth  of  Focus 

Many  authors  (e.g.,  Brace  and  Byerlee,  1970;  Cook  and  Smith,  1967;  Tocher, 
1958)  have  affirmed  that  all  seismicity  in  the  Southwestern  United  States  is 
shallow.  All  accurately  located  earthquakes  have  been  within  the  crust,  the 
majority  at  depths  of  less  than  15  km.  No  present  seismic  evidence  exists 
for  continued  motion  of  the  subduction  zone,  which  has  been  postulated  to  be 
active  on  the  West  Coast  during  the  Cenozoic  Era.  In  the  B-R  province,  Smith 
and  Sbar  (1974)  reported  a  rather  abrupt  cutoff  in  seismicity  related  to  a  crus¬ 
tal  low-velocity  zone  starting  at  depth  of  roughly  15  km  (Braile  et  al.,  1974; 


Healy,  J.  H.,  W.  W.  Rubey,  D.  T.  Griggs,  and  C.  B.  Rayleigh  (1968).  The  Den¬ 
ver  earthquakes.  Science,  161 ,  1301-1310 

Wyss,  M. ,  and  P.  Molnar  (1972).  Efficiency,  stress  drop,  apparent  stress, 

effective  stress,  and  frictional  stress  of  Denver,  Colorado,-  earthquakes. 
J.  Geophys.  Res.,  77,  1433-1438. 

Brace,  W.  F. ,  and  J.  D.  Byerlee  (1970).  California  earthquakes,  why  only  shal¬ 
low  focus.  Science .  156,  1573-1575. 

Cook,  K.  L.,  and  R.  B.  Smith  (1967).  Seismicity  in  Utah,  1850  through  June 
1965.  Bull.  Seism.  Soc,  Am.,  57,  689-718. 

Tocher,  D.  (1958).  Earthquake  energy  and  ground  breakage.  Bull  Seism.  Soc. 
Am.,  48,  147-153. 


Mueller  and  Landisman,  1971).  Shurbet  and  Cebull  (1971)  elaborated  on  this 
idea,  postulating  that  such  a  crustal  layer  is  undergoing  extension  and  that 
it  is  the  origin  of  the  graben-and-horst  character  of  the  B-R  province.  It 
is  most  reasonable  that  the  extension  be  due  to  an  anomalous  upper  mantle, 
and  that,  if  seismicity  indeed  abates  below  15  km,  then  the  lower  crust  must 
be  deformed  plastically  along  with  the  upper  mantle. 

In  sum,  the  shallowness  of  seismicity  in  the  Western  United  States  means 
that  location  alone  is  of  little  importance  in  discriminating  events  in  this 
region  unless  master  events  were  widely  employed  to  narrow  the  depth  confi¬ 
dence  limits  produced  by  a  high-quality  network  of  teleseismic  stations.  How¬ 
ever,  use  of  master  events  may  not  be  the  solution  because  Blandford  (1975) 
found  a  20-km  standard  deviation  for  depth  estimates  of  master-controlled 
NTS  explosions  when  well-distributed  networks  of  teleseismic  stations  were 
used. 


Mueller,  S.,  and  M.  Landisman  (1971).  An  example  of  the  unified  method  of 
interpretation  for  crustal  seismic  data.  Geophys,  J.,  23,  365-371. 

Shurbet,  D.  H.,  and  S.  E.  Cebull  (1971).  Crustal  low-velocity  layer  and 
regional  extension  in  the  Basin  and  Range  province.  Bull.  Geol.  Soc. 
Am.,  82,  3241-3244. 

Blandford,  R.  R.  (1975).  Use  of  source-region-station-time  corrections  at 
NTS  for  depth  estimation.  Report  No.  SDAC-TR-75-4,  Teledyne  Geotech, 
Alexandria,  Virginia. 
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The  twenty-three  seismic  events  described  below  were  selected  from  known 
events  in  the  Southwestern  United  States.  Eleven  were  explosions  and  twelve 
were  earthquakes.  This  sample  was  considered  sufficient  to  produce  clear 
results  and  also  to  supply  many  important  and  interesting  particular  features 
for  study.  In  order  to  take  advantage  of  the  continuous  recording  of  several 
sites  in  the  Long  Range  Seismic  Measurements  (LRSM)  network  of  the  VELA-Uniform 
program,  the  events  were  restricted  to  the  years  1966-1969.  The  magnitude 
range,  4.8  to  6.5  m^,  of  these  events  insures  that  most  signals  recorded  by  the 
North  American  LRSM  network  are  of  good  signal/noise  ratio.  Low-magnitude 
events  were  intentionally  rejected  since  this  study  is  not  one  of  detection 
and  signal  enhancement  but  one  of  multi-dimensional  discrimination  based  upon 
the  best  available  data.  However,  results  attained  here  for  medium-magnitude 
events  are  expected  to  apply  in  most  aspects  to  events  of  lower  magnitude. 

The  only  significant  deviations  at  lower  magnitudes  would  result  from  dif¬ 
ferent  scaling  of  source  dimensions  and  time  functions  between  earthquakes 
and  explosions  for  lower  magnitudes  versus  medium  magnitudes. 

Data  for  the  eleven  explosions  listed  in  Table  III  were  taken  from 
Springer  and  Kinnaman  (1971).  Eight  of  these  explosions  were  at  the  Nevada 
Test  Site  and  represented  various  yields,  depths,  locations,  and  mediums. 
Several,  such  as  PILE  DRIVER,  GREELEY,  and  BENHAM,  had  a  considerable  com¬ 
ponent  of  tectonic  strain  release  (Toksoz  and  Kehrer,  1972a);  these  explosions 
were  deliberately  included  to  empirically  assess  the  degree  of  adverse  effect 
on  discrimination  entailed  in  this  phenomenon.  Only  three  nuclear  explosions — 
FAULTLESS,  GASBUGGY,  and  RULISON —  were  detonated  outside  NTS  in  the  desired 
time  period;  they  were  included  to  spread  the  explosion  sources  over  a  wide 
area  of  the  Southwestern  United  States.  The  locations  of  all  these  explosions 
have  been  plotted  on  Figures  9  and  11. 


Springer,  D.  L. ,  and  R.  L.  Kinnaman  (1971).  Seismic  source  summary  for 
U.  S.  underground  nuclear  explosions,  1961-1970.  Bull.  Seism.  Soc., 
Am.,  61,  1073-1098. 
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The  twelve  earthquakes  selected  are  also  listed  in  Table  III;  the  epicen¬ 
ter  data  is  from  the  NEIS  list.  The  symbols  associated  with  the  earthquakes 
and  explosions  in  this  Table  are  used  throughout  this  report.  As  seen  in 

Figures  9  and  11,  the  earthquakes  are  reasonably  well  distributed  over  the 

region  of  study.  Four  of  these  earthquakes  are  in  Nevada  and  will  provide 
the  most  direct  comparison  with  explosions.  Each  of  these  events  have  them¬ 
selves  been  the  subject  of  separate  studies  or  part  of  a  group  of  events 
studied  or  examined  in  another  published  work.  Their  presence  in  such  studies 

was  one  criterion  for  inclusion  in  this  study.  Another  criterion  used  for 

earthquake  selection  was  that  they  should  only  slightly  separate  from  NTS 
explosions  in  published  Mg  -  m^  plots.  All  of  this  previous  information  fits 
into  this  regional  study  and  aids  in  sustaining  its  conclusions. 

A  subset  of  the  LRSM  network,  namely  WH-YK,  NP-NT,  RK-ON,  SV3QB,  HN-ME, 
BE-FL,  and  PG-BC,  plotted  in  Figure  11,  was  chosen  for  this  study.  Each  sta¬ 
tion  of  this  subset  was  operational  for  most,  or  all,  the  time  encompassing 
the  selected  events.  The  stations  were  chosen  because  they  are  teleseismic 
or  nearly  teleseismic  from  the  selected  events.  This  situation  eases  prob¬ 
lems  associated  with  multiple  arrivals  due  to  travel-time  triplications  or 
problems  inherent  in  head-wave  interpretations.  Teleseismic  distances 
tend  to  simplify  magnitude  determination  compared  to  regional  distances. 
Although  this  selected  subset  provides  nearly  as  wide  an  azimuthal  distri¬ 
bution  of  teleseismic  stations  as  possible  over  North  America,  relative  to 
events  in  the  Southwestern  United  States,  the  coverage  still  is  barely  more 
than  one  quadrant.  All  the  seismograms  were  re-examined  to  produce  the 
amplitude  data  for  this  report.  Other  data,  in  the  form  of  spectral  ratios 
and  complexities,  was  obtained  from  digitized  recordings.  The  data  base 
pertinent  to  discrimination  of  these  events  is  listed  in  Table  IV.  The  mag¬ 
nitude  data  of  this  list  was  computed  according  to  the  common  relations 

=  log (A/T)  +  B 

Mg  =  log (A/T)  +-1.66  logA  +0.3 
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TABU  IV  (Continued) 


A  -  one-half  the  peak-to-trough  maximum  amplitude  on  the  recording,  re¬ 
duced  to  run, 

T  -  period  in  seconds, 

A  -  distance  in  degrees, 

B  -  Gutenberg-Riehter  distance-correction  terms  for  P  and  S  phases. 

For  the  event  magnitudes  to  be  shown  in  later  plots,  the  magnitude  esti¬ 
mation  technique  of  Kingdal  (1976)  has  been  employed.  This  technique  tends 
to  eliminate  upward  bias,  relative  to  the  true  magnitude,  of  an  event  magni¬ 
tude  which  is  averaged  over  only  a  few  recording  stations  for  small  events. 

The  technique  will  be  explained  in  detail  in  a  later  section.  If  no  signals 
of  a  particular  type  were  recorded  for  an  event,  then  an  event  magnitude  will 
be  plotted  as  the  average  of  the  station  magnitudes  based  on  noise  measurements. 
This  approach  tends  to  overestimate  the  true  magnitude,  but  there  is  no  sim¬ 
ple  and  satisfactory  alternative. 

The  complexity  values  in  Table  IV  were  computed  in  the  manner  shown  by 
Lambert  et  al.  (1969).  Essentially,  they  represent  the  inverse  of  the  power 
in  the  first  5  seconds  of  P  arrival  divided  by  that  in  the  following  30  sec¬ 
onds  of  coda  (higher  complexity  value  means  more  coda).  The  spectral  ratios 
were  computed  on  the  raw  spectrum  of  the  P  and  LR  signals.  For  P  waves,  the 
first  6.4  seconds  of  signal  data,  sampled  at  20  pts/sec,  were  tapered  and 
transformed.  For  LR  waves,  the  data  in  the  group-velocity  window  from  roughly 
3.7  to  2.8  km/sec,  sampled  at  1  pt/sec,  were  tapered  and  transformed.  The 
P-wave  spectral  ratio  was  computed  as  the  sum  of  the  Fourier  amplitude  coef¬ 
ficients  from  1.56  to  1.87  Hz  over  that  from  0.47  to  0.78  Hz,  and  the  LR-wave 
ratio  used  0.0459  to  0.0718  Hz  over  0.0283  to  0.0449  Hz.  In  both  cases,  the 
spectral  ratio  represents  high  over  low  frequencies.  Spectral  ratios  regarded 
as  inaccurate  due  to  interfering  noise  in  any  of  these  bands,  are  indicated 
in  Table  IV. 

The  focal  mechanisms  of  the  selected  events  are  plotted  in  Figure  12, 
and  the  source  of  each  focal  mechanism  determination  is  listed  in  Table  V. 
Several  earthquakes  had  more  than  one  proposed  mechanism,  and  the  mechanism 
that  seemed  best  was  chosen.  Note  that  some  published  mechanisms  adopted 
here  rested  on  scant  data  and  could  be  significantly  in  error.  Three  events 
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Figure  12  Focal  mechanism  of  the  earthquakes  selected  for  this  study 
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Figure  12  (cont.)  Focal  mechanism  of  the  earthquakes  selected  for  this  study 
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Figure  12  (cont.)  Focal  mechanism  of  the  earthquakes  selected  for  this  study 
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Figure  12  (cont.)  Focal  mechanism  of  the  earthquakes  selected  for  this  study 


TABLE  V 


Sources  of  Information  on  the  Focal  Mechanisms 
of  Earthquakes  Selected  for  this  Study 

Name  Reference 


Dulce,  N.M.  EQ 
Caliente  EQ 
Truckee  EQ 
Caliente  AS 
Denver  EQ 
Ba j a  EQ 

Utah  EQ 

Borrego  Mtn.  EQ 
Nevada  E0 

BENHAM  AS 

29  Palms  EQ 


Coyote  Mtn.  EQ 


Smith  and  Sbar  (1974) 

Tsai  and  Aki  (1971) 

Tsai  and  Aki  (1970) 
assumed  to  be  same  event  B 
Major  and  Simon  (1968) 
assumed  or  basis  of  nearby 
fault  motion 
Smith  and  Sbar  (1974) 

Hanks  and  Wyss  (1972) 
assumed  on  the  basis  of  nearby 
fault  mechanism  solutions 
Savino  et  al.  (1971) 
assumed  on  the  basis  of  fault 
orientation  and  ISC  first 
motions 

Thatcher  and  Hamilton  (1973) 


have  no  published  mechanisms,  and  the  mechanics  of  faulting  were  inferred 
from  association  of  the  events  with  known  faults  or  proximity  to  earthquakes 
with  known  mechanisms.  First-motion  data  plotted  in  these  figures  is  from 
our  analysis  of  the  seven  LRSM  stations  and  will  be  discussed  later  in  this 
report. 
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SOURCE  EFFECTS 


In  this  section  aspects  of  the  source  responsible  for  identifiable  dif¬ 
ferences  in  the  recorded  signals  of  our  selected  events  are  discussed  within 
a  context  of  discrimination  which  employs  commonly  regarded  measures  such  as 
Mg  -  m^,  spectral  ratios,  shear-wave  generation,  and  complexity.  The  variety 
of  tectonic  processes  occurring  within  the  source  region  implies  that  the 
selected  earthquakes  cannot  be  characterized  uniformly.  In  regard  to  the 
explosions,  a  more  nearly  uniform  source  behavior  is  anticipated. 

Polarity  of  First  Motion 

In  the  focal-mechanism  diagrams  of  Figure  12,  data  from  the  seven  LRSM 
stations  are  added.  Except  for  larger  events,  where  long-period  P  waves  were 
often  recorded,  this  first  motion  data  is  from  short-period  instruments.  In 
most  cases,  this  data  is  consistent  with  the  known  or  assumed  fault  mechanisms, 
and  the  exceptions  merely  demonstrate  the  well-known  inadequacy  of  first  motion 
taken  from  short-period  recordings.  Because  of  epicenter-to-station 
distance  and  because  of  the  signal-to-noise  ratio,  few  of  the  observed 
first  motions  would  have  satisfied  the  criteria  set  forth  by  the  Technical 
Working  Group  II  of  the  Geneva  Conference  (U.S.  Joint  Committee  on  Atomic 
Energy,  1960).  For  ten  of  the  twelve  earthquakes,  at  least  one  apparent  dilata- 
tional  first  motion  was  found  among  the  stations;  such  data  cannot  constitute 
positive  discrimination  of  source  type,  but  merely  a  strong  diagnostic. 
Furthermore,  an  examination  of  P  signals  on  these  LRSM  recordings  of  the 
eleven  explosions  revealed  several  P  signals  that  could  be  described  objectively 
as  dilatational  first  motions. 

S  Waves 

Earlier  work  of  von  Seggern  (1972)  already  demonstrated  with  real  data 
that  measurement  of  the  excitation  of  S  waves  relative  to  P  waves  should  pre¬ 
sent  a  reliable  discriminant.  In  Figures  13  and  14  are  short-period  S  and 
long-period  S  magnitudes  plotted  against  routine  Many  of  the  events  are 


U.S.  Joint  Committee  on  Atomic  Energy  (1960).  Technical  aspects  of  detection 
and  inspection  controls  of  a  nuclear  weapon  test  ban  (Summary  analysis  of 
hearings  before  the  Special  Subcommittee  on  Radiation  and  the  Subcommittee 
on  Research  and  Development,  April  19-22,  1960).  Superintendent  of 
Documents,  U.S.  Government  Printing  Office,  Washington,  D.C. 
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Short-period  S  magnitude  versus  short-period  P  magnitude  for 
selected  Southwestern  United  States  events. 
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represented  by  nois.e  estimates  (see  Table  IV  and  VI),  and  the  effective  threshold 
for  application  of  this  discriminant  is  m^  ^5.0  in  the  short-period  mode  and  m^ 

^  4.5  in  the  long-period  mode.  These  figures  show  that,  when  measurable,  the 

short-period  and  long-period  S  waves  are  generally  smaller  relative  to  P  for 
the  explosions  than  the  earthquakes.  There  are,  however,  several  exceptions: 
earthquakes  A  and  B  have  low  short-period  S  and  event  L,  which  is  difficult 
to  identify  on  the  basis  of  Mg  -  m^  also,  is  not  well  identified  on  the  basis 
of  S  waves.  The  Coyote  Mountain  earthquake,  event  L,  is  especially  trouble¬ 
some  because,  while  its  short-period  S  magnitude  in  Figure  13  is  based 
upon  noise  measurements,  the  lack  of  visible  S  waves  on  short-period  components 
for  a  ^  6  earthquake  is  certainly  anomalous.  Von  Seggern  (1972)  found  no 
S/P  amplitude  ratios  <  0.1  on  LRSM  recordings  of  ii^  >  6.0  earthquakes  world¬ 
wide;  undetected  S  waves  were  not  treated  in  von  Seggern' s  study  though. 


Excitation  of  LQ 

Figure  15  outlines  relative  LQ/LR  excitation  of  the  selected  events  and 
shows  that  no  distinct  separation  between  the  earthquakes  and  explosions  exists. 
PILE  DRIVER  (2)  and  GREELEY  (4)  are  the  events  Toksoz  an(j  Kehrer  (1972a)  assigned 
the  largest  tectonic  component  in  a  study  of  strain  release  by  underground  explo¬ 
sions  at  the  Nevada  Test  Site.  Note  that  these  two  explosions  were  not  included 
in  the  earthquake  population  in  terms  of  short-period  and  long-period  S  excita¬ 
tion  (Figures  13  and  14.)  It  is  postulated  that  the  long-period  S-wave  ampli¬ 
tude,  which  is  due  to  tectonic-strain  release,  is  also  cancelled  by  the  sur¬ 
face  reflection  in  the  same  way  as  is  shallow  explosion  P  by  np.  In  addition,  the 
authors  think  that  the  tectonic-strain  release  component  is  deficient  in 
short-period  S  waves  due  to  a  long  time  or  space  dimension  for  the  source. 

Neither  of  these  arguments  prevents  the  degree  of  LQ  excitation  required  to 
explain  transverse-component  recordings  of  explosions.  RULISON,  Event  11,  in 
the  more  stable  Colorado  Plateau,  has  appreciably  less  LQ  than  Nevada  explo¬ 
sions;  and  GASBUGGY,  Event  7,  which  lies  close  to  seismic  activity  in  north¬ 
ern  New  Mexico,  appears  to  excite  LQ  to  the  same  degree,  that  is,  lower  than 
at  NTS.  The  low  LQ  amplitudes  of  RULISON  and  GASBUGGY  provide  further  evidence 
that  the  high  LQ  excitation  of  certain  NTS  explosions  is  a  source-related 
phenomenon  and  is  not  due  to  propagation-path  effects,  such  as  mode  conversion. 
Overall,  excitation  of  shear  modes  as  represented  by  short-period  S,  long- 


-110- 


3  0 


4  0 


SO 


ID 


7  0 


10 


LO 


Figure  15  Rayleigh-wave  magnitude  versus  Love-wave  magnitude  tor  selec 
Southwestern  United  States  events. 
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period  S,  and  LQ,  relative  to  P  and  LR,  betrays  whether  the  source  type  is  double 
couple  or  uniform  compression.  Note  that  when  attempts  are  made  to  confuse  iden¬ 
tification  of  explosions  by  the  firing  of  shot  arrays  to  enhance  Mg  without 
similarly  increasing  m^  and  to  produce  pP-like  phases  (Kolar  and  Pruvost,  1975), 
then  the  low  relative  excitation  of  S  waves  and  LQ  waves  by  an  explosion  will 
not  be  changed,  and  careful  analysis  of  recordings  would  suggest  the  possibil¬ 
ity  of  such  an  evasion  attempt. 

P-Wave  Spectrum  and  Frequency 

Figure  16  shows  excitation  of  long-period  P  waves  relative  to  short-per¬ 
iod  P  waves,  as  measured  on  their  respective  LRSM  components.  A  fairly  clear 
separation  exists  between  the  two  classes  of  events  whenever  long-period  P 
was  detected.  Event  L,  the  Coyote  Mountain  earthquake,  is  somewhat  low  and 
also  appeared  as  rather  inefficient  in  long-period  S  excitation  in  Figure  14. 

The  known  depth  of  this  earthquake,  10-15  km,  precludes  pP  cancellation  of  the 
longer-period  portion  of  the  P-wave  spectrum  (Mo.lnar,  1971;  Helmberger  and 
Harkrider,  1972),  which  is  evident  for  explosions  in  the  data  set.  The  explo¬ 
sions  with  highest  relative  long-period  P  amplitude  have  the  largest  yields 
in  the  group,  and  the  relative  increase  of  long-period  P  is  due  to  the  decrease 
in  corner  frequency  with  yield.  This  causes  from  the  short-period  P  to  be 
measured  at  frequencies  beyond  this  corner  in  the  spectrum  and  thus  at  a  lower 
level  than  that  of  a  smaller  explosion. 

Study  of  the  spectra  of  the  P  waves  provides  a  fuller  understanding  of 
their  generation.  The  P-wave  spectra  of  the  selected  events  from  recordings 
at  NP-NT  and  RK-ON  are  shown  in  Figures  17  and  18,  respectively.  More  record¬ 
ings  were  recovered  for  these  stations  than  for  the  other  stations  in  our 
group,  so  comparative  analysis  was  facilitated  for  them.  As  shown  in  Figure  11, 
NP-NT  is  teleseismic  and  somewhat  equidistant  to  all  the  events  while  RK-ON  has  a 
large  range  of  epicentral  distances,  some  of  which  entail  multiple  arrivals. 

The  good  S/N  ratios  at  RK-ON  for  high  frequencies  is  an  advantage,  however. 

To  estimate  the  source  spectrum  Us(w),  corrections  have  been  made  to  the 
observed  spectrum  Uo(ai)  according  to 

U  (w)  -  U  (m)e1Tft*/I(w) 
s  o 

Kolar,  0.  C. ,  and  N.  L.  Pruvost  (1975).  Earthquake  simulation  by  nuclear 

explosions.  Nature,  253,  242-245. 
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Figure  17  Spectra  of  P  waves  from  short-period  recordings 
selected  Southwestern  United  States  events. 
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Figure  17  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  NP-NT 

of  selected  Southwestern  United  States  events. 
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Figure  17  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  NP-NT 

of  selected  Southwestern  United  States  events. 


Figure  17  (cont.)  Spectra  of  P  waves  from  short-period  recordings 

of  selected  Southwestern  United  States  events. 
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Figure  17  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  NP-NT 

of  selected  Southwestern  United  States  events. 
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Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 
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Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 
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Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 
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Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 
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Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 
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Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 


HALFBEAK 


o  I 

v:  m 

h 


ot  o'e  o*i  o 
C  S0NDD3£-H3l3UDNdN  mi 


"ItiNSIS 


3SI0N 


-149- 


i 


Figure 


NOISE  AMPLITUDE  SPECTRUM 


RULISON 


Figure  18  (cont.)  Spectra  of  P  waves  from  short-period  recordings  at  RK-ON 

of  selected  Southwestern  United  States  events. 
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where  C*  is  defined  as  travel  time  divided  by  average  Q  along  the  travel  path, 
and  I(w)  is  the  displacement  response  of  the  LRSM  seismographs.  The  chosen 
values  of  t*  varied  according  to  source  region  (Der,  1976);  they  were  0.2  for 
events  east  of  the  Wasatch  Front  (see  Figure  9)  and  0.4  for  events  to  the  west 
of  it,  except  for  events  G  and  F,  which  were  arbitrarily  assigned  an  inter¬ 
mediate  value  of  0.3.  Spectra  of  noise  samples  preceding  the  first  arrival 
of  the  P  wave  were  also  computed  and,  to  retain  the  comparative  S/N  ratio, 
these  spectra  were  also  adjusted  by  the  factor  exp  [nf t*] /l(u>) . 

The  corner  frequencies  and  low-frequency  asymptotic  levels  of  these  spec¬ 
tra  were  determined  in  a  visual,  but  objective,  manner.  To  conform  to  theo- 

rectical  predictions  of  von  Seggern  and  Blandford  (1972),  high-frequency  asymp 
-2 

totic  slopes  of  f  were  uniformly  drawn  for  the  explosion  spectra.  This 

-3 

slope  is  a  good  fit  in  almost  all  cases,  but  f  would  be  acceptable  for  some 

spectra,  such  as  for  GREELEY  at  RK-ON  or  BENHAM  at  NP-NT.  For  the  earthquakes 
—2  —3 

either  an  f  or  f  slope  was  used  to  fit  the  high-frequency  falloff;  either 
slope  can  be  predicted  from  various  earthquake  source  theories.  Diagrams  of 
the  long-period  level  versus  corner  frequency  are  shown  in  Figure  19,  where 
a  fairly  clear  separation  of  earthquakes  from  explosions  occurs.  Event  1 
(DURYEA)  fails  to  separate  for  RK-ON  even  though  the  parameters  seem  to  be 
well-determined  from  its  spectrum  in  Figure  18.  Event  7  (GASBUGGY)  fails  to 
separate  at  NP-NT,  but  its  spectrum  in  Figure  17  is  close  to  the  noise  level 
and  possibly  yields  an  inaccurate  spectral  shape.  Event  L  (Coyote  Mtn.  EQ) 
lies  in  the  earthquake  population  for  both  stations.  (Previously  it  was 
determined  that  it  had  rather  low  S-wave  and  long-period  P-wave  excitation, 
and  in  the  next  section  we  show  that  its  Mg  is  low.)  In  the  final  analysis, 
the  data  in  Figure  19  could  be  interpreted  to  reflect  one  physical  basis  for 
discrimination:  that  for  a  given  low-frequency  level  indicative  of  source 
size,  the  characteristic  rise  time  for  the  earthquake  displacement-time  func¬ 
tion  is  longer  than  that  for  an  explosion. 


Der,  Z.  (1976).  On  the  existence,  magnitude,  and  causes  of  broad  regional 
variations  in  body-wave  amplitudes.  SDAC  Report  No.  TR-76-8,  Teledyne 
Geotech,  Alexandria,  Virginia. 
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Figure  19  Long-period  spectral  level  uQ  versus  corner  frequency  for 
waves  recorded  at  NP-NT  and  RK-ON. 


Figure  19  (cont.)  Long-period  spectral  level  uQ  versus  corner  frequency  for 

P  waves  recorded  at  NP-NT  and  RK-ON. 
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Figure  20  illustrates  the  utility  of  the  Mg  -  m^  plot  for  distinguishing 
source  type.  However,  in  our  group  three  earthquakes  at  lower  magnitude  (E, 

F,  and  J)  are  definitely  near  the  explosion  population  and  a  fourth,  event  L, 
is  somewhat  anomalous  at  a  higher  magnitude.  Because  of  the  importance  of  the 
-  m^  discriminant,  possible  causes  of  such  anomalous  behavior  will  be  exam¬ 
ined  in  detail. 

Thatcher  and  Hamilton  (1973)  have  previously  noted  the  unusual  character¬ 
istics  of  Event  L,  the  Coyote  Mountain  earthquake.  They  reported  a  small 
source  dimension  (^3  km)  and  a  relatively  large  stress  drop  (^80  bars).  Source 
dimensions  of  the  largest  explosions  in  the  study  sample  have  been  variously 
estimated  to  be  on  the  order  of  1  km,  so  that  the  small  difference  in  source 
dimension  is  compatible  with  the  proximity  of  the  Coyote  Mountain  earthquake 
to  the  explosion  population  in  Figure  20.  The  relatively  small  source  dimen¬ 
sion  for  this  earthquake  causes  the  corner  frequency  in  Figure  19  to  differ 
little  from  the  largest  explosions  in  our  sample. 

Thatcher  and  Hamilton  calculated  the  seismic  moment  of  the  Coyote  Moun- 

24 

tain  earthquake  from  Rayleigh  waves  to  be  roughly  7x10  dyne* cm,  while  Aki 

9  / 

and  Tsai  (1972)  calculated  that  of  event  9,  BOXCAR,  to  be  roughly  2x10 
Thus,  if  moment  rather  than  Ms  were  plotted  in  Figure  20,  the  Coyote  Mountain 
earthquake  would  separate  clearly  from  the  explosions.  The  observed  Mg  for 
this  earthquake  is  low  due  probably  to  the  focal  depth's  location  near  the  node 
of  the  Rayleigh-wave  displacement-depth  function  in  combination  with  predomi¬ 
nantly  strike-slip  motion.  To  illustrate  this  condition,  Figure  21  shows  the 
relative  Rayleigh-wave  excitation  of  various  infinitesimal  sources  in  a  lay¬ 
ered  halfspace  chosen  to  represent  the  Basin-Range  structure.  The  zeroes  for 
spectral  excitation  at  .05  Hz  will  change  only  slightly  with  various  continen¬ 
tal  structures.  Note  that  if  the  Coyote  Mountain  earthquake  is  a  strike-slip 
fault  at  depths  of  10-13  km  (Thatcher  and  Hamilton,  1973),  then  the  20-second 
LR  excitation  would  be  roughly  five  times  less  than  an  explosion  of  equal 
moment  near  the  surface.  In  fact,  however,  the  moment  ratio  is  roughly  three 
for  the  Coyote  Mountain  earthquake  versus  BOXCAR,  and  thus  the  20-second  exci¬ 
tation  should  be  nearly  equal,  as  the  Mg  values  in  Figure  20  indicated. 
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Figure  20  versus  ra^  for  selected  Southwestern  United  States  events. 
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Figure  21  Relative  Rayleigh-wave  excitation  at  a  period  of  20  seconds  for 
various  source  types  in  a  Basin-Range  structure. 


Event  J,  the  BENHAM  aftershock,  is  the  most  important  earthquake  in  our 
study  because  of  its  proximity  to  the  underground  explosions  at  NTS.  This 
fact  enables  us  to  ignore  path  effects  in  discrimination  and  to  focus  on  the 
source  contribution  itself.  Points  discussed  earlier  in  the  review  of  dis¬ 
crimination,  and  the  relation  of  LR  and  P  excitation  to  focal  depth  for  a 
Basin-Range  structure  (shown  in  Figure  21  and  22),  suggest  that  for  small 
sources  discrimination  of  shallow  earthquakes  from  explosions  should  be 
rather  difficult  on  the  basis  of  Mg  -  m^  alone.  Basham  et  al.  (1970),  Peppin 
and  McEvilly  (1974),  and  Savino  et  al.  (1971)  have  all  noted  that  the  BENHAM 
aftershock,  reported  to  have  occurred  at  3.5  km  depth  (Hamilton  et  Healy, 

1969),  had  low  Mg.  The  focal  mechanism  in  Figure  12  is  dominated  by  dip-slip 
motion  on  a  high-angle  fault.  Figure  21  shows  that  excitation  tends  to  zero 
in  dip-slip  sources  near  the  surface  as  the  fault  plane  steepens.  Thus,  the 
cause  of  low  Mg  for  this  event  has  probably  been  identified  but  m^  must  still 
be  considered.  According  to  Figure  22,  the  of  the  BENHAM  aftershock  should 
be  low  relative  to  an  explosion  of  equal  moment  in  the  upper  layers  at  NTS. 

If  a  correction  for  source  depth  were  made  to  observed  m^,  then  the  BENHAM 
aftershock  would  move  farther  right  into  the  explosion  population  in  Figure  20. 
Part  of  the  physical  basis  which  holds  this  earthquake  among  the  explosions  is 
a  high  corner  frequency  and  small  source  dimension,  a  fact  represented  in  the 
spectral  parameters  illustrated  in  Figure  19. 

Event  E,  the  Denver  earthquake,  is  interesting  because  it  was  probably 
induced  by  hydraulic  lubrication  of  pre-existing  fractures.  It  took  place  in 
an  environment  considerably  altered  from  the  natural  stress  state  by  forced 
pumping  of  liquid  wastes  at  the  Rocky  Mountain  Arsenal.  This  earthquake  was 
at  a  relatively  shallow  depth  of  5  km,  and  Major  and  Simon  (1968)  determined 
that  it  was  predominantly  strike-slip  on  a  nearly  vertical  fault.  Figure  21 
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Basham,  P.  W.,  D.  H.  Weichert,  and  F.  M.  Anglin  (1970).  An  analysis  of  the 
BENHAM  aftershock  sequence  using  Canadian  recordings.  J.  Geophys,  Res., 
75,  1545-1556. 

Major,  M.  W. ,  and  R.  B.  Simon  (1968).  A  seismic  study  of  the  Denver  (Derby) 
Earthquakes.  Quarterly  of  the  Colorado  School  of  Mines,  63,  9-56. 


-159- 


RELATIVE  EXCITATION 


P  WAVE  SYNTHESIS  FROM 

II  STEP  SOURCE  TIME  FUNCTION 

2)  EARTH  MODEL  YFI11  (GLOVER  &  ALEXANDER  1970) 

3)  NO  ATTENUATION 

4)  BENIOFF  SHORT  PERIOD  RESPONSE 
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Figure  22  Relative  amplitude  of  teleseismic  P  waves  from  dilatational 
sources  of  identical  moment  at  various  depths  in  a  Basin- 
Range  structure. 
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indicated  that  such  a  mechanism  at  such  depth  should  not  be  low  on  the  M 

s 

scale.  Therefore,  an  anomalously  high  must  be  the  reason  why  this  event 
falls  close  to  the  explosions.  Wyss  and  Molnar  (1972)  inferred  a  fault  dimen¬ 
sion  of  ^3  kin,  which  is  above  average  for  an  earthquake  at  this  magnitude. 

Indeed,  Major  and  Simon  concluded  that  the  fault  dimension  may  be  as  large  as 
10  km.  The  P-wave  corner  frequency  data  in  Figure  19  supports  these  interpre¬ 
tations  that  a  source  dimension  much  larger  than  that  of  explosions  at  this 
magnitude  exists.  Thus,  the  anomalous  Mg  -  m^  of  this  event  is  not  compatible 
with  known  source  characteristics. 

Event  F,  the  Baja  California  earthquake,  is  one  of  many  that  Thatcher 
(1972)  studied  in  the  surrounding  area.  The  focal  depth,  taken  as  33  km  from 
the  NEIS  list,  is  not  accurate  and  it  could  actually  be  anywhere  within  the 
crust.  In  the  next  section,  which  deals  with  complexity,  event  F  is  shown  to 
have  a  high  P  coda,  suggesting  a  focus  well  below  the  surface.  Also,  the 
mechanism  of  this  event  could  not  be  specified  and  it  was  arbitrarily  assumed 
to  be  strike-slip  (Table  V).  Thatcher  observed  that  most  of  the  Baja  events 
in  his  data  set,  including  this  one,  had  low  long-period  spectral  level  for 
their  local  magnitude  M^,  measured  at  high  frequencies.  He  concluded  that 
source  dimensions  for  these  events  are  relatively  small,  nearly  like  explo¬ 
sions.  On  this  basis  the  relatively  low  Mg  of  this  event  could  be  explained. 
Furthermore,  if  the  depth  were  such  (=10km)  that  the  LR  excitation  was  con¬ 
siderably  diminished  for  a  strike-slip  source  as  shown  in  Figure  21,  then  low 
Mg  would  certainly  be  expected.  Note  also  that  Figure  21  implies  that  very 
shallow  dip-slip  sources  would  also  have  low  Mg. 

Rayleigh  waves  from  events  E,  F,  and  J  were  well  recorded  on  the  network 
used  for  this  study,  with  at  least  five  observations  each.  Therefore,  low 
Mg  is  probably  not  due  to  the  stations'  all  being  near  a  node  of  the  radiation 
pattern.  Examination  of  LQ/LR  radiation  patterns  reveals  that  for  cases  where 
LR  excitation  is  small,  LQ  excitation  is  often  large,  a  fact  making  the  sum 
of  Love-wave  and  Rayleigh-wave  Ms's  an  improved  discriminant.  Figure  23  shows 
the  results  from  this  combination  for  all  events.  Although  now  a  better  separa¬ 
tion  of  the  two  populations  exists,  it  is  not  sufficient  to  insure  with  high  con¬ 
fidence  identification  of  Southwestern  United  States  explosions  and  earthquakes. 


-161- 


Complexity 

The  complexity  of  the  P  signal  is  considered  a  source  discriminant  because 
it  stems  either  from  depth  of  focus,  which  would  entail  additional  phases  from 
the  upgoing  wave  (the  free-surface  reflections  will  be  most  pronounced)  or  to 
prolonged  signal  generation  from  complex,  spatially  extended  fault  ruptures. 

The  complexity  measure  has  been  computed  from  the  vertical  short-period  tape 
recordings  available  for  the  event  set.  The  complexity  computation  method 
was  to  window  the  recordings  from  30  sec  before  the  P  arrival  to  35  sec  after, 
square  the  data  samples  in  these  windowed  traces,  integrate  with  a  3-sec  mov¬ 
ing  average  (moved  1/20  sec  for  each  output  point),  and  take  the  square  root 
of  the  integrated  trace.  The  complexity  value  is  formed  from  this  final  trace 
and  results  from  subtracting  a  noise  estimate  based  on  the  10  sec  preceding  P 
from  the  trace  comprising  35  sec  after  it  and  then  ratioing  the  5-35  sec  por¬ 
tion  after  P  to  the  0-5  sec  portion  after  P.  Event  mean  complexities  are  shown 
in  Figure  24  versus  P-wave  m^.  The  range  in  complexity  values  is  rather  narrow, 
except  for  three  events,  which  are  the  Baja  California,  Coyote  Mountain,  and 
Caliente  earthquakes.  Each  of  these  three  earthquakes  clearly  shows  its  complex 
character  in  the  NP-NT  and  RK-ON  recordings  of  Figures  17  and  18  respectively. 
Source  depth  alone  cannot  fully  explain  the  large  complexity,  and  a  preferable 
interpretation  would  be  a  prolonged  rupture  duration  or  extended  source  time 
function. 

Figure  24  shows  that  there  is  some  separation  between  earthquakes  and 
explosions  but  whether  or  not  it  is  due  to  the  generally  higher  explosion  m^'s 
is  not  clear.  The  data  in  Figure  24  must  also  be  questioned  because  the  data 
points  are  event  means  over  varying  sets  of  stations.  Lambert  et  al.  (1969) 
and  Davies  and  Julian  (1972)  showed  that  complexity  is  highly  dependent  on  the 
receiver  location  and  should  exhibit  scatter  among  stations  equivalent  to  any 
proposed  source  differences.  Because  of  data  recovery  problems  and  varying 
station  operation  times,  the  network  is  different  for  each  event  of  Figure  24, 
and  significant  biases  are  possible.  It  may  be  preferable  to  consider  com¬ 
plexity  on  a  single-station  basis  as  shown  for  NP-NT  and  RK-ON  in  Figure  25. 
Because  of  data  unavailability  or  poor  quality,  several  events  are  missing 
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Figure  25  (cont.)  mj,  versus  complexity  of  signals  recorded  at  NP-NT  and 

RK-ON  from  selected  Southwestern  United  States  events. 
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from  this  plot,  but  there  are  enough  to  indicate  that  complexity  fails  to  separ¬ 
ate  the  two  populations.  Note,  however,  that  no  explosion  has  a  high  complexity 
and  that  the  range  of  explosion  complexities  is  much  less  than  that  of  earth¬ 
quakes.  The  large  complexity  values  at  RK-ON  for  RULISON  and  GASBUGGY  are  not 
reliable  due  to  the  regional  recording  distance  and  the  low  S/N  ratios  (see 
Figure  18).  Generally,  the  larger  explosions  have  higher  complexity,  which 
may  be  due  to  delayed  tectonic  strain  release. 

P-Wave  Spectral  Ratios 

A  short-period  spectral  ratio  has  been  computed  according  to 
SPSR  -  ,1'87  ,0-78 

{.56  A<f>«£'4.47  A<£>d£ 

where  sums  over  discrete  Fourier  amplitude  coefficients  replace  the  integrals. 
The  spectrum  A(f)  is  computed  on  the  first  6.4  seconds  of  the  P  wave.  The 
network  average  spectral  ratios  are  plotted  versus  m^  in  Figure  26.  At  a 
given  m,  the  explosion  spectral  ratios  are  generally  higher  than  earthquake 

“  _2  _3 

ones.  This  observation  agrees  with  the  id  and  to  high-frequency  slope  mod¬ 
els  for  explosions  and  earthquakes,  respectively,  outlined  theoretically  in 
an  earlier  section  of  this  report  and  supported  by  real  spectra  shown  in  Fig¬ 
ure  17  and  18.  The  observed  decrease  of  spectral  ratio  as  a  function  of  yield 
is  also  in  agreement  with  the  theory.  Still,  no  complete  separation  is  made  on 
the  basis  of  spectral  ratio  for  our  events.  Event  H  in  Figure  26,  the  Borrego 
Mountain  earthquake,  falls  farthest  into  the  explosion  population,  but  because 
of  being  overdriven  on  LRSM  recordings,  its  spectral  ratio  is  from  one  station 
only  (WH2YK)  and  it  could  be  heavily  biased  by  a  station  effect.  Attenuation 
corrections  were  made  to  A(f)  (see  P-wave  spectrum  section  above)  with  the 
aim  of  removing  variable  path  effects,  thus  determining  a  quantity  more  repre¬ 
sentative  of  the  source.  The  resulting  spectral  ratios  (not  shown)  separated 
no  more  clearly  than  those  shown  in  Figure  26. 

Following  Anglin  (1971),  P-wave  spectral  ratio  vs.  complexity  was  plotted 
in  Figure  27.  Anglin  obtained  a  fair  separation  of  his  populations,  composed 
of  Eurasian  events,  mostly  on  the  basis  of  spectral  ratio  ("third  moment  of 
frequency"  in  his  paper)  with  little  classification  power  stemming  from 
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Figure  26  ra^  versus  P-wave  spectral  ratio  for  selected  Southwestern 
United  States  events. 
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complexity.  No  attenuation  corrections  were  made  in  his  work,  and  the  "third 
moment  of  frequency"  may  be  measuring  a  difference  in  source-path  effects 
either  in  addition  to  or  rather  than  an  actual  difference  in  source  spectra 
between  explosions  and  earthquakes.  For  data  shown  in  Figure  27,  there  appears 
to  be  more  discrimination  power  in  complexity  than  in  the  spectral  ratio, 
which  is  corrected  for  attenuation,  but  there  is  no  satisfactory  separation 
of  events  in  this  plot. 

Rayleigh-Wave  Spectral  Ratios 

Von  Seggern  and  Lambert  (1970)  studied  the  discrimination  power  of  spec¬ 
tral  ratio  from  long-period  Rayleigh  waves  for  a  sample  of  global  events. 

Their  results  show  that  while  explosions  generate  LR  waves  with  a  fairly  con¬ 
stant  observed  spectral  ratio,  earthquakes  can  exhibit  values  in  the  same 
range,  but  also  much  higher  or  much  lower.  The  LR  spectral  ratio  is  computed 
for  our  events  according  to 

=  .071  .045 

(046  A^df/,{o28  A(f)df 

where  sums  of  discrete  transform  points  replace  the  integrals.  The  time  win¬ 
dow  was  dependent  on  path  length,  from  the  expected  arrival  time  of  T  =  50+ 
sec  to  that  of  T  =  10  sec.  This  definition  is  close  to  von  Seggern  and 
Lambert's,  and  the  resulting  network  average  values  are  plotted  versus  Mg  in 
Figure  28.  This  plot  shows  no  separation  between  earthquakes  and  explosions 
and,  in  fact,  shows  the  same  general  picture  as  in  von  Seggern  and  Lambert, 
with  explosions  confined  to  median  values  and  earthquakes  ranging  on  both 
sides.  These  results  agree  with  theoretical  predictions  of  Rayleigh-wave 
excitation  from  double-couple  sources  which  generate  a  variety  of  Rayleigh-wave 
spectral  shapes,  dependent  on  source  depth,  fault  orientation,  and  the  azimuth 
of  the  receiver. 
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MULTIPLE  DISCRIMINATION 


In  the  previous  section,  the  capability  of  various  seismic  measurements 
to  classify  the  explosions  and  earthquakes  in  our  event  set  was  examined. 

These  discriminants  were  largely  discussed  independently  and  more  qualitatively 
than  quantitatively.  In  this  section,  by  applying  statistical  techniques, 
the  power  of  the  various  discrimination  measures  will  be  examined  quantitatively 
in  a  unified  approach.  A  multiple  discrimination  technique  is  used  to  identify, 
rank,  and  fully  utilize  meaningful  discrimination  parameters.  A  list  of  the 
parameters  used  in  this  experiment  is  given  in  Table  VI;  all  of  these  parameters 
were  displayed  in  one  or  more  figures  in  the  previous  section,  and  their  values 
for  the  individual  stations  can  be  found  in  Table  IV. 

Multiple  discrimination  is  used  because  significant  classification  infor¬ 
mation  is  contained  in  almost  all  parts  of  the  seismogram  and  in  both  the  fre¬ 
quency  and  time  domains.  Although  the  main  effort  is  directed  toward  network 
discrimination,  the  discriminating  capability  of  single  sites — as  opposed  to 
the  entire  network — will  be  examined.  In  this  experiment  our  data  set  is 
treated  as  a  training  set,  where  correct  classification  is  known  a  priori. 

Booker  and  Mitronovas  (1964)  did  some  pioneering  work  on  seismic  multiple 
discrimination  based  upon  finding  discriminant  functions  for  known  groups,  but 
investigators  have  performed  little  subsequent  work  using  more  than  a  few  para¬ 
meters. 

Description  of  Discrimination  Experiments  to  be  Performed 

The  data  base  shown  in  Table  IV  was  used  in  three  separate  experiments 
aimed  at  evaluating  the  effectiveness  of  several  possible  discriminants.  In 
the  first  experiment,  data  from  all  network  stations  were  combined  using  Ring- 
dal's  (1976)  method  to  give  estimates  of  the  event  magnitudes.  These  magni¬ 
tude  estimates  were  used  to  calculate  linear  discriminant  functions  for  se¬ 
lected  subsets  of  the  ten  variables  defined  as  discriminants  here.  In  the  sec¬ 
ond  experiment,  discriminant  functions  were  calculated  as  before,  but  using 
data  from  one  station  in  place  of  the  estimates  determined  by  the  whole  net¬ 
work.  In  the  third  experiment  single-station  discriminant  functions  were 


Ringdal,  F. ,  (1976).  Maximum-likelihood  estimation  of  seismic  magnitude. 
Bull.  Seism.  Soc.  Am.,  66,  789-802. 
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TABLE  VI 


Network  Discrimination  Parameters  for  the  Multiple  Discrimination  Exper iment 


Event 

SPP 

SPS 

LPP 

LPS 

Ul 

LR 

log  SPSR 

log  LPSR 

log  comp 

1st  mtn 

1 

4.97 

*4.97 

*4.96 

*4.10 

3.78 

4.03 

-1.453 

0.412 

0.652 

0. 

2 

5.53 

4.23 

3.55 

4.06 

4.62 

4.17 

-1.322 

0.221 

0.635 

0. 

3 

6.22 

*5.18 

3.68 

4.89 

5.23 

4.92 

-1.838 

0.  308 

0.604 

0. 

4 

6.37 

5. 11 

4.07 

5.17 

5.86 

5.36 

-2.187 

0.290 

0.705 

0. 

5 

5.77 

4.26 

3.39 

4.42 

4.71 

4.68 

-1.863 

0.221 

0.596 

0. 

6 

5.53 

*4.73 

3.61 

4.76 

5.13 

4.67 

-1.444 

0.248 

0.588 

0. 

7 

4.64 

*5.07 

*4.28 

3.35 

3.71 

3.77 

-1.237 

0.180 

0.730 

0. 

8 

6.34 

5.14 

4.37 

5.08 

5.25 

5.34 

-1.995 

0.240 

0.822 

0. 

9 

6.37 

5.20 

4.26 

5.40 

5.27 

5.58 

-1.884 

0.103 

0.799 

0. 

10 

6.43 

4.79 

4.72 

5.61 

6.11 

5.63 

-2.023 

0.189 

0.858 

0. 

11 

4.74 

4.68 

*4.18 

4.67 

3.79 

3.87 

-1.183 

0.212 

0.830 

0. 

A 

5.04 

*4.51 

4.13 

4.75 

4.86 

4.82 

-1.628 

0.172 

0.720 

1. 

B 

5.37 

*4.56 

4.26 

5.32 

5.72 

5.35 

-1.697 

0.043 

0.821 

1. 

C 

5.51 

5.05 

5.02 

5.85 

5.46 

5.65 

-1.771 

0.053 

1.094 

1. 

D 

4.45 

4.20 

4.03 

5.09 

5.47 

5.01 

-1.125 

0.131 

0.578 

1. 

E 

4.89 

*4.89 

3.94 

4.44 

4.55 

4.19 

-1.447 

0.088 

1.008 

1. 

V 

4.75 

*4.67 

*4.31 

4.04 

4.45 

4.14 

+-1.344 

0.248 

1.362 

0. 

C 

4.84 

4.28 

4.32 

5.11 

5.09 

4.78 

-1.711 

0.517 

0.945 

1. 

11 

6.33 

6.14 

5.63 

6.08 

6.54 

6.54 

-1.095 

+-0.118 

1.065 

1. 

I 

4.32 

*4.32 

3.58 

4.39 

4.48 

4.08 

-1.460 

-0.012 

0.866 

1. 

J 

4.44 

*4.65 

*4.29 

4.05 

4.20 

3.61 

-1.686 

0.055 

0.981 

0. 

K 

4.39 

4.18 

*4.41 

4.08 

4.82 

4.00 

-1.396 

0.445 

0.815 

1. 

L 

5.63 

*4.36 

4.52 

5.34 

6.11 

5.39 

-1.836 

tO. 361 

1.565 

1. 

*  no 

detection;  magnitude 

is  average  of 

noise 

levels 

+  spectral  r 

atio  at 

all  detecting 

stations  was 

affected 

by  noise 

SPP : 

Short- 

period 

P\ 

SPS: 

Short- 

period 

Smb 

LPP : 

Long-p 

eriod  P 

mb 

LPS : 

Long-p 

eriod  S 

"b 

LO:  Love-wave  M 

s 


LR:  Rayleigh-wave  M^ 

log  SPSR:  log  short-period  P-v;ave  spectral  ratio 

log  LPSR:  log  long-period  LR-wave  spectral  ratio 

log  comp:  log  complexity 

1st  mtn:  polarity  of  first  motion 
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averaged  to  perform  discrimination  based  upon  the  network.  Now,  before  dis¬ 
cussing  the  experimental  results,  Ringdal's  method  is  presented. 

Estimation  of  Event  Magnitudes 

The  magnitude  m  of  a  seismic  event  is  conventionally  estimated  by  com¬ 
puting  the  arithmetic  mean  nK  of  the  magnitudes  that  each  detecting  station 
measured.  During  this  process  certain  information  is  lost,  namely  that  the 
signal  was  less  than  the  noise  level  at  non-detecting  stations.  Ignoring 
those  stations  where  the  signal  was  below  the  threshold  of  detection  elimi¬ 
nates  what  would  generally  be,  for  a  perfectly  detecting  (noise-free)  network, 
the  smaller  values  of  iik.  The  arithmetic  mean  would  thus  tend  to  be  biased 
so  that  nn  >  m.  The  maximum-likelihood  estimator  that  Ringdai  (1976)  proposed 
is  a  method  of  incorporating  into  the  computation  the  constraint  that,  for 
non-detecting  stations,  the  station  magnitude  is  less  than  some  specified 
noise  level. 


For  an  event  of  magnitude  m,  assume  that  the  set  of  magnitudes  nn,  which 


would  be  measured  by  seismic  stations  of  a  noise-free  network,  comprises  a 
Gaussian  distribution  having  mean  p  and  standard  deviation  o.  When  the  num¬ 
ber  of  stations  in  the  network  becomes  infinite,  the  estimated  p  =  m.  If  all 
stations  detect  the  signal,  then  the  best  estimate  of  p  is  the  arithmetical 


mean  m,.  If,  however,  the  instantaneous  noise  level  a.  at  some  stations  is 
1  1 


too  high  for  the  signal  to  be  detected  (i.e.,  cu  >  itu,  assuming  S/N  =  1  is 


sufficient  for  detection),  then  Ringdal's  method  permits  calculating  the  most 
likely  Gaussian  distribution  such  that  the  detecting  stations  measure  magni¬ 
tudes  m.  and  such  that  m.  <  a.  at  the  non-detecting  stations.  The  mean  value 
i  l  i 

of  this  most  likely  distribution  is  the  maximum-likelihood  estimate  of  m, 

For  the  purpose  of  calculating  the  instantaneous  noise  levels 

have  been  measured  on  the  LRSM  recordings  at  the  expected  arrival  times  of 
short-period  P  and  S  and  long-period  P,  S,  LQ,  and  LR  when  these  signals 
could  not  be  detected  on  the  seismograms. 


Consider  a  family  of  Gaussian  distributions  where  each  is  characterized 


by  the  parameters  (hj »  Oj).  For  any  given  member  of  this  family  the  probabil¬ 


ity  P^(m^)dm^  that  a  detecting  station  would  measure  the  magnitude  to  be  in 

2* 


the  range  dnK  about  m^  is 


P1(m1;M j ,Oj)dmi 


/27o.eXP 

J 


/ID. 

4  1 


h 


a . 
J 


dm . 


(1) 
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For  a  station  with  a  sharp  detection  threshold  a^,  the  probability  P2(nu  <  a.) 
that  the  signal  would  be  less  than  a ^  and  therefore  undetected  is 


p2<mi  <  a±’>  °j) 


/“£  exp 


dm . 


(2) 


Ringdal's  maximum-likelihood  estimate  of  m  is  thus  the  value  of 
the  function 


for  which 


(3) 


n  P  (m  ;u  ,a  )  dm  •  n  P„(m.  <  a  ;  y.,  a.) 

J  J  1  Z  X  1  J  J 

detecting  stations  non-detecting  stations 

attains  the  greatest  value.  One  may  calculate  simply  by  evaluating  the 

likelihood  function  L  for  all  values  of  y.  and  o.  incremented  over  a  suitable 

J  1 

range  and  choosing  the  pair  (Uj,cO  for  which  L  attains  a  maximum. 


In  several  instances  (short-period  S  for  events  1,  3,  6,  7,  E,  F,  I,  J, 
and  L;  long-period  P  for  events  1,  7,  11,  F,  J,  and  K;  long-period  S  for  event 
1)  neither  m^  nor  y^^  could  be  calculated  since  all  stations  in  the  network 
measured  only  noise.  For  these  cases  the  arithmetic  mean  of  the  noise  level 
was  computed;  this  value  was  used  as  an  approximation  to  m,  although  in 
most  cases  this  mean  is  expected  to  be  higher  than  the  true  m.  For  many  events 
LQ  could  not  be  measured  at  most  stations  because  the  long-period  horizontal 
instruments  were  overdriven.  In  four  instances  (events  10,  B,  C,  and  H) ,  sup¬ 
plemental  information  was  used  to  infer  the  magnitude  of  LQ  from  the  measured 
value  of  LR  (Table  IV).  Arithmetic  means  were  calculated  for  the  logarithms 
of  the  short-period  and  long-period  spectral  ratios  and  for  the  logarithms  of 
the  complexities.  For  calculating  the  network  means  of  these  last  three  vari¬ 
ables,  all  values  judged  unreliable  because  of  interfering  noise  were  deleted 
because  in  these  cases  the  Ringdal  estimate  y  ^  Is  meaningless.  Finally,  each 
event  was  assigned  a  tenth  variable,  whose  value  was  unity  if  any  network  sta¬ 
tion  detected  a  clear  dilational  first  motion  and  which  was  otherwise  zero. 


Discriminant  Fuctions 

The  data  for  the  first  discrimination  experiment  consist  of  the  twenty- 
three  vectors  7=  (X^,  X2,...,  Xp)  (listed  in  Table  VI)  where  Xi  is  one  of 


the  ten  variables  whose  network  mean  values  have  been  calculated  using  proce¬ 
dures  described  above.  Because  it  is  known  ji  priori  which  of  these  data  vec¬ 
tors  represent  explosions  and  which  represent  earthquakes,  this  set  of  vectors 
can  be  used  as  a  "training  set"  in  order  to  provide  discriminant  functions  for 
classifying  future  data  vectors.  However,  because  these  discriminant  functions 
could  not  be  correctly  applied  to  events  outside  the  Western  United  States, 
this  study  will  not  focus  on  the  numerical  coefficients  of  these  functions  but 
rather  on  the  subset  of  variables  X^,  X^  which  gives  the  "best"  dis¬ 

criminant.  Thus,  for  each  discriminant  tested,  the  a  posteriori  probabilities 
of  correctly  identifing  each  event  of  the  training  set  are  examined  by  apply¬ 
ing  the  calculated  discriminant  function  to  it.  In  this  manner  one  finds 
which  subsets  of  variables  yield  the  best  separation  between  earthquakes  and 
explosions  in  the  training  set. 

This  report  assumes  that  the  eleven  explosion  and  twelve  earthquake  data 
vectors  are  members  of  two  p-dimensional  multinormal  populations  with  mean 
vectors  and  and  with  a  common  covariance  matrix  Z.  Assumption  of  equal 
covariance  for  the  two  populations  of  seismic  data  leads  only  to  small  errors 
(Shumway  and  Blandford,  1970).  The  probaoility  of  classifying  a  data  vector 
X  as  an  explosion,  P^(X),  and  the  probability  of  classifying  it  as  an  earth¬ 
quake,  P2(X),  are  assumed  a  priori  to  be  0.5,  although  for  this  data  set  they 
are  actually  11/23  and  12/23,  respectively.  Under  these  assumptions  the  clas¬ 
sification  probabilities  are  given  by  (Anderson,  1958) 

P.CT)  =  - 577 - 777  exp  [-%OT-~T)  Z  1  OT-~T, )  ]  j  =  1,2  (4) 

3  (2tt)H/Z|z|  3  ~  3 

Here  the  "-1"  exponent  implies  the  inverse  matrix,  and  the  prime  implies  the 

transpose  matrix. 

Taking  the  ratio, 

Pl(“T)  _  _ 

— »'  =  exp  (X  •  D  +  C},  (5) 

P2(  X  ) 

Shumway,  R.  H. ,  and  R.  R.  Blandford  (1970).  A  simulation  of  seismic  discrimin¬ 
ant  analysis.  SDL  Report  No.  261,  Teledyne  Geotech,  Alexandria,  Virginia. 

Anderson,  T.  W.  (1958).  An  Introduction  to  Multivariate  Statistical  Analysis. 

New  York,  NY.,  John  Wiley  and  Sons, 
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where 


(6) 


D  =  £_1  (  M1  -  W2) 

and 

c  =  -h  (  +  v2)'  £_1  (  (7) 

The  criterion  upon  which  discrimination  is  based  is  thus 

P1(-jT )  >  P2(T)iffT-  D  +  C  >  0.  (8) 

We  shall  refer  to  X  *  D  +  C,  sometimes  known  as  the  group  separation  function, 
simply  as  the  discriminant  function.  We  may  rewrite  the  equation  (5)  as 

P.OT)  __  __  __  __ 

— —  — -  =  exp  {X  •  D  +  C.  -  X  .  D„  -  C„  },  (9) 

P2(  X  ) 

where 


=  f1  y  j  «  1,2  (10) 

and 

CJ  J-1.2  W1> 

We  shall  refer  to  X  •  and  X  *  D2  +  as  group  classification  functions, 

although  they  also  are  sometimes  called  discriminant  functions,  a  name  we  re¬ 
serve  for  X  '•  D  +  C.  Adding  unity  to  both  sides  of  equation  (9),  we  have 

Pj^dT)  p2(lT)  exp (  ~x~ •  Dx  +  expdT*  D2  +  C2) 

p2nr>  p2nr)  expor-"^  +  c2>  expor*-^  +  c2) 

Taking  reciprocals,  we  find  for  the  a  posteriori  probability  of  classifying  X 
as  an  earthquake 

-  exp  (~x"  *  "^2  +  C2^ 

p2(  x  )  -  - zr“ - - 

exp (  x  '  +  c^)  +  exP(  X  *  &2  +  C2)  (12) 


or  as  an  explosion 


__  exp  (X  ♦  +  C^) 


PX(X) 


exp  (X  •  +  C^)  +  exp  (X 


d2  +  c2) 


(13) 


In  actual  practice  the  population  parameters  y2»  and  I,  are  unknown  and 
must  be  estimated  by  the  parameters  ,  X2  and  S,  which  are  calculated  from 
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the  sample  of  twenty-three  events.  For  every  set  of  variables  X  selected  as 
a  possible  discriminant,  we  thus  find  the  linear  discriminant  function  and 
group  classification  functions  (evaluated  for  arbitrary  data  vector  X)  to  be 


and 


D  •  X  +  C  =  X' 


D.  •  X  +  C. 
J  J 


S  1  (  xx  -  x2)  -hi  \  +  x2r 

=  Tr  s_1  IT.  -h  Y'.  s-1  X 

j  j  -  j 


j 


1.2 


(14) 


(15) 


respectively.  A  measure  of  the  effectiveness  of  the  discriminant  function  is 
then  given  by  the  a  posteriori  probability  for  the  k-th  event 

~  x  (16) 

exp(  xk  •  D  +  C  ) 

p.(  xT)  - 5 J - — — XI -  j  =  k  =  1>23 

3  exp(  xk  *  D2  +  C2)  +  exp (  Xfc  •  D2  +  C2) 

Another  measure  of  the  separation  of  the  two  groups  is  provided  by  Wilks' 
A-criterion,  defined  as 


A  (  X  ) 


det  W 
det  T 


(17) 


where  W  is  the  within  -group  cross-product  matrix  and  T  is  the  total  cross- 
product  matrix.  Letting  denote  the  i-th  component  of  the  k-th  vector  of 

the  j-th  group,  these  matrices  are  given  by 


and 


W 


il 


2 

=  Z 

j=l 


n. 

3 

Z 

k=l 


\ 


-  Xj't) 


T 

i£ 


2  nj 

A  A  «na  -  S..1>V  “••I1 

3=1  k=l  J 


(18) 


(19) 


where  the  dot  subscript  notation  indicates  averaging  over  that  subscript  which 
has  been  replaced  by  a  dot.  Since  Wilks'  A-statistic  is  poorly  tabulated,  we 
shall  instead  measure  group  separation  (actually  the  hypothesis  tested  is  "hat 
of  equality  of  group  means)  by  means  of  a  statistic  which  is  derived  from  A 
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and  which  has  approximately  an  F-distribution  (Jennrich,  1977).  The  A-cri- 
terion  may  also  be  used  to  test  which  variable  Xp+1  gives  the  best  discrimina¬ 
tion  when  added  to  the  set  (X^,  X2»...,X  ).  The  question  of  choosing  the  pro¬ 
per  set  of  variables  is  discussed  in  the  next  section. 


Stepwise  Discriminant  Analysis 

From  the  ten  variables  listed  in  Table  VI,  there  exist  1013  different 
subsets  (X^,  X  )  of  two  or  more  variables  for  which  the  discriminant 

function  could  be  calculated.  Clearly,  only  a  few  of  these  subsets  should  be 
examined.  Some  subsets  (e.g.,  short-period  P,  LK)  were  selected  for  testing 
because  they  are  well-known,  commonly-applied  discriminants;  others  were  select¬ 
ed  because  they  contain  significant  composite  information  from  a  single  seismic 
phase,  e.g.,  short-period  magnitude,  short-period  spectral  ratio,  and  complex¬ 
ity  from  P  waves.  However,  to  determine  the  subset  which,  in  accordance  with 
some  suitable  criterion,  is  the  "best"  subset  for  discrimination,  an  analytic 
procedure  described  below  is  employed. 


The  stepwise  discriminant  analysis  program  used  here,  BM07M  (Jennrich), 
1977),  is  a  "step-up"  procedure.  After  the  discriminant  analysis  has  been 
performed  for  a  subset  of  variables  (X^,  X^,...,  X^),  analysis  is  then  perform¬ 
ed  for  another  subset  obtained  by  adding,  and  sometimes  deleting,  a  variable 
in  the  next  step.  The  method  of  selecting  variables  to  be  added  or  deleted 
is  based  upon  the  partial  A-statistic,  defined  by  the  multiplicative  increment 


A'(X,u)  = 

A  (X) 


(20) 


in  Wilks'  A-criterion,  obtained  by  adding  the  variable  u  to  the  set  X  =  (X^, 
X2,...,  X  )  or  by  deleting  it  from  the  set  X=  (X^, 
set  the  corresponding  F-statistic  is  (Jennrich,  1977) 


■ » Xp , u ) . 


For  our  data 


F  =  (23-2-p)  1-A  (21) 

A'(X,u) 

This  statistic  is  known  as  either  the  "F-to-enter"  or  the  "F-to-remove"  statis¬ 
tic  depending  on  whether  the  variable  u  is  being  added  to  or  deleted  from  the 
set  of  variables  used  for  discrimination  in  the  previous  step.  The 


Jennrich,  R.  I.  (1977).  Stepwise  discriminant  analysis,  in  Statistical  Meth¬ 
ods  for  Digital  Computers,  K.  Enslien  et  al.  (eds.).  New  York,  NY  John 
Wiley  and  Sons. 
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numbers  of  degrees  of  freedom  are  1  and  21-p  for  F-to-enter  and  1  and  22-p  for 
F-to-remove.  To  move  from  one  step  to  the  next,  the  variable  with  the  smallest 
F-to-remove  is  deleted,  if  that  statistic  is  less  than  some  specified  thres¬ 
hold.  If  the  F-to-remove  values  are  such  that  all  variables  currently  entered 
are  retained,  then  that  one  of  the  remaining  variables  which  has  the  largest 
F-to-enter  is  added,  provided  that  the  F  value  is  higher  than  another  specified 
threshold  and  that  the  tolerance  is  greater  than  some  suitable  cut-off  value. 

(The  tolerance  is  given  by  one  minus  the  square  of  the  variable's  within-group 
multiple  correlation  with  the  variables  already  entered.)  As  a  result  of  this 
procedure  the  variables  are  ranked  at  each  step,  according  to  their  ability  to 
enhance  the  discrimination  capability  of  those  variables  already  found  signi¬ 
ficant.  Variables  which  become  superfluous  after  adding  others  are  lowered  in 
rank  and  deleted  if  their  informational  content  is  low  enough.  By  specifying 
satisfactory  values  of  the  F-to-enter  and  F-to-remove  thresholds,  one  may  find 
the  "best"  subset  of  variables;  that  is,  the  smallest  set  yielding  the  greatest 
group  separation  measured  by  Wilks'  A-criterion.  Since  the  degrees  of  freedom 
of  the  denominator  of  the  F-statistic  must  lie  in  the  range  21  to  12,  satisfac¬ 
tory  thresholds  of  F-to-enter  and  F-to-remove  were  somewhat  arbitrarily  taken 
to  be  3.0  and  4.0  respectively.  Note  that  the  procedure  for  selecting  the 
"best"  subset  of  variables  is  arbitrary  and  that  other  investigators  may  choose 
different  criteria  for  selecting  the  "best"  subset. 

Results  of  the  Stepwise  Discrimination  Analysis 

For  the  chosen  values  of  the  F-to-enter  and  F-to-remove  thresholds,  four 
of  the  ten  variables  were  entered  into  the  discrimination  algorithm:  first  motion, 
complexity,  short-period  P,  and  short-period  spectral  ratio  (see  Table  VII-A). 

For  this  subset  of  variables  the  lowest  a  posteriori  probability  of  correct 
classification  was  0.916,  which  was  found  for  one  of  the  two  earthquakes  without 
dilatational  first  motion  (Table  X).  As  expected,  examination  of  the  F-to-remove 
values  of  the  retained  variables  shows  that  first  motion  was  the  key  factor  in 
the  discrimination.  When  thresholds  were  set  low  enough  for  all  ten  variables 
to  be  entered  and  retained  (Table  VII-B) ,  the  a  posteriori  probabilities  were 
even  greater  than  before,  with  the  lowest  value  becoming  0.973.  Note  that 
this  situation  does  not  necessarily  mean  that  the  addition  of  six  other 
variables  would  significantly  enhance  discrimination  capability  outside  of  our 
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"training  set".  This  is  true  because  high  a  posteriori  probabilities  are,  to 
a  certain  extent,  an  artifact  of  the  large  number  of  variables  and  the  fact 
that  each  new  variable  is  forced  to  work  for  the  training  set.  The  F-statis- 
tic  approximation  to  Wilks'  A-criterion  actually  corresponds  to  a  greater 
separation  of  group  means  in  the  case  of  the  four  variables  than  in  the  case 
of  all  ten,  even  though  in  both  cases  F  greatly  exceeds  the  value  of  the  99.95 
percentage  point  for  the  appropriate  numbers  of  degrees  of  freedom. 

Again  note  that  the  determination  of  first  motions  was  frequently  rather 
subjective,  so  one  should  exercise  caution  in  accepting  the  validity  of  any 
discriminant  function  that  this  potentially  unreliable  variable  influences. 
Therefore,  the  stepwise  discriminant  analysis  was  repeated  using  the  same 
data  base  as  before,  but  this  time  omitting  the  tenth  (first  motion)  variable. 
Using  the  same  F-to-enter  and  F-to-remove  thresholds  as  before,  the  program 
first  entered  the  same  two  variables  previously  entered  immediately  after  the 
first  motion  — complexity  and  short-period  P —  but  then  it  added  Love  waves 
into  the  discriminant  function  and  deleted  complexity  (Table  VII1-A).  Thus, 
of  all  possible  subsets  of  the  nine  variables,  the  pair  of  short-period  P 
and  LQ  is  the  best  discriminant.  However,  since  the  F-to-enter  value  of  the 
short-period  spectral  ratio,  3.994,  was  so  close  to  the  arbitrarily  chosen 
threshold  of  4.00,  the  discriminant  analysis  was  repeated  using  a  threshold 
that  allowed  short-period  spectral  ratio  to  be  entered  as  the  third  variable. 
The  lowest  a_  posteriori  probabilities  were  quite  high  (Table  X) ,  but  as  in  the 
previous  cases  this  does  not  necessarily  indicate  a  significant  improvement 
in  separating  the  two  groups. 

The  lowest  a  posteriori  probability  found  in  both  the  three-  and  nine- 
variables  cases  was  that  of  event  E,  the  Denver  earthquake.  This  event  was 
deleted  from  the  data  base  and  the  three-variable  discriminant  function  was 
recomputed  on  the  basis  of  the  remaining  twenty-two  events.  When  this  new 
discriminant  was  applied  to  event  E,  the  probability  of  this  event's  being 
an  earthquake  was  calculated  to  be  0.703.  Thus,  the  SPP-LQ-SPSR  discrim¬ 
inant  was  capable  of  correctly  classifying  the  most  anomalous  event,  even 
when  that  event  was  not  used  in  the  training  set  to  compute  the  discriminant 
function. 
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log  comp  1.1098 
LPS  0.8812 


TABLE  VIII  (B) 
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Note  again  that  several  of  the  values  of  the  short-period  S  and  long- 
period  P  magnitudes  were  actually  upper  bounds  of  true  values  too  small  to 
be  detected  at  any  station.  To  judge  how  seriously  discrimination  was 
affected  by  the  use  of  these  upper  bounds,  a  seven-variable  discriminant  func¬ 
tion  was  calculated  which  did  not  use  short-period  S,  long-period  P,  or  first 
motion.  The  discrimination  capability  of  the  seven-variable  subset  was  about 
the  same  as  that  of  the  nine-variable  subset.  Therefore,  the  two  especially 
noisy  variables  did  not  contribute  any  significant  information  — or  misinfor¬ 
mation —  to  the  rest  of  the  data  base. 

The  magnitudes  of  the  Love  waves  could  not  always  be  measured  directly 
and  sometimes  had  to  be  inferred  from  Rayleigh-wave  magnitudes.  A  discrimi¬ 
nant  function  without  the  unreliable  Love-wave  data  had  to  be  examined,  and, 
so,  the  stepwise  discriminant  analysis  was  repeated,  suppressing  the  Love-wave 
and  first  motion  variables.  Table  IX  indicates  that  the  significant  variables 
in  this  case  were  short-period  P,  Rayleigh  waves,  and  short-period  spectral 
ratio.  This  particular  discriminant  correctly  classified  all  twenty-three 
events,  although  it  did  yield  a  rather  low  posteriori  probability  for  the 
Denver  earthquake. 

In  summary,  stepwise  discriminant  analysis  was  used  to  calculate  three 
different,  potentially  useful,  discriminant  functions.  The  first  of  these 
discriminants  amounts  to  little  more  than  the  criterion  that  events  exhibit¬ 
ing  dilatational  first  motions  are  earthquakes.  If  this  first-motion  crite¬ 
rion  is  judged  unreliable,  then  the  most  powerful  discriminant  is  one  based 
upon  short-period  P,  Love  waves,  and  short-period  spectral  ratio.  If  Love- 
wave  data  is  unavailable,  then  the  next-best  set  of  variables  is  short-period 
P,  Rayleigh  waves,  and  short-period  spectral  ratio.  Thus,  by  adding  SPSR  the 
discrimination  capability  of  the  classical  Mg  -  relation  is  somewhat  en¬ 
hanced,  and  the  discriminant  is  even  more  effective  when  M  is  measured  by 
Love  waves  rather  than  by  Rayleigh  waves.  Note,  however,  that  several  of 
the  twenty-three  events  studied  were  included  in  the  training  set  specifically 
because  of  their  failure  to  be  clearly  classified  on  the  basis  of  the  classi¬ 
cal  Mg  -  m^  criterion. 

Evaluation  of  Subsets  of  Discriminants 

Although  results  from  the  stepwise  analysis  revealed  which  subset  of 


SPS  0.0142 

LPP  0.0106 


variables  was  the  "best"  subset  for  discrimination,  it  is  useful  to  evaluate 
the  ability  of  certain  other  discriminants  to  correctly  classify  the  twenty- 
three  events  of  the  training  set.  The  results  are  given  below  of  attempts  at 
classification  where  these  other  subsets  of  discriminants  were  used. 

SPP,LR:  The  classical  Mg  -  m^  discriminant  correctly  classified  all  twen¬ 
ty-three  events,  although  the  a  posteriori  probability  was  rather  low  for  both 
event  E  (the  Denver  earthquake)  and  event  J  (BENHAM  aftershock).  Comparing  the 
results  of  this  discriminant  (Table  XI)  with  those  of  the  SPP-LR-SPSR  case  (Ta¬ 
ble  X)  shows  that  adding  the  short-period  spectral  ratio  significantly  widened 
the  separation  of  event  J  from  the  explosions  but  that  it  had  only  minimal  ef¬ 
fect  on  event  E  (Figure  26  indicates  this  result  was  expected).  Two  explosions, 
events  7  and  11,  also  had  low  a_  posteriori  probabilities.  Figure  20  shows  that 
although  these  events  lie  close  to  the  hypothetical  least-squares  regression  line 
passing  through  the  explosion  population,  they  also  lie  close  to  outlying  earth¬ 
quakes  that  are  widely  scattered  about  that  population  for  small  m^. 

SPP,  LPS:  This  discriminant  misclassif led  the  Baja  earthquake  and  strongly 
misclassif ied  RULISON.  Note  that  the  very  large  magnitude  of  LPS  for  RULISON 
(Figure  14)  is  suspect,  however,  because  this  value  was  based  upon  two  detec¬ 
tions  and  two  noise  levels,  and  the  two  detections  themselves  were  very  noisy, 
if  not  spurious.  Omitting  RULISON  from  the  training  set  would  lead  to  correct 
classification  by  this  discriminant  of  all  the  other  events,  including  the 
Baja  earthquake. 

SPP,  LPP:  For  the  twenty-three  events  as  a  whole,  this  was  a  poor  discri¬ 
minant,  misclassifying  four  events.  Six  of  the  events  in  the  training  set  had 
no  detectable  long-period  P  waves,  however,  and  were  characterized  by  upper 
bounds  for  LPP  (Figure  16).  When  these  six  events  were  deleted  from  the  ana¬ 
lysis,  this  discriminant  became  the  most  effective  one  tested;  even  the  out¬ 
liers  in  each  of  the  two  groups  were  widely  separated  from  the  other  group.  Thus 
SPP-LPP  is  a  valuablr  discriminant  when  long-period  P  waves  can  be  detected,  but 
note  that  LPP  cannot  be  measured  for  explosions  with  <  5.2. 

SPP,  SPS:  For  the  full  data  set,  including  the  nine  upper  bounds,  this 
pair  of  variables  was  useless  as  a  discriminant  (Figure  13).  However,  when  con¬ 
sidering  only  the  fourteen  events  actually  detecting  SPS,  this  discriminant  mis- 
classified  only  RULISON,  which  was  misclassif ied  strongly.  However,  like  the 
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TABLE  X  (Continued) 


_A  Posteriori  Probabilities — Stepwise  Discriminant 
*Test #  Discriminant  Tested 


1  F-to-Enter  =  4.00,  F-to-Remove  =  3.00 
(SPP,  log  SPSR,  log  comp,  1st  mtn) 

2  F-to-Enter  =  0.01,  F-to-Remove  =  0.005 
(all  10  variables) 

3  no  1st  mtn,  F-to-Enter  =  4.00,  F-to-Remove 
(SPP,  LQ) 

4  SPP,  LQ,  log  SPSR 

5  no  1st  mtn,  F-to-Enter  =  0.01,  F-to-Remove 
(all  other  9  variables) 

6  SPP,  LQ,  log  SPSR;  22  events  (no  event  E) 

7  7  variables  (no  SPS,  LPP ,  1st  mtn) 

8  no  LQ  or  1st  mtn,  F-to-Remove  =  3.00 
(SPP,  LR,  log  SPSR) 
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Analysis 


3.00 


0.005 


TABLE  XI  (Continued) 


A  Posteriori  Probabilities — Selected  Discriminants  Using  Network  Estimates 
*Test  #  Discriminant  Tested 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


SPP,  LR 
SPP,  LPS 
SPP,  LPP 

SPP,  LPP  (no  noise  levels) 

SPP,  SPS 

SPP,  SPS  (no  noise  levels) 

LR,  log  LPSR 

SPP,  log  SPSR 

SPP,  log  comp 

log  SPSR,  log  comp 

SPP,  log  SPSR,  log  comp 

SPP,  LQ,LR 

SPP,  LQ  +  LR 

LQ  LR 

SPP,  SPP 

- - -  .  (SPP,  LQ,  LR) 

'SPP2'"+  LQ2  +  LR2 
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case  of  long-period  S,  the  RULISON  SPS  magnitude  is  uncertain  because  it  is 
based  upon  three  noise  levels  and  a  single  noisy  detection.  If  this  borderline 
"detection"  had  been  considered  a  noise  measurement,  then  this  variable  pair 
would  have  correctly  classified  the  remaining  thirteen  events.  This  discrimi¬ 
nant  may  be  useful  for  events  whose  magnitude  (m^)  is  large  enough  for  detec¬ 
tion  of  SPS.  Figure  13  suggests  m^  >  5  1/2  for  explosions  and  m^  >  5  for 
earthquakes  for  SPS  detection. 

LR,LPSR:  This  pair  of  variables  yields  no  discrimination.  Thus  long- 
period  surface  waves  alone  are  insufficient  for  discrimination,  as  expected 
on  theoretical  grounds. 

SPP.SPSR  and  SPP,  complexity;  These  two  pairs  of  variables  are  ineffec¬ 
tive  as  discriminants.  The  former  misclassified  three  events,  and  the  latter 
misclassif led  five  events.  When  combined  to  form  the  triad  SPP-SPSR-complex- 
ity,  however,  they  misclassified  only  the  Borrego  Mountain  earthquake  (although 
the  Caliente  earthquake  is  a  marginal  case).  This  three-variable  discriminant 
is  noteworthy  because  it  used  only  short-period  P-wave  data. 

SPP,  LQ,  LR;  This  discriminant  was  tested  in  two  different  ways.  In 
the  first  test  the  discriminant  function  for  these  three  variables,  calculated 
in  the  manner  previously  described,  was  found  to  be 

13.098-SPP  -  5.985*LQ  -  4.674’LR  -  18.122  =  0.0 

In  the  second  test  the  sum  of  the  Rayleigh-wave  and  Love-wave  magnitudes  was 
treated  as  one  variable  (Figure  23),  and  the  discriminant  function  was  found 
to  be 

10.879'SPP  -  4.488(LQ  +  LR)  -  14.626  =  0.0 

Although  controversy  may  arise  over  what  is  physically  more  meaningful,  the 
individual  values  of  LQ  and  LR  or  their  sum,  there  is  still  value  in  treating 
the  two  surface-wave  magnitudes  as  one  variable  on  account  of  the  facility 
with  which  a  discriminant  line  in  a  plane  may  be  graphically  manipulated.  A 
discriminant  plane  in  a  three-dimensional  space  must  be  treated  by  analytic 
methods  alone.  One  would  ordinarily  expect  the  transformation  from  the  three- 
dimensional  space  to  the  plane  to  involve  a  loss  of  information,  but  since  LQ 
versus  LR  is  itself  a  poor  discriminant,  the  pair  SPP,  LQ  +  LR  gives 


satisfactory  results  when  compared  with  those  of  the  triad  SPP,LQ,LR.  Because 
Mg  -  was  a  poorer  discriminant  for  the  data  set  than  the  SPP,LQ  pair  (Table 
X) ,  neither  of  the  two  discriminants  involving  LR  gives  results  quite  go  good 
as  SPP,LQ  alone.  Again  note  that  poor  separation  on  Mg  -  m^  plots  was  one 
criterion  for  the  inclusion  of  several  of  the  events  in  the  training  set. 


The  values  of  F-to-enter  and  F-to-remove  for  short-period  P,  which  were 
calculated  in  the  stepwise  experiments,  indicate  that  various  discriminant 
functions  gave  significance  to  the  difference  in  the  mean  of  m^  for  the  two 
populations  (m^  =  5.72  for  explosions  and  5.00  for  earthquakes).  Thus,  m^  was 
acting  as  a  weak  discriminant  by  itself.  Since,  for  the  purposes  of  discrim¬ 
ination,  what  is  important  is  not  the  magnitude  of  short-period  P  itself  but 
rather  its  relation  to  the  other  variables,  two  discriminants  were  tested  that 
were  in  some  sense  "normalized",  to  correct  for  m^  variations.  In  the  first 
test  the  magnitudes  of  LQ  and  LR,  each  divided  by  m^,  were  used  as  a  discrim¬ 
inant  pair.  This  test  resulted  in  four  misclassifications.  The  second  test 
was  a  three-variable  discrimination  experiment  in  which  the  magnitudes  of 
SPP,LQ,  and  LR  were  each  divided  by  the  square  root  of  the  sum  of  their  squares. 
This  test  is  equivalent  to  the  earlier  three-variable  discriminant  except  that 
now  each  of  the  twenty-three  data  points  lies  at  unit  distance  from  the  origin 
of  the  three-dimensional  space.  This  test  resulted  in  one  misclassif ication. 


Thus,  discrimination  based  on  magnitude  ratios  such  as  etc.,  is  less 


effective,  perhaps  justifiably,  than  discrimination  based  on  the  magnitudes 
themselves. 


If  the  SPP,LQ  and  SPP,LR  discriminant  lines  have  the  form  LQ  =  a«SPP  +  b 
and  LR  =  c'SPP  +  d,  then  the  earthquake  and  explosion  populations  will  lie 
within  different  regions  on  a  plot  of  (LQ  -  a*SPP)  versus  (LR  -  c*SPP).  Be¬ 
cause  the  slopes  a  and  c  are  close  to  unity,  a  useful  discriminant  not  affected 
by  m^  differences  in  the  two  populations  is  (LQ-SPP)  versus  (LR-SPP).  For  our 
data  set  this  discriminant  (not  tabulated)  resulted  in  no  misclassifications. 

Single-Station  Discrimination 

Several  of  the  discrimination  tests  that  had  been  performed  using  event 
magnitudes  from  network  data  were  repeated  using  data  from  only  one  station. 
Since  more  data  were  available  from  NP-NT  and  RK-ON  than  from  the  other 
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stations,  data  from  each  of  these  stations  were  used  separately  in  an  attempt  to 
evaluate  the  effectiveness  of  single-station  discrimination  functions.  Tables 
XII  and  XIII  indicate  that  there  were  many  gaps  in  the  available  data,  even 
for  those  two  most  complete  station  records.  Some  discriminants  could  not  be 
evaluated  for  certain  events  because  of  missing  data  (expecially  due  to  over¬ 
driven  long-period  recordings  and  unavailable  digital  data) ;  and  others  were 
inapplicable  because  of  failure  to  detect  all  variables  used  in  the  discrim¬ 
inant  function,  especially  long-period  P  and  short-period  S. 

Results  of  the  single-station  tests  (Tables  XII  and  XIII)  show  that  even 
those  discriminants  that  worked  well  in  the  previous  experiments  now  misclas- 
sified  a  number  of  events.  The  large  scatter  of  single-station  measurements 
of  event  magnitudes  about  the  "true"  values  (i.e.,  the  network-based  estimates) 
causes  significant  overlap  of  the  earthquake  and  explosion  populations  for  any 
subset  of  the  ten  variables.  Thus,  discrimination  based  upon  a  network  is 
significantly  more  powerful  than  discrimination  based  on  a  single  station 
because:  1)  the  training  set  data  base  is  more  nearly  complete;  and  2)  the 
magnitude  estimates  are  more  accurate,  hence  the  population  variances  are  smaller. 
The  question  remains  of  whether  a  method  of  handling  network  data  exists  which 
is  preferable  to  the  method  already  outlined.  A  case  could  be  made  that  aver¬ 
aging  measurements  from  the  entire  network,  eiLher  by  the  conventional  arith¬ 
metic  mean  or  by  the  maximum- likelihood  method,  blurs  out  station  effects  and 
source-path  effects  which  might  affect  discrimination.  Therefore,  an  alter¬ 
nate  process  is  examined  in  which  discriminant  functions  are  calculated  indi¬ 
vidually  at  each  station  and  are  then  combined  to  perform  classification  based 
upon  the  entire  network. 

Multiple-Station  Discriminant  Functions — Theory 

In  calculating  the  discriminant  functions  detailed  thus  far,  we  assumed 
that  the  covariance  matrices  of  the  earthquake  and  explosion  populations  were 
identical.  To  calculate  multiple-station  discriminant  functions,  we  further 
assume  that  the  covariance  matrices  are  the  same  for  the  data  measured  at  each 
network  station.  This  common  covariance  E  is  approximated  by  the  sample  pop¬ 
ulation  covariance  S  obtained  by  considering  the  pooled  data  from  the  entire 
network.  In  order  to  obtain  S  it  is  necessary  to  reject  all  incomplete  data 
vectors,  i.e.,  data  for  a  given  event  and  a  given  station  for  which  one  or 
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TABLE  XII 


I 


TABLE  XII  (Continued) 

A  Posteriori  Probabilities — Selected  Discriminants  Using  NP-NT  Estimates 


*Test  if  Discriminant  Tested 


1  SPP,  LQ 

2  SPP,  LR 

3  SPP,  LPS 

4  SPP,  SPS 

5  LR,  log  LPSR 

6  SPP,  log  SPSR 

7  SPP,  log  comp 

8  log  SPSR,  log  comp 

9  SPP,  log  SPSR,  log  comp 

10  SPP,  LQ,  LR 
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TABLE  XIII  (Continued) 


A  Posteriori  Probabilities — Selected  Discriminants  Using  RK-ON  Estimates 


*Test  #  Discriminant  Tested 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


SPP,  LQ 

SPP,  LR 

SPP,  LPS 

SPP,  LPP 

SPP,  SPS 

LR,  log  LPSR 

SPP,  log  SPSR 

SPP,  log  comp 

log  SPSR,  log  comp 

SPP,  log  SPSR,  log  comp 

SPP,  LQ,  LR 
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more  of  the  variables  used  in  the  discriminant  could  not  be  measured,  either 
due  to  missing  data  or  to  noise  thresholds.  It  is  to  be  expected  that  the 
value  of  S  calculated  in  this  manner  from  the  pooled  data  is  less  accurate 
than  the  value  which  was  calculated  from  the  network-averaged  data,  because: 

1)  many  measurements  of  signal  amplitudes,  spectral  ratios,  and  complexities 
had  to  be  deleted  from  the  calculation  since  they  were  components  of  incomplete 
data  vectors;  and  2)  it  was  impossible  to  utilize  the  maximum-likelihood  esti¬ 
mator  for  the  magnitudes  of  weak  signals.  (An  alternative  procedure  that  also 
was  tested  would  be  to  use  the  same  S  matrix  which  previously  was  calculated 
for  the  network  discriminant  case.)  Although  the  covariance  at  each  station  is 
the  same,  the  population  means  are  not;  so  for  the  m-th  station  of  the  sta¬ 
tion  network,  the  means  of  the  earthquake  and  explosion  populations  as  meas¬ 
ured  at  that  station  are  calculated  to  be  X^m  and  Equation  (14)  and  (lb) 

indicate  the  discriminant  function  and  classif ication  function  for  the  m-th 
station  alone  to  be 


and 
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Jm 
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j  =  1,2;  m  ■  »Ng“ 


(22) 

(23) 


Using  these  functions,  the  a  posteriori  probability  P-j^X^)  that  the  k-th 
event  of  the  training  set  belongs  to  the  j-th  group,  as  determined  by  the  m-th 
station,  can  be  calculated,  where  denotes  the  data  vector  of  the  k-th  event 
measured  at  the  m-th  station.  However,  there  is  more  interest  in  classifica¬ 
tion  probabilities  p^ (X^)  determined  by  the  entire  network,  where  X^  denotes 
the  vector  of  all  observations  of  the  k-th  event  (X^,  X^  ).  In 

order  to  determine  these  probabilities,  the  discriminant  function  is  employed 
for  the  network  as  a  whole,  which  is  simply  (Shumway  and  Blandford,  1972) 


«  v 


N 


*  *1™ 


N 

+ 


(24) 


ra=l  m“l 

To  classify  a  given  event,  then,  take  the  sum  of  the  values  of  the  individual- 
station  discriminant  functions  applied  to  that  event  and  classify  it  as  an 
explosion  or  an  earthquake  based  upon  whether  that  sum  is  positive  or  negative. 
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By  introducing  the  classification  functions  for  the  network,  as  given  by 


N 

£S 

m“l 


N 

+  I8 
m-1 


j  *  1,2 


(25) 


the  a  posteriori  probabilities  that  an  event  is  an  explosion  or  an  earthquake, 
as  determined  by  the  network,  are  obtained  from  equations  13  and  12  as 


P,  (X.  )  -  exp[dl(*k)] _ _  _  1  _ 

expLdj_(X^)J  +  exp[d2(XJp]  1  +  exp[-d(X^)] 

Po(X.)  _  exp[d2(xk)  ] _ _ 1  _ 

exp  [d^O^p  ]  +  exp[d2(X^)]  1  +  exp[d(X^)  ] 


It  should  be  known  whether  the  individual  station  discriminant  functions  are 
less  reliable  for  some  stations  than  for  others  before  combining  them  to  form 
the  network  discrimination  functions.  Unreliable  stations  could  then  be 
deleted  from  the  sum  over  the  network  (equation  24).  To  assess  station  reli¬ 
ability,  a  criterion  for  estimating  the  probability  of  misclassif ication  by  a 
single  discriminant  function  must  be  examined. 


For  the  discriminant  function  with  known  parameters,  given  by  (6)  and 
(7),  the  probability  of  misclassifying  an  explosion  as  an  earthquake  is  given  by 

P  (2  1 1)  =  H-  |)  (28) 

2 

where  $  (Z)  denotes  the  normal  distribution  function  and  6  is  the  Mahalanobis 
distance  (Morrison,  1976),  defined  as 

62  =  cq  -T2r  f1  (7X  -T2)-  (29) 

Since  the  parameters  are  unknown  in  practice'',  one  must  estimate  P ( 2  1 1 )  by 

1>(2|1)  =  *(-§)  (30) 

where 

D2  =  (  -  -x2r  S  1  (  ^  -  ~x2)  •  (31) 


Morrison,  D.  F.  (1976).  Multivariate  Statistical  Methods.  McGraw-Hill 
Publ.  Co.,  New  York. 
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Dunn  (1971)  and  McLachlan  (1974)  demonstrated  the  inadequacy  of  P (2 1 1 )  and 
investigated  the  conditions  under  which  other  functions  are  preferable  as 
estimators  of  P ( 2 | 1) .  For  the  purpose  of  this  study,  however,  the  signifi¬ 
cant  result  is  not  the  exact  value  of  the  misclassif ication  probability  cal¬ 
culated  by  means  of  equation  (30)  or  some  other  estimator,  but  rather  the 

2 

simple  criterion  that  discriminant  functions  for  which  D  is  small  have  a 

2 

higher  misclassification  probability  than  do  those  for  which  D  is  large. 

The  Mahalanobis  distance  for  multiple-station  discrimination  is  given 
2 

by  the  sum  of  D  for  each  of  the  Ng  single-station  discriminant  functions. 

Therefore,  one  should  evaluate  (30)  for  each  station  in  the  network  and  iden- 

2 

tify  stations  with  the  smallest  D  as  stations  contributing  the  least  to  the 

network  discriminant  (24)  and  a  posteriori  probabilities  (26,27).  Stations 

2 

with  sufficiently  small  D  may  than  be  deleted  from  the  analysis.  Identifi¬ 
cation  of  stations  that  yield  poor  discrimination  for  a  given  subset  of  var¬ 
iables  may  also  relate  valuable  information  about  the  source-path  effects  of 
these  variables. 

Evaluation  of  Selected  Multiple-Station  Discriminant  Functions 

Several  particular  subsets  of  variables  tested  as  discriminants  in  the 
network  and  single-station  experiments  were  also  tested  as  multiple-station 
discriminants.  The  test  results  showed  that  the  subsets  which  were  most 
effective  as  network  discriminants  were  also  the  most  effective  as  multiple- 
station  discriminants.  However,  in  using  multiple-station  discriminant  func¬ 
tions,  some  problems  did  develop  because  of  missing  data  or  because  noise 
levels  exceeded  weak  signal  amplitudes;  the  problems  reduced  the  number  of 
stations  detecting  the  chosen  variables.  Thus,  few  single-station  discriminants 
were  used  to  calculate  the  network  discriminant  for  a  given  event,  and  the 
single-station  discriminants  themselves  were  calculated  from  the  data  of  only 
a  few  events. 

For  the  ten  selected  discriminants,  a.  posteriori  probabilities  for  each  event 
2 

and  both  D  and  the  number  of  detections  for  each  station  are  listed  in 


Dunn,  0.  J.  (1971).  Some  expected  values  for  probabilities  of  correct  clas¬ 
sification  in  discriminant  analysis.  Technometrics ,  13,  345-353. 

McLachlan,  G.  J.  (1974).  Estimation  of  the  errors  of  misclassification  on  the 
criterion  of  asymptotic  mean  square  error.  Technometrics ,  16,  255-260. 
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Table  XIV.  Note  the  number  of  detections  when  evaluating  the  reliability 
of  a  single-station  discriminant  function.  This  is  necessary  because  sta¬ 
tions  detecting  only  a  few  events  from  either  or  both  of  the  populations 

2 

would  likely  measure,  by  means  of  D  ,  a  large  separation  between  the  two 

groups.  Thus,  while  such  stations  seem  to  have  reliable  discriminant  func- 
2 

tions  as  far  as  D  is  concerned,  in  reality  they  are  probably  less  reliable 

2 

than  stations  with  lower  values  of  D  whicli  have  more  detections.  If  one 

such  station  strongly  misclassif ies  a  certain  event,  then  the  sum  (24)  may 

also  classify  it  incorrectly.  Note,  for  instance,  that  in  the  case  of  the 

2 

common  Mg  -  m^  discriminant,  D  is  low  for  NP-NT  and  WH-YK.  If  the  analysis 

is  repeated  omitting  data  from  these  two  stations,  the  misclassif ication  of 

RULISON  is  eliminated  but  also  the  only  station  detecting  both  Mg  and  m^  for 

the  Borrego  Mountain  earthquake  is  deleted.  The  reliability  of  BE-FL  is  also 

questionable  because  the  observed  explosion  population  consisted  of  only  one 

event.  The  M  -  m,  discriminant  functions  for  the  seven  stations  are  shown 
s  b 

in  Figure  29. 

By  eliminating  all  data  vectors  incomplete  because  some  variables  used 
in  a  given  discriminant  function  were  undetected,  a  bias  is  introduced  into 
the  pooled  data  used  to  compute  the  covariance  matrix.  Therefore,  the  multi¬ 
ple-station  discriminant  analysis  was  repeated  using  instead  the  covariance 
matrices  calculated  for  the  network  maximum-likelihood  magnitude  estimates. 
While  this  process  seemed  to  improve  the  results  somewhat  for  discriminants  not 
involving  digital  data  (Table  XV),  the  discriminants  involving  SPSR,  LPSR, 
and/or  complexity  behaved  erratically,  a  situation  expected  because  the  maxi¬ 
mum-likelihood  estimator  for  these  three  variables  also  deleted  noisy  observa¬ 
tions.  The  principle  difference  between  the  two  estimates  of  the  covariance 
for  the  digital  data  is  that,  for  the  maximum-likelihood  estimator,  means  are 
computed  for  each  event  and  then  the  covariance  is  computed  using  those  means 
with  equal  weight,  whereas  for  the  pooled  data,  the  covariance  is  computed 
directly  from  the  individual  observations.  While  the  latter  method  is  pro¬ 
bably  preferable  for  complexity  and  spectral  ratios,  the  maximum-likelihood 
estimates  may  yield  more  accurate  covariances  for  the  other  variables.  For 
evidence  of  this  accuracy,  note  that  the  Mg  -  m^  discriminant  correctly  clas¬ 
sifies  RULISON  when  the  covariance  matrix  from  the  network  case  is  used. 
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Figure  30  Single-station  Mg-m^  discriminant  functions  with  covariance 
matrices  computed  from  network  estimates  of  event  magnitudes. 


Also  note  that  the  value  of  D  was  greater  for  every  station  when  the  network 
covariance  was  applied  to  each  of  the  discriminants  not  using  digital  data. 
The  effect  of  using  the  network  covariance  in  the  Mg  -  discriminant  func¬ 
tions  is  shown  in  Figure  30. 


Summary  of  the  Three  Experiments 

In  the  first  experiment,  the  maximum-likelihood  estimator  was  used  to 
calculate  network  averages  of  each  variable's  magnitude  for  each  event,  exclud¬ 
ing  spectral  ratios  and  complexity.  A  process  of  stepwise  discriminant  anal¬ 
ysis  selected  the  most  effective  subsets  of  the  ten  variables.  The  discrimin¬ 
ation  capabilities  of  these  and  certain  other  subsets  were  evaluated  according 
to  three  criteria:  1)  the  number  of  training  set  events  which  the  discriminant 
misclassif led;  2)  the  a  posteriori  probabilities,  which  measured  how  strongly 
a  given  event  was  classified  into  one  of  the  two  populations;  3)  the  F-approxi- 
mation  to  Wilks'  A-statistic,  which  measured  how  well  the  discriminant  separa¬ 
ted  the  means  of  the  two  groups.  The  polarity  of  the  first  motion  was  found 
to  be  the  best  discriminant.  The  results  also  showed  that  the  effectiveness 
of  the  classical  Mg  -  m^  discriminant  was  enhanced  when  SPSR  was  added  as  a 

third  variable  and  that  it  was  enhanced  still  further  when  M  was  measured 

s 

from  LQ  rather  than  from  LR.  Other  discriminant  pairs  were  found  effective, 
though  only  when  applied  to  events  large  enough  to  permit  the  magnitudes  of  all 
variables  to  be  determined.  Among  them  were  SPP-SPS,  SPP-LPS,  and,  particu¬ 
larly,  SPP-LPP.  The  discriminant  SPP-SPSR-complexity  proved  of  some  value  in 
the  absence  of  long-period  data.  Results  also  showed  that  forming  ratios  of 
magnitudes  gave  poorer  discrimination  than  the  magnitudes  themselves.  Finally, 
the  A-statistic  indicated  a  large  separation  of  population  means  for  almost 
every  discriminant,  even  those  which  misclassif ied  some  events.  Thus,  the 
percentage  point  corresponding  to  the  value  found  for  the  F-statistic  tends  to 
be  misleadingly  high  in  its  evaluation  of  a  discriminant's  effectiveness. 


The  network  discrimination  method  has  two  advantages  over  single-station 
discrimination.  First,  since  the  magnitude  of  each  variable  is  calculated 
as  an  event  mean,  incomplete  observations  may  be  retained  in  the  data  set. 

For  example,  if  the  discriminant  in  question  is  Mg  -  m^,  then  a  station  meas¬ 
urement  of  m^  may  be  used  in  the  calculation  of  the  mean  m^  for  that  event 
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despite  the  absence  of  data  for  Mg  from  that  station.  Second,  because  event 
magnitudes  based  upon  network  detections  are  more  accurate  than  are  those  based 
on  individual  station  detections,  the  scatter  of  events  in  the  two  populations 
is  reduced  and  there  is  less  overlap.  This  is  particularly  true  of  events  for 
which  some  stations  were  unable  to  measure  all  the  variables  on  account  of 
noise  and  for  which  magnitudes  were  then  calculated  by  means  of  the  maximun- 
likelihood  estimator. 

The  principal  disadvantage  of  the  network  method  was  that  source-path 
effects  and  station  effects  were  ignored,  so  the  magnitudes  of  some  events 
may  have  been  systematically  influenced  by  discordant  values  from  particular 
stations.  Another  disadvantage  was  that  a  different  number  and  configura¬ 
tion  of  stations  was  usually  used  to  estimate  the  means  of  different  vari¬ 
ables  for  a  given  event. 

The  second  experiment  consisted  of  evaluating  certain  previously  tested 
discriminants  using  data  only  from  NP-NT  or  RK-ON.  Results  showed  that  no 
discriminant  was  very  effective  because  for  many  events  no  one  station  detected 
all  necessary  variables.  Thus,  the  size  of  the  training  set  was  reduced,  and 
the  event  magnitudes  retained  were  unreliable,  causing  the  two  populations  to 
overlap.  The  single-station  discriminant  functions  were  more  poorly  determined 
in  this  experiment  than  in  the  multiple-station  experiment  because  in  the  lat¬ 
ter  case  the  covariance  matrices  were  calculated  with  pooled  data  from  all  sta¬ 
tions  rather  than  from  each  station's  data  alone. 

The  third  experiment  used  the  covariances  of  pooled  data  to  calculate 
single-station  discriminant  functions  which  were  summed  to  give  a  discriminant 
function  for  the  entire  network.  Results  showed  that,  for  discriminants  not 
involving  spectral  ratios  or  complexity,  covariances  calculated  using  the  meth¬ 
od  of  the  first  experiment  were  preferable  to  those  based  upon  pooled  data. 

For  the  most  part,  the  effective  discriminants  were  the  same  as  those  in  the 
first  experiment.  However,  there  were  more  misclassif ications  by  multi-station 
discriminants  than  by  those  based  on  network  averages.  In  spite  of  this  situ¬ 
ation,  average  a_  posteriori  probabilities  were  sometimes  higher  for  the  multi¬ 
ple-station  case,  indicating  that  this  method  tends  to  classify  — or  misclas- 
sify —  events  strongly.  A  bad  observation  or  a  bad  discriminant  function  for  a 
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single  station  can,  therefore,  strongly  affect  the  classification  of  an  event. 

The  multiple-station  discriminant  function  has  at  least  two  advantages. 

First,  single-station  discriminants  can  be  ranked  according  to  their  values 
2  2 

of  D  ,  and  stations  with  low  D  can  be  deleted  from  the  data  base.  One  might 
thus  accommodate  certain  source-path  effects  by  using  multiple-station  discri¬ 
minants.  Second,  by  not  separately  computing  network  averages  of  each  vari¬ 
able  and  by  retaining  only  those  observations  which  detected  all  the  necessary 
variables,  a  uniform  data  base  is  used  for  each  of  the  variables  for  a  given 
event.  Discarding  stations  with  incomplete  observations  creates  a  smaller 
volume  of  data  than  was  retained  in  the  first  experiment;  nevertheless,  for 
a  given  event,  the  means  for  each  variable  are  computed  on  the  basis  of 
observations  from  the  same  set  of  stations,  and  thus  one  may  avoid  possible 
biases  in  the  means  of  different  variables  stemming  from  different  station 
configurations  and  different  variances  due  to  nonuniform  sample  sizes. 

The  problem  of  missing  data,  already  discussed  in  connection  with  single¬ 
station  discrimination  functions,  causes  many  events  to  be  deleted  from  the 
training  set  for  many  stations,  resulting  in  poorly  determined  discriminant 
functions.  The  scatter  in  the  magnitudes  is  also  large  for  individual  sta¬ 
tions,  and  the  maximum-likelihood  estimator  of  event  magnitudes  is  inappli¬ 
cable.  In  the  final  analysis,  therefore,  discrimination  based  upon  network 
magnitude  estimates  is  preferable  to  that  based  on  multiple-station  discrimi¬ 
nant  functions  unless:  1)  there  is  a  large  data  base  of  complete  observations 
for  events  in  both  populations  of  the  training  set  for  each  station;  2)  the 
magnitudes  measured  at  each  station  are  not  widely  scattered  around  the  true 
values;  3)  the  magnitudes  are  large  enough  that  a  bias  is  not  introduced  by 
the  rejection  of  data  vectors  which,  because  of  noise,  did  not  include  certain 
variable(s) . 

Recommendations  on  Classification  of  Unknown  Events 

While  applying  discriminant  functions  to  events  whose  classification  is 
unknown  a  priori  is  a  simple  process,  it  requires  that  an  adequately  large 
training  set  of  known  explosions  and  earthquakes  be  in  the  same  source  region 
as  the  unknown  events.  For  both  known  and  unknown  events,  magnitudes  of  the 
short-period  and  long-period  P  waves  and  S  waves  and  those  of  the  surface  waves. 
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LQ  and  LR,  should  be  measured,  and  the  short-period  spectral  ratio  should  be  calcu¬ 
lated.  (Discrimination  based  upon  polarity  of  first  motions  can  be  carried  out 
without  reference  to  any  discriminant  function  or  training  set).  Ringdal's 
method  should  be  used  to  compute  event  magnitudes  using  noise  levels  at  the  ex¬ 
pected  arrival  times  for  all  the  stations  in  the  network  that  do  not  detect 
the  signal.  The  training  set  should  be  used  to  compute  these  seven  discrimi¬ 
nant  functions:  SPP-SPS,  SPP-LPP,  SPP-LPS,  SPP-LQ,  SPP-LR,  SPP-LQ-SPSR,  and 
SPP-LR-SPSR.  Since  the  relative  effectiveness  of  these  discriminants  may  not 
be  the  same  for  events  outside  the  Western  United  States  as  for  events  in  this 
study,  the  number  of  a  posteriori  event  misclassif ications  for  each  discrimi¬ 
nant  in  the  training  set  should  be  noted  so  that  arbitrary  weights  can  be  as¬ 
signed  to  the  "best"  discriminant  functions.  Classification  based  on  the  seven 
weighted  discriminants  (or  as  many  as  are  applicable  to  the  given  region)  is 
assumed  more  reliable  than  classification  based  on  any  one  discriminant  alone. 

Each  of  the  seven  discriminant  functions  can  be  evaluated  for  the  data 
vector  consisting  of  magnitudes  measured  for  a  given  unknown  event.  For  the 
unknown  events,  however,  we  introduce  the  a  priori  probabilities  h  and  1-h 
for  an  explosion  or  an  earthquake  respectively  and  loss  function  £(i,j)  for 
classifying  an  event  from  the  j-th  population  as  one  from  the  i-th.  The 
criterion  for  discrimination,  formerly  equation  (8),  becomes  (Morrison,  1976) 

px(  T )  >  p2<  T  )  if  T'D~+c>ioge  t-(h~^(2.i)?2)  1  .  (32) 

In  addition,  the  appropriate  changes  are  made  in  the  computation  of  the  clas¬ 
sification  probabilities  given  by  (13). 


The  final  step  in  the  classification  process  is  to  take  the  sum,  with 
arbitrary  weights  assigned  to  the  "best"  discriminants,  of  the  classification 
probabilities  calculated  for  each  of  the  N  discriminants  (N  <_  7)  and  compute 
the  mean  probabilities 
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One  may  assume  for  small  events  that  discriminants  SPP-SPS,  SPP-LPP,  and  pos¬ 
sibly  SPP-LPS  will  be  inapplicable.  If  an  event  is  so  small  that  no  surface 
waves  can  be  measured,  then  it  may  be  classified,  though  not  confidently, 
using  the  discriminant  SPP-SPSR-complexity . 


CONCLUSION 


Results  for  seismic  discrimination  of  events,  selected  from  the  South¬ 
western  United  States,  were  in  good  agreement  with  theoretical  considerations 
and  with  past  studies  on  other  event  groups  in  this  region.  Specifically, 
agreement  existed  between  expectations  and  results  for: 

1)  SP-LR.  Theory  predicts  that  only  a  few  tenths  separate  smaller  mag¬ 
nitude  events  on  Mg  -  m^  plots  if  the  earthquakes  are  shallow.  Our 
data  included  three  that  were  only  marginally  separated  from  explo¬ 
sions.  Two  of  these  events  were  strike-slip  (Baja  California  and 
Denver  earthquakes)  and  the  third  was  a  mixture  of  strike-slip  and 
dip-slip  (BENHAM  aftershock). 

2)  SPS-SPP  or  LPS-LPP .  Shear  waves  from  explosions  are  predicted  to  be 
relatively  small.  Our  data  shows  that  short-period  or  long-period 
S-wave  amplitude  is  a  good  discriminant  if  it  can  be  measured  and 
that  it  should  definitely  be  included  in  multiple-discrimination 
analysis. 

3)  SPP-LQ.  For  LQ  as  for  shear  waves,  theory  and  past  experience  predict 
relatively  less  generation  for  explosions.  Our  data  showed  that  LQ, 
if  measurable,  enhances  discrimination  because  explosions  have 
smaller  LQ  amplitudes  and  earthquakes  often  have  large  LQ  when 

LR  is  small. 

4)  LPP-SPP.  Theory  predicts  cancellation  of  LPP  for  shallow  explosions. 
Data  in  this  study  supports  this  theory,  showing  that  LPP,  if  measur¬ 
able,  is  a  good  discriminant. 

5)  Corner  frequency  of  P-wave  spectrum.  Because  of  more  prolonged 
source  time  functions,  earthquakes  will  have  lo'-er  corner  frequen¬ 
cies  for  a  given  seismic  moment.  This  discriminant  was  studied  at 
only  two  stations,  but  results  were  still  encouraging.  However, 
much  higher  P-wave  S/N  ratio  is  required  than  for  m^  measurement  in 
order  to  estimate  these  spectral  parameters. 

6)  P-wave  spectral  ratio.  Theory  predicts  no  clear  results  for  the  band 
used  in  this  study  —  roughly  0.5  -  2.0  Hz.  Likewise,  our  data  showed 
little  or  no  separation  between  earthquakes  and  explosions,  even 
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after  application  of  presumed  attenuation  corrections  appropriate  to 
each  event's  source  area. 

7)  Complexity  of  P-wave  signal.  Only  small  enhancement  of  the  coda 
level  is  possible  because  of  the  deeper  crustal  focus  of  earthquakes. 
For  earthquakes  with  simple,  short-duration  rupture  mechanisms,  no 
increase  in  the  coda  level  should  occur.  Although  data  from  our  event 
set  showed  that  complexity  was  of  some  discrimination  value,  it  is 
highly  path-receiver  dependent,  making  it  unreliable. 

8)  Rayleigh-wave  spectral  ratio.  This  parameter,  in  the  .02  -  .05  Hz 
band  used  here,  is  known  to  vary  in  theory  over  a  larger  range  of 
values  for  earthquakes  than  for  explosions.  Our  data  supported 
this  prediction,  and  no  separation  was  found  between  the  two  source 
types . 

In  addition  to  these  discriminants,  first  motion  was  studied.  In  theory, 
first  motion  is  an  almost  infallible  method  of  identifying  earthquakes,  if 
clear  dilatations  are  present.  Although  results  from  our  events  favored  the 
first-motion  discriminant,  accuracy  of  first-motion  data  in  practice  is  so 
debatable  and  the  threshold  so  high  for  positive  reading  that  these  results 
are  of  limited  importance.  The  theory  of  higher-mode  surface  waves  was  also 
discussed.  However,  examination  of  long-period  recordings  for  our  earthquakes 
did  not  reveal  sufficient  positive  identifications  of  these  arrivals  to  util¬ 
ize  them  as  another  discrimination  parameter. 

A  stepwise  multiple-discrimination  program  provided  a  means  of  ranking 
the  importance  of  the  discriminants  in  classifying  the  two  groups  of  the  data 
set  (explosions  and  earthquakes)  used  in  this  study.  Ignoring  first  motion, 
they  were  ranked  in  the  network  experiment  in  order  of  value  as  discriminants 
thus : 

1  -  short-period  P  magnitude 

2  -  short-period  P  spectral  ratio 

3  -  LQ  magnitude 

4  -  complexity 

5  -  long-period  LR  spectral  ratio 

6  -  LR  magnitude 

7  -  short-period  S  magnitude 
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8  -  long-period  P  magnitude 

9  -  long-period  S  magnitude 

This  ranking  is  tenuous  because  of  the  large  number  of  missing  data  values, 
forcing  reliance  on  poorly  estimated  values  for  several  discriminants  for  many 
of  the  events.  Changes  in  the  ranking  would  undoubtedly  occur  witli  a  new 
"training  set."  Note  also  that  this  is  a  single-dimension  viewpoint  and  that 
the  worth  of  each  variable  is  highly  dependent  on  the  sampling  procedure. 

For  instance,  although  long-period  S  magnitudes  were  distributed  similarly  for 
the  earthquakes  and  explosions  of  this  study  (see  Figure  14),  it  was  shown  that 
in  conjunction  with  other  data,  long-period  S  waves  are  a  valuable  discriminant. 

Results  of  numerous  experiments  with  the  stepwise  discriminant  analysis 
program  have  been  summarized  adequately  in  the  test  and  tables.  Reiterated 
here 'are  these  important  results: 

1)  The  lowest  a  posteriori  probability  of  correct  classification  for 
the  23  events  was  .983  when  all  ten  variables  were  used.  Removing 
first  motion  decreased  this  to  .955. 

2)  A  particularly  good  discriminant  pair  was  versus  Mg  when  was 
estimated  from  the  combined  LQ  +  LR  amplitude.  The  lowest  a  poster¬ 
iori  probability  for  this  pair  was  .683,  the  value  for  the  Denver 
earthquake. 

3)  A  combination  of  short-period  discriminants,  namely  m^,  short-period 
spectral  ratio,  and  complexity,  misclassif ied  only  the  Borrego  Moun¬ 
tain  earthquake.  (Recall  that  m^  alone  gives  some  separation  between 
the  t  ;o  groups.)  Spectral  ratio  and  complexity  as  a  pair  misclassif ied 
four  events,  with  several  others  nearly  misclassif ied.  Thus,  discrimina¬ 
tion  using  short-period  ?-wave  data  alone  is  probably  not  adequate. 

4)  Stepwise  discrimination  analysis  using  the  data  from  only  one  station 
showed  significantly  smaller  average  a  posteriori  probabilities  than 
analysis  using  network  data  for  an  identical  parameter  set. 

5)  Discrimination  using  discriminant  functions  averaged  over  separate 
functions  computed  for  each  station  alone  is  inferior  to  starting  with 
network-averaged  parameters  and  using  a  network  discrimination  function. 

Corner-frequency  and  moment  parameters  were  not  tested  in  the  multiple-dis- 
discrimination  work  because  of  the  small  amount  of  data  accumulated  and  the 
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rather  subjective  way  of  estimating  parameters.  Still,  where  digital  data 
is  available  and  when  a  more  objective  approach  to  estimating  these  parameters 
is  implemented,  they  should  be  included  because  even  the  limited  results  here 
demonstrated  their  potential  classification  power. 

The  data  set  for  this  study  contained  only  single  explosions  and  no  shot 
arrays.  Traditional  discrimination  tests  involving  complexity  and  possibly 
those  involving  first  motion  would  be  inapplicable  to  multiple  explosions. 

Tests  based  upon  ratios  of  short-period  to  long-period  data,  such  as  Mg  -  m^, 
would  also  be  invalid  for  such  events.  Perhaps  the  best  teleseismic  discri¬ 
minants  applicable  to  shot  arrays  are  the  ratio  of  short-period  S  to  m^  and 
of  long-period  P  or  S  to  Hg;  at  regional  distances,  the  ratio  of  crustal  P  phases 
to  Lg  could  be  used  as  a  short-period  substitute  for  Mg  » 

There  are  no  reasons  to  think  that  path  or  receiver  effects  peculiar  to 
this  study  would  negate  the  conclusions  regarding  the  worth  of  the  various 
discriminants  studied.  Note  that  the  conclusions  are  robust  because  network 
estimates  of  parameters  were  made.  Less  positive  conclusions  would  have  been 
reached  if  only  data  from  one  station  were  used.  Statistics  of  the  multiple 
discrimination  experiment  clearly  demonstrated  this  test.  In  addition,  using 
a  network  obviously  lowers  the  threshold  for  application  of  the  discriminants 
and  allows  magnitude  estimates  to  be  made  according  to  Ringdal's  (1976)  method, 
which  is  a  powerful  tool  in  overcoming  magnitude  bias  due  to  a  combination  of 
detection  threshold  and  path  and  station  effects. 

Results  of  this  study  show  that  some  earthquakes  in  the  Southwestern 
United  States  are  only  marginally  separated  from  explosions,  even  with  the 
application  of  a  number  of  discrimination  parameters.  The  failure  of  the 
Mg  -  m^  pair  to  separate  events  clearly  is  an  important  finding  here.  No  com¬ 
mon  thread  seems  to  link  those  three  events  (Denver  earthquake,  Baja  earth¬ 
quake,  and  BENHAM  Aftershock)  that  were  particularly  hard  to  discriminant. 

They  were  in  distinct,  separate  areas,  had  no  common  source  mechanism,  and 
had  different  depths  of  focus. 

This  study  indicates  that  confident  seismic  discrimination  requires  both 
long-period  and  short-period  recordings  and  that  horizontal  motion  recordings 
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are  quite  helpful  in  both  passbands.  Though  pP  analysis  has  been  virtually 
ignored  because  of  the  generally  shallow  event  depths,  the  authors  believe 
that  efforts  to  extract  depth  information  could  lead  to  some  additional  aid 
for  classification  of  Southwestern  United  States  events. 
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APPENDIX  I 


Mean  Values  of  Seismic  Source  Radiation  Patterns 

In  order  to  evaluate  the  root  mean  square  average  of  the  P-wave  radia¬ 
tion  pattern  R  ,  the  square  of  this  function  is  integrated  in  a  spherical  co¬ 
ordinate  system  over  all  values  of  the  polar  angle  0  and  the  meridian  angle 
the  result  is  divided  by  the  area  of  a  unit  radius  spherical  surface,  and 
the  square  root  is  taken.  The  expression  for  the  P-wave  pattern  from  a  double 
couple  with  no  resultant  moment  is: 

R.  ,  =  sin  20*cos<(> 
u  9 

Thus  the  ras  value  is  given  by : 

~=  [/  2v  1 71  R.J-sineded*//  2v  1 17  sin0d0d<}>]  1/2 

O  0  t)  <p  o  o 


e<j> 


The  integrals  are  completed  with  the  result  that 

rTT  =  2//15 


The  Rayleigh-wave  radiation  pattern  is  a  function  of:  1)  0,  the  azimuth 
(measured  counterclockwise)  of  the  observer  with  respect  to  the  strike  of  the 
fault;  2)  A,  the  angle  between  the  slip  vector  and  the  horizontal,  measured 
counterclockwise  in  the  fault  plane;  3)  6,  the  dip  of  the  fault  plane  from  the 
horizontal;  4)  h,  the  depth  of  the  focus.  The  expression  for  the  modulus  of 
the  complex  radiation  pattern  is  taken  from  Harkrider  (1970)  as: 


XR(0,A,6,h)  =  [(Dq  +  D3  T3  +  T4)2  +  (Dx  Tx  +  D,,  T^2] 


1/2 


where 


Tl  =  sin0 
T2  =  cos0 
T3  =  sin  20 

T.  =  cos  20 

4 

Dq  =  1/2  sinA  sin  26  B(h) 
Dj  =  -sinA  cos  26  C(h) 

D2  *  cosA  sin6  A(h) 

=  cosA  sin6  A(h) 

*  -1/2  sinA  sin  26  A(h) 
Ah  =*  -U(h) 


AI-1 


B  (h)  =  - 


23  ‘ 


u(h)  +  (  ^  z(h) 

P«  /  h 


C(h) 


X(h) 


U(h)  =  radial  velocity  at  depth  h  normalized  by  W 


W(h)  =  vertical  velocity  at  depth  h  normalized  by  W 


Z(h)  =  vertical  stress  at  depth  h  normalized  by  (W  /CD) 

O  K 


X(h)  =  horizontal  stress  at  depth  h  normalized  by  (W  /CD) 

O  K 


“R 


=  W(0)  =  vertical  velocity  at  surface 
=  Rayleigh-wave  phase  velocity 


The  root  mean  square  value  of  xT D(9»X,6,h)  is  found,  for  a  given  focal  depth 
2 

h,  by  evaluating  x  D  f°r  all  possible  orientations  in  space  of  the  slip  and 
R 


dip  vectors  and  averaging  the  resulting  quantity  over  all  azimuths  0.  A 
spherical  coordinate  system  is  determined  by  allowing  6  to  take  on  values 
from  0  to  2rr  and  X  from  0  to  it  (see  Figure  Al).  Thus  we  have 


VW)b  -{tl’tl'f'o  X2r  sinA  dX  d6  it/fl 

>}  1/2 


2tt  j.  2tt  it 

o  o 


sinX  dX  d6  d0 


From  (A2),  we  have 


x  R(e,x,«.h) 


'V  +  D32t32  +  “dV  +  2DoD3T3  +  2DoD4T4  + 


2D3D4T3T4  +  »lV  +  2D1D2T1T2  "  D22t22> 


Performing  first  the  integral  over  0,  we  find  that  terms  of  the  form  / T 


equal  n  while  terms  of  form  /  T  T  or  /T  equal  zero,  so  that 
1  n  m  n 


/  ^x2^6 » X  ,6  ,h)d0  =  it  [2Dq2  +  Dx2  +  +  V*  +  D^] 


Integrating  (A3)  next  over  6,  we  find  that  all  the  terms  of  (A4)  are  non-zero, 
with  the  result  that 


/  2it/  2V  (0,X,6  ,h)d  6  d  0  =  TT2 

O  O  A  R 


sin2X  +  A2(h)cos2X  +  C2(h)  ]  (A5) 


AI-2 


Integrating  (A3)  finally  over  A,  we  obtain  from  (A5) 

2 

2ir/  2-u  f  it  2  (e  A,6,h)sinX  dX  d6  d8  =  tt2  [A2(h)  +  2B_(h)  +  2C2(h)] 
o  o  o  R  -5 

Using  the  fact  that 

.  j  2tt  y  71  gx  d6  d0  =  8tt 

o  o  o 

with  the  result  (A6),  equation  (A3)  reduces  to: 


=  4  A2(h)  ,  B2(h)  .  „2„  s 

XR(9,X,6)h  2  — 2 —  +  — 3 —  +  c 


The  corresponding  radiation  pattern  for  Love  waves  is  given  in  Harkrider 
(1970)  as 

.  ?  O  1/2 

|xL(e,X,6,h) I  =  [(Eo  +  E3T3  +  E4T4)  +  (EyTy  +  E^n 

where  EQ  =  0 

E^  =  cosX  cos6  G(h) 

E2  =  -sinX  cos26  G(h) 

E3  =  1/2  sinX  sin26  V(h) 

E4  =  cosX  sln<5  V (h) 

G(h)  = 

Uh 

V(h)  =  tangential  velocity  at  depth  h,  normalized  by  Vq 
Y(h)  =  tangential  stress  at  depth  h  normalized  by  (V  /C  ) 

O  L 

V  =  V(0)  =  tangential  velocity  at  surface 
o 

C  =  Love-wave  phase  velocity 

L 

Performing  the  integrals  as  for  the  Rayleigh-wave  pattern,  we  obtain  for  the 
Love-wave  pattern 

/  fxj  (e,X,6,h)d0  =  TT[E2  +  e2  +  e3  +  E2] 

Following  through  with  the  integration  as  for  the  Rayleigh-wave  pattern,  we 
obtain 


XL<e,A,6>h  -  i  ,=2(W  +^2 


V2(h) 


AI-4 


Some  special  cases  are  of 
the  surface  would  have 


practical 
A(0)  -  e0 
B(0)  =  (3 
C(0)  =  0 


interest. 


a 

o 


e 

o 


For  instance,  an 

G(0)  = 
V(0)  = 


earthquake 

0 

1.0 


at 


where  eq  is  the  ellipticity  of  the  particle  motion  at  the  surface, 
these  values  in  (A7)  and  (A9)  we  obtain 


Xl(6,A,6) 

h=o 


1 


I 


2/2 


Using 


2  2 

For  a  Poisson  medium,  a  =  33  ,  and 

XR(e-X'5>h-o-£ 

Also  of  interest  are  special  cases  such  as  pure  dip-slip,  purfe  strike-slip, 
and  thrust  faults,  for  which  the  expressions  simplify  for  all  depths.  The 
rms  radiation  pattern  for  these  cases  is  derived  by  changing  (A3)  to  inte¬ 
grals  over  0  only 

‘  “o’  (6-X-6>h>  d8//o*  d6,1/2 


The  integral  in  the  numerator  is  (A4)  while  the  denominator  is  rimply  2tt,  so 
that 


XR(9)h,M  L 


2D 


2  2  2  2 
+  D,  +  D.  4-  D_  +  D. 

1 _ 2 _ 3 _ 4 

2  J 


1/2 


Similarly,  using  (A8)  we  obtain  for  the  LQ  radiation  pattern 

1/2 


XL(0)h,A,6 


r  2  2  2  2  *1 

E1  +E2  +E3  +E4 

L  2  J 


Using  the  expressions  for  Dn  and  terms  given  above  with  the  appropriate  X 
and  6  values,  we  obtain: 


I.  Strike-slip  fault  (A  =  0°,  A  =  90°I 


TToT,  -  = 

K  77 

— for  .  li<±oJ 

V'Vi.s  /-j 

II.  Dip-slip  faults  (A  =  270°,  6  =  90°) 

V^.m  -  IxM/Ii! 


xL<»)hlij<  -  IIMIIt) 

n 

III.  Thrust  faults  (a  =  90°,  6  =  45°) 

V^.s-T  JB<h)2+I  4<h)2 

tigt  .te) 

L  h,X’6  2  /~2 

In  order  to  calculate  the  complete  theoretical  excitation  versus  fre 
quency,  these  results  mu:  t  be  multiplied  by  the  other  frequency-dependent 
and  medium-dependent  factors  given  in  Harkrider  (1970)  for  double-couple 
sources,  thus 


R(w)  =  kRX/  •  AR(u>)xR(w,e). 
L(w)  =  kl/2  *  Aj^uOXlGo.0)* 


where 


k  ,  k  -  wavenumber 
K  L 

An,  At  -  medium  amplitude  response. 


APPENDIX  II 


Derivation  of  the  Explosion  Spectrum 


a 


APPENDIX  II 

Derivation  of  the  Explosion  Spectrum 

From  Rodean  (1971) ,  the  relation  between  the  spectrum  of  the  reduced  dis- 

A  A 

placement  potential  4>(to)  and  spectrum  of  the  pressure  P (to)  at  the  cavity 
radius  (or  equivalent  elastic  radius  a)  is 

f(w)  =  ~2 - — - 2  P(w) 

to  +  2i£to  uj  -  a) 

,  o  o 

where 

p  =  density 

to  =  20/a 
o 

0  =  shear  velocity 
a  =  compressional  velocity 
5  =  e/a 


Also,  from  Rodean  (1971)  the  well-known  relation  between  the  displacement 

A  A 

spectrum  u(r,to)  and  V  (to)  is 


u(r,u>)  =  ~  (~+  — "j)  T(u) 
r 

where  r  is  the  distance  from  the  center  of  the  cavity.  Combining  (A-10)  and 
(A-ll),  we  have 


u(r,ui)  =  ~  (“  + 

r 


a/p 


2  2 
to  +  2i£to  a)  -  to 
o  o 


P(to) 


Equation  (A-12)  can  be  evaluated  at  r  =  a  in  the  near  field  and  at  r  in  the 

2 

far  field  where  we  neglect  the  1/r  term,  so  that 


G(a,U)  =  -  (^  +  V  2"  l/P - 2  P(o)) 

a  to  +  2i£u>  to  -  to 
o  o 

and 


u(r,w)  z 


-  (,—) 
ar 


a/p 


2  2 
<o  +  2i to  -  to 
o  o 


P(u>). 


The  ratio  of  the  near-field  and  far-field  displacement  spectra  is  then  simply 


u(r,to) 

* 

u(a,to) 


lUL  /  +  1_) 

ar  '  (aa  l’  * 
a 


All— 1 


For  a  step  function  at  the  cavity  radius  with  final  displacement  Do>  u(a,u>) 
D  /iw  so  Chat 

°  D  . 

,  .  «,  o  ,  ,iw  ,  1  v 

u(r’w)  „  —  /  (—  +  -2  ) 

a 

The  modulus  of  this  expression  is 


-  aD  „  2  -1/2 

|u(r,oj)  |  :  (w  + 

a 


The  same  result  is  also  obtained  by  initially  solving  the  equation  of  motion 
with  a  boundary  condition  of  a  step  displacement  at  the  cavity  radius. 
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